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Uber die Excistenz kanonischer Transformationen 
bei mehreren unabbdngigen Veranderlichen 


HELMUT RUSSMANN 


Vorgelegt von J. L. ERICKSEN 


1. Einleitung 
Entsprechend der Herleitung der Hamiltonschen Differentialgleichungen aus 
Variationsproblemen in einer unabhangigen Verdanderlichen gelangt man von 
Variationsproblemen in mehreren unabhangigen Verdanderlichen zu Systemen 


partieller Differentialgleichungen der Gestalt ‘ 
(1) 2 Paix ot ide Vou, x; =H, 


(aT re ING AD yt) 


wobei H=H(q,,.-., Gy, Pir» --+» Prn> %1> +++» %) eine in einem gegebenen Gebiet A 
zweimal stetig differenzierbare Funktion ihrer »+Nu+N Argumente ist. Die 
Indizes hinter dem Komma bedeuten dabei wie iiblich partielle Differentiation. 
Analog zum klassischen Fall »=1 existiert auch fiir ~ >1 die Gruppe der 
,Kanonischen“ Transformationen, welche jedes System (1) wieder in ein solches 
iiberfiihren. Die notwendige und hinreichende Bedingung dafiir, daB eine einmal 


stetig differenzierbare und umkehrbare Transformation - 
Yaa =JaQrs ++ Ovi Piss ee) Pyrns ay ++ Mn), 
(2) Pai = Pai (Qrs +++» Ovi Pras) Pym My +++) Hn) 


(ogee Oe a St Anes 500) 
kanonisch ist, hat D. G. B. EDELEN angegeben?. 
Zur Formulierung dieser Bedingung bezeichne 
M+; M,., eeey My 44 
M=(M,,) =| : 
M41 5 ae Meine 


die Funktionalmatrix der Transformation (2), wobei also 
M,, = (Ja.,05) , My 41= (Ju, Ps) , 
M11 =F (Pxi,0,) , M1441 = (Pai, Pa 4) 


1 EDELEN, D. G. B.: The Invariance Group for Hamiltonian Systems of Partial 
Differential Equations I, Analysis. Arch. Rational Mech. Anal. 5 (1960). 
Arch. Rational Mech. AnaJ., Vol. 8 25 
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Beweis. Aus den Gleichungen (5) erhalten wir durch spezielle Wahl der 
Indizes fiir 7 >1, k >1 


(12) M,;M;,—™Mj,;M.= 9, 

, ’ 1 
(13) M,,M;,— fi, My, =~ OnE, 
(14) Mi1.M,,—MiM,=9. 


Nach Lemma 1 ist det M,,-+-0. Daher folgt aus (13) 
Me Mia + Mi;M;,Mz3 


1) 
und aus (43) und (14) fiir7+-k 
M,,=M,,MiiM,,. 
Setzen wir die beiden Ausdriicke fiir Mj; und M;, in (12) ein, so ergibt sich 
“Mit M,,= 9, also 
(15) M,,—O ESA); 
are.d: 


Lemma 1 und Lemma 2 zusammen bilden gerade die Aussage des in der 
Einleitung angekiindigten Satzes. 


3. Die explizite Darstellung der L6sungen von (3) 
und der kanonischen Transformationen 
Um sdmtliche Lésungen des Systems (5) zu bestimmen, beachten wir zu- 
nachst (13). Wegen (15) ergibt sich 


(16) Misses = in Maa" (6, b= 4,2, a», 2): 
Setzen wir noch 
(17) M,4,—M,7°S; (Ree noe 


so folgt aus (14) S;=S,. 

Die Gleichungen (15), (16) und (17) mit symmetrischen Matrizen S; sind 
also notwendig fiir das Bestehen der Gleichungen (5). DaB sie auch hinreichen, 
verifiziert man durch Einsetzen. 


Die allgemeinste Lésung des Gleichungssystems (3) lautet also fiir  >1 


Mi, O one Oo 


we. 1 hex 
(18) M— Myi"S, -M;;* O 


re 4 De 
Mi;'S,, O 4 Ma 


mit beliebiger nicht ausgearteter Matrix M,, und beliebigen symmetrischen 
Matrizen S;. 
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Zum Schlu8 wollen wir noch alle kanonischen Transformationen fiir n >1 
bestimmen. Seien hierzu 


(19) Gi = Gal Ors Hoch Onetis ta, | (cee Beet cst) 


in einem gewissen Gebiet ihrer N+ Argumente stetig differenzierbare Funk- 


tionen mit det (¢,, 9 A) =-0. Wir bezeichnen die Spaltenvektoren mit den Kompo- 
nenten F;,...,.,; mit P, und erhalten aus (18) 


(20) p= Mii" (7 B+ a) (¢=1,2,...,n), 


wobei die Komponenten a,; der Spaltenvektoren a; stetig differenzierbare Funk- 
tionen in Q),..., Qy, %1,---, %, sind. Differentiation von (20) nach Q, und 
Multiplikation mit Mj, ergibt dann 


, , foe 4 
mM), Pi, Os me M,, My, 105 P+ a) a 4,08 


= Mh1,0, b; + Gi, Og, 


wobei b;= — (TP +a) ist. Hieraus folgt 


S;= My, Mj... = (= By iI, Oop 04) (4.5, Op): 


Die erste Matrix in der Summe ist symmetrisch. Die Symmetrie der zweiten 
ist gleichbedeutend mit der Existenz von in den Q, zweimal und in den x; einmal 
stetig differenzierbaren Funktionen c,;, so daB 


(21) Ce —- C;, OP, 


gilt. Durch (19), (20) und (21) werden also alle kanonischen Transformationen 
fiir » >1 dargestellt. Diese Formeln finden sich auch bei EDELEN‘, doch geht 
aus der dort gegebenen Ableitung nicht hervor, da8 keine weiteren kanonischen 
Transformationen existieren. 


4 Siehe1 S. 154. 
Mathematisches Institut 
der Freien Universitat 
Berlin-Dahlem 


(Eingegangen am 30. Juni 1961) 


Surface Waves in Deformed Elastic Materials 


M. Haves & R.S. RIVLIN 


1. Introduction 


In the present paper we apply the theory [/]* of the superposition of in- 
finitesimal deformations on finite deformations in an isotropic elastic material 
to the study of the propagation of surface waves in a semi-infinite body which 
is subjected to a static, pure homogeneous deformation. 

Waves of the type discovered by RAYLEIGH [2] are first examined. It is 
seen that the motion does not differ essentially from that formulated by RAYLEIGH 
on the basis of classical elasticity theory. There is no dispersion, as might be 
expected, since no fundamental length is involved either in the dimensions of 
the body or in the deformations it undergoes. Also, the motion is two-dimensional 
in the sense that there is no displacement transverse to the direction of propa- 
gation of the waves. A form of strain-energy function, equivalent to that used 
by MuRNAGHAN [3], is employed in obtaining the correction, due to an initial 
small pure homogeneous deformation, in the frequency of propagation of waves 
of this type. The method, outlined by RIvLIN [4], of deducing the results for 
an incompressible body in second-order elasticity theory, from the corresponding 
second-order results for a compressible one, is used to obtain in explicit form 
the displacement components for the case of an incompressible body. 

The possibility of having surface waves with transverse horizontal movement 
in a semi-infinite body, subjected to a pure homogeneous deformation, is next 
discussed. We find that such waves may not be propagated in a principal 
direction unless stringent conditions are imposed on the strain-energy function. 
Then, following Love [4], we assume that there rests on the semi-infinite body 
a layer of material with properties differing from those of the semi-infinite body 
and subjected to a different pure homogeneous deformation. The principal 
directions are taken to be the same in the layer as in the subjacent body. An 
infinitesimal deformation is superimposed on the composite body and conditions 


for the existence of Love-type waves are sought. The conditions obtained are 
similar to those found by Love. 


2. The Equations of Motion and Boundary Conditions 


We take the origin of a rectangular cartesian coordinate system (X;) on the 
surface of a semi-infinite body which occupies the region X,>0. This body is 
subjected to a pure homogeneous deformation defined Dy = 


Ny =A, Xy. (2.1) 


* Numbers in square brackets refer to the references at the end of the paper. 
** The summation convention will not apply when Greek subscripts are employed. 
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The stresses o;; associated with this deformation are given [6] by 

O55 = 0 (a = B) , 

Ong = 215 {AW + Ag (T — Ab) We, + Ty WG}, (2.2) 
h=h+h+4%, L=HR+BR+R, L=HRR, W,=eW/al,, 


W being the strain-energy function. 


Finally, we superimpose on the deformation (2.1) an infinitesimal deformation 
such that the particle whose position is x; is displaced to a point %;, where 


X, = %;, + €u;. (2.3) 


Here ¢ is a constant such that its squares and higher powers may be neglected 
in comparison with unity. The equations of motion are found to be [7] 


00;; = 2 
ax,  § ap” (2.4) 
where 6;,, the incremental stress associated with (2.3), is given by 
Gag = 215 *{W, + (Ty — A — A) Wo} (An Mae + AB Ma.) («+ 8), 
(2.5) 


3 
xe = D3 Cap UB B> 
p=1 


and g is the density of the material in the state of pure homogeneous deformation. 
Differentiation with respect to x, is denoted by |, and cy, is given by 


Cup = 215 *{[W, + (Ly — Aa — 45) We] Az(2 6,8 — 1) + Ag Ag We + 5 Wy + 
+ 20a [Maa + (Fy — Ag) Wer) + 1g War Ap + 2a [Me + (Lr — Ag) Wosl+ (2.6) 
+ I Wao} Ap (Ly — Ap) + 2{0a [Ms + (Li — 2g) Was] +15 Mea} Lo}, 


where @?W/0I,,0I, is denoted by Wy ,. 
The surface tractions 7;, in the state of deformation (2.1), are given by 
T;=0,,1, (2:7) 


a jG? 


where l; are the direction-cosines of the normal to the surface in that state. 
The surface tractions, 7;+e7;, in the state of deformation (2.3) are given by 


T, + ¢T, = (o;; + €G;;) (+ el), (2.8) 


where o;;+€6,;, J+ el; are the stress components and direction-cosines of the 
normal to the surface in that state. It has been shown [6] that 


i; == zh; Lin L, (Um, n = Un, m) re L, Up 5 : (2.9) 
Neglecting powers of ¢ higher than the first and using (2.7), (2.8) reduces to 
T=6,;1;, +6,,1;. (2.10) 
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3. The Frequency Equation 
We assume the surface x,—0 to be disturbed in some manner and suppose 
the displacement components u; to be given by 
uy = f,(%) exp (pm +4972), 
Uy = fa (%2) exp e(p%y+ 2), (3.1) 
Uz =0, 
where f, and f, are to be determined. Substituting these values in (2.5) and 
using (2.4), we find that 
{(0 g? — G3 B®) + AOD fh + (42+ A) pD i, =0, 
{(@?— Ay bp?) + C9 D*} fo + (cor + 450) epDA =O, 
where 
Diddy, b=21, 4, 1, — #— Aw). 
The third equation of motion is satisfied identically. To solve equations (3.2), 
we write 
h=— (q2+4,6) epDV, 
f= (09? — 41 0? + 430 D)V, 


where V(x.) is to be determined. Evidently, these values of f,, f, satisfy the 
first of equations (3.2) whatever may be the value of V. Substituting them in 
the second of (3.2) we find 


{03 b cy. D* + [co2(@ q? — Cy, p?) + A b(og? — Ai bp) a 


(3-3) 


; (3.4) 
+ (Cra + Az) (Cor + Ag 2) £2] D? + (0.9? — 11 B*) (0g? — Aid P»)} V =O. 
The general solution of (3.4) is 
V=Aerm* + Bent % ot Ae ems 4 B; ema *s (3.5) 


where A, B, A,, B, are constants and -+-m,, +m, are the roots of the equation 
Az 5 Cog m4 te {Co2(0 9? — C11 p?) + A3 b (0 gq? — Ay bp?) + 
+ (Cre + AZ) (Cor + Aj b) p?} m? + (O02 "G4P) (oy Ai b p?) =". 


We may now find ,, u, from (3.1), (3.3) and (3.5). If «,, w, are to tend to 
zero aS X» tends to infinity, then, taking m,,m,>0, we must have 


Ay = Dyes On 


(3.6) 


Then 
My = (Cro + ALD) op (m, A e~™* + my Be-™*) exp u(px, +), 
Uy = {(0 9? — cy p2 + mi AZ b) A e~™* (3-7) 
+ (99? — G1 f? + m3 Az b) Be-™*} exp u(p%, + 92). 
These displacement components correspond to waves, propagated parallel to the 


surface of the body, with their amplitudes decreasing with distance from the 


surface. They are analogous to Rayleigh surface waves in the classical theory 
(see, for example, [2]). 
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We shall now specify the boundary conditions to be satisfied on the surface 


_ of the body and from them determine an equation for the ratio Plog. The 


values of J; on the surface may be found from (2.9) by replacing J,,/,, 1; by 


_ 0, —1,0 and taking note of the fact that ~,, , are not functions of x,. We 


find eae 1,=l,=0. Thus, if we take the surface to be free of traction both 
in the state of pure homogeneous deformation and in the state of the super- 
imposed infinitesimal deformation, the boundary conditions 


O22=0, 031 M21 — 2=0, O22 =0, 032 =0, (3.8) 


must be satisfied on the surface x,=0. The last of these conditions is satisfied 
identically. The first condition gives 


Wy + (Ly — 4a) We + Ap Ag = 0. (3.9) 
The second condition gives, using (2.2) and (2.5), 
{aw + AL, — Mi) W, + Is W3} Uys 


(3.10) 
{W, + (L, — At — 23) We} (Aj e+ Ag Mr, 2) = 0- 
With (3.9), this may be written 
Ay Ag (We + 23 Ws) (ita, 2 + Ma,1) = 0. (3.11) 
Thus, provided AA} (W, +23 W,) +0, we have 
Uy,2 + Ue 1 =0. (3.12) 
The third boundary condition (3.8) gives, using (2.5), 
C21 Uy,1 + C22 U2 9 = 0. (3.13) 
Substituting u,, u, from (3.7) into (3.12), we find 
PIf(o 9? — ¢11 £2) + mi{(Ag — Aj) b — Gy a}} A + (3.14) 
+ {(0 9? — cy B®) + me {(Ag — Aj) b — ey 5}} B] = 0. 
Thus we have either bs (3.15) 
“G {(9. 9? — 11 B?) + my [(Ag — Az) 8 — Gel} A + (3.16) 
+ {(@ 9? — ey £?) + mg [(Ag — Az) b — ey 2]} B= 0. 
From (3.13), again using (3.7), we find 
my [ (Cre Cox — 611 C22 + C91 AZO) £? + Co 2 0G + Coop b my] A + (3.17) 


+ mg [(Cr2 C1 — &11 622 + C21 Ai b) B® + Co: 09? + Coe Az b ms] B=0. 


In order that (3.16) and (3.17) have non-vanishing solutions for the integration 
constants A, B we must have 


(09? — 412?) + ma {(As — Az) b — eo}, 
My {(Cr2 C21 — 611 C22 + Co1 Aj b) p® + Con 0g? + Cae Ag b my} ry 
(09% — G1?) + ms { (Ap = Aa) Be Cra}, 
My { (C2 C21 — C11 S22 + Cor Aj b) p® + C2209? + Coo Ag b m3} 


0. (3.18) 
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This equation may be written 


(Mm, — Me) {(@ gq? — (1 *) [(Cy2Ca1 — C11 C22 + Cor A b) £2 + Co20 97] + 
+ Cg 5 b(o gq? — Cy, p?) (my sis ms) — Co9 Ap bl ere i (Ay as d3) b] my; Ms a (3.19) 
— My My (Cy + Ay b) [ (Cre Co1 — C11 S22 + a1 (At == 43) 5) p* Co20q"} =), 


Inserting the values of m?+m}, m{mz obtainable from (3.6), in (3.19), we find 


either 
M, = Mz, (3.20) 


or 
(Cry + A? b) (0 g? — 1 B®) {0 9? — (At — AG) b B?} 3.24) 
= (44. + 44d) {[¢11 C22 — C9 Cy — C94 (Ai — 43) b] p? — Co20 q*} My, Mz - 


Assuming ¢,,+Aj{b-++0, squaring both sides of (3.21) and again using (3.6), we 
have 


(0.9? — cy 8°)? {0g — (Ay — Ag) OP? Ag b Cao 


. - ; s 4 (3.22) 
= (0 g?—¢c, P*) (@ g—At b p?) {es C99—Cy 9 Coy—C91 (Aq —Ag) 0] P? Con Og } : 
Thus either 
eg=%, Pp, (3.23) 
or 
(0.9? — e11 ) {0 42 — (A — 28) b 8} 28 cag nae 
=(0¢?— rr b p?) {Les Co — C129 C94 — Cai (AG = 43) b] p? — Ces 0 Gy? 
Writing 
Pleg?=«a, (3.25) 
(3.24) may be written 
(1 — ey «) [1 — (Aj — Ap) ba]? Ap 5 co, (3.26) 


eet ii ba) {[¢11 Coe — C19 Cy — Cay (Ay — AQ) b] a — Gaal 


In the next section we shall show that the solutions given by (3.15), (3.20), 
(3.23) are degenerate. We refer to (3.26) as the frequency equation since it governs 
the frequency with which Rayleigh surface waves of given wave-length are 
propagated. We note that, since the coefficients of « in (3.26) do not involve 
either p or g, there is no dispersion, as indeed is expected since there is no funda- 
mental length involved either in the dimensions of the body or in the deformation 
it undergoes. 


4. The degenerate solutions 


In this section we examine the solutions (3.15), (3.20), (3.23) and the case 
when ¢,.-+Azb=0. 


If p=0 (cf. (3.15)) then the equations (3.2) become 


(99? + Ab D*) , =0, 


(0.97 + Cy. D?) jo =0. (4.1) 


| 


| 
, 
) 
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Solving for /, and f/, and substituting the expressions obtained in (3.1), we find 
for u, and w,, 

tt, = [A cos (9 42/22)! xy + B sin (9 g2/420)! x9] exp et, Hes 

My = [Ay 60S (9 g*/cg)* vq + By sin (0 g?/cy9)4 x9] exp eg, | 


where A, B, A, and B, are constants. In order that the boundary condition 
(3.12) be satisfied we must have either g=0, or B=0. In the former case we 
do not have a wave. We conclude that B=0. In similar fashion, from the 


- condition (3.13) we conclude that B,=0. Thus we have 


u, = A cos(o g2/A2 b)? x, expugt, 
1 (0 9?/Az 5)? %2 exp vg (43) 


Us, = A, COS (0 g?/Cy5)? x, expigt, 


where A and A, are arbitrary constants. The expression for u, represents a 
standing wave in the x-direction, polarized in the x,-direction, with wave-length 
22 (A3b/oq?)? and angular frequency g. The expression for w, represents a standing 
wave in the x-direction, polarized in the x-direction, with wave-length 2 (cy5/0q?)! 
and angular frequency gq. 


Again, if m= my (cf. (3.20)), from (3.6) we have 


425 b Con (0g? — 1 P*) (0 — Ai b p?) 


(4.4) 
= {co (0 G@— ¢, p*)+- Ap b(og? — Ai b p?) + (42+ Ai b) (Cx + ag b) poh 


In this case the general solution for V is 
V =(A+ 8x,) e~”* + (A, + B, x2) e”*, (4.5) 


where A, B, A, and B, are constants and m is the positive root of (3.6). We 
may find the displacement components “,, “, from (3.1) and (3.3). In order 
that u, and u, tend to zero as x, tends to infinity we must have 4;=B,=0. 
Then 


Uy = (449 + Ab) ep{m(A + Bx) — Bhe~™™* exp u(p +92), 


(4.6) 
i = {(0 q?— C fp? + 13 bm?) (A + B x.) — 2436 m BY e—™* exp Up x,+qt). 


If we substitute these values of u, and uw, in the boundary conditions (3.12) 
and (3.13), we find two equations replacing (3.16) and (3.17). In order that 
these equations have non-trivial solutions for the constants A and B, we must 


have 
Pte Ge CaP) a [(As ii) b — cy,] m*}, 2p m[(At -: 3) b+ C9] 
aa M{Co1 (Cre +45 b) 2? + C.2(0— ah? + Ag bm*)}, =0. (4.7) 


{C21 (Cre + Ai b) £? + Con(o g? — C1 2? + 3.136 m*)} 


If we expand this determinant and, bearing in mind the fact that m,=mg, 
substitute the values of m? and m* from (3.6), we find an equation for #? in 
terms of g?. However, the resulting equation is different from (4.4). Accordingly 
we conclude that we do not have a solution if m,=mzy. 
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If og?=Cy,P? (cf. (3.23)), then the equations (3.2) may be written 


D{2RbD fy + (Cre + 18) up f.} =0, 


(4.8) 
{(414 — Aj b) p? + Coo D*} fy + (coy + 43 8) tp Df, =0. 


Solving these equations for f, and f, and using the fact that the displacements 
must tend to zero as x, tends to infinity, we find 


Uy = (Co + ALD) ep (Ag b m,) 1A e~-™* exp L(p % + 98), 


(4.9) 
Ug = Ae ™* exp t(p 4+ 94), 


where A is a constant and m, is the positive root of the equation 
Ag b Cog m® + {(G1 — Ai b) Ag b + (C2 + Ai b) (Cor + 22 b)} p?=0. (4.10) 


In order that the boundary condition (3.12) be satisfied we must have 


P (cre + Atb — Ag b) A =0. (4.11) 
Thus, either 
A=0; (4.12) 
Or 
P=; (4.13) 
or 
bya A? b= 2b. (4.14) 


If (4.12) is valid, then, from (4.9), w4,=u,=0. If (4.13) is valid, then, from the 
relation 0q7?=c,,p?, we have g=0. Hence, from (4.9) 4,=0, u.=Ae ™™. The 
boundary condition (3.13) then yields 


and equation (4.10) reduces to 


22d co, m? = 0. (4.16) 


Both (4.15) and (4.16) are satisfied if #0, in which case wz, is a constant, or 
if C3 =0 when we have 


= 0 Mo A eWee, (4.17) 


where A is arbitrary and m is positive. This is a trivial static solution. 


If (4.14) is valid we find c,,;=0, from (2.6). To satisfy the boundary condition 
(3.13) we again arrive at (4.15). Then, from (4.10), we find either p=0, or 
(cy, — ALD + 12d) 22b=0. The case P=0 has already been discussed. We now 
have to consider the two subcases (i) c,;= (Ai — 23) b and (ii) 225=0. In subcase 
(i) we have ¢,;= —cy,= (Aj — 43) b. Then (4.9) gives a solution, provided 72 -£0. 
In subcase (ii) we have #,=0, u,=Ae~™*expu(p x,+ qt). The boundary con- 
dition (3.13) is not satisfied unless m,A=0. This leads to Uy =U,=0. 

Thus, we conclude that if oq?=c,,f2 and also if c,= — Cy9= (Ai — /3)d the 
solution given by (4.9) is valid provided 72-0. In this case A is arbitrary 
and m, is an arbitrary positive number. We note that (4.9) gives a solution 
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with the property that there is no finite depth at which the components uw, and 
uy Of the displacement vanish. 


Finally, we discuss the case when c,,+#{/=0. From (2.6) and the definition 
of b, it follows that 


Cd Cen Oe (4.18) 
Hence, in this case, the equations of motion (3.2) reduce to 
{(o@— 110%) + BODY =0, {(0@—2bf%) +0,D%f,=0. (4.19) 
The solutions of (4.19) are 
f= Ae Ae | fe = Been B, ee, (4.20) 
where m, is the positive root of 
156 mi = 4, f?— 09, (4.21) 
and mz, is the positive root of 
Coo My = A bf? — 0 @?. (4.22) 


A, B, A, and B, are constants. Since the displacement components must tend 
to zero as x, tends to infinity, we obtain from (3.1) and (4.20), 


wy =Ae ™*expi(px%,+gqt), w=Be™*expi(px%tgt). (4.23) 
The boundary conditions (3.12) and (3.13) are satisfied, provided 
WA — Up B= Oye bP Cy, A Co5 Mab =O). (4.24) 
Hence, eliminating A and B, we obtain 
Cag tly Ne == 054 P (4.25) 
Using (3.25), (4.21) and (4.22), we obtain from (4.25), 
chy (cy. % — 1)* (At ba — 1)§ = — cy, (AR D)P a. (4.26) 


We note that /3d/o is the square of a velocity with which a plane wave may 
be propagated along the x,-axis in an infinite solid of the material when it is 
subjected to the same deformation (2.1). Hence, since c,,+/3b=0, we have 
Cy,<0. Also, ¢,/@ and Ajb/@ are each the square of a velocity of propagation 
of a plane wave along the x-direction in an infinite body of the material when 
it is subjected to the deformation (2.1). Accordingly, we conclude from (4.26) 
that if Rayleigh-type surface waves can be transmitted in the x,-direction when 
C19 +Aj{b=0, then the velocity of propagation of these waves must be less than 
two of the velocities of propagation of plane waves along the x-direction when 
the body is infinite. This result is also true when Cg +A?b+0. This follows 
from the fact that the equation (3.6) must give positive values for both mj and 
mz, for otherwise the displacement components (3.7) would not tend to zero as 
x, tends to infinity. Since c,, and A3b are positive, we conclude from (3.6) that 
necessary conditions for m} and mj to be positive are 9g?<c,,p 2and og?< Aj bp. 
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5. Second-Order Theory 
In this section we derive the changes in the velocities of propagation of 
Rayleigh surface waves, which are produced by an initial small pure homogeneous 
deformation of the body, by using the second-order elasticity theory. This is 


done in the following manner. 

The form 
W=ahtaht+ahht+uh +4sh (5.1) 
is assumed for the strain-energy function, where 


Ah=h—3, peal is. eleven 1; (5.2) 


Then, evaluating W,, W,, and taking 
Al ep (5-3) 


we use (2.6) and the expression for 6 following (3.2), to obtain expressions for 
C11» Coa» C19» Ca1, Arb, AZb in which we neglect terms of order higher than the first 
in the e,’s. The coefficients in (3.26) may then be found. Inserting them, we 
find the first approximation to the classical frequency equation. 


From (5.1) we obtain 


ow aW ow aw 2 
| — | | 
he Mp ake ads As — 20, + 2(a,— 43) f+ 43 ),+ 3a. Jz, 
ow = OW SOW 
aT aT hoy i eee 
ow ooW 
oi eae 9 ees (5.4) 
aw 
ar? = 2 (a, — 2a3) + 6a, jy, 
ne, AWS. CVE ML Ue 
él, el, a” OL, Ol, (Olgele ela Ok ; 


Using (2.2) and (5.3) we find, on neglecting terms of order higher than the 
first in the e,’s, that 


L,=3+2(4+ e+ 63), Iz=3+4(+e+6s), Tg= 1+ 2(@ + € + é 3), (5.5) 
and hence, from (5.2), that 

A=2(4+¢-¢), ~=0, “f=0. (5.6) 
Inserting these values of J, in (5.4) we obtain 


W, = as — 2a, + 4 (ag — as) (€, + € + 6), 

Wy = a, — a5 + 2a (& + e+ 65), W,=4;, 

Wy = 2 (a, 2a3) + 12a, (e + e+ 6), Wi. = 43, 
Ws =W,3=W.2 = =0. 
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Bearing in mind the last relations of (5.7), ¢y1,Co9, C19, Co, defined by (2.6), 


are given by 
| 


q 


) G1 = 24 Is*{W 4+ (— A) WRB + 20W +420 (L Wor}, 
Coo = 2g 15 *{W, + ee Wat 4h ite ai} 

bya = 2A TS *{— Wy + (Ag — Ag) We + AB AR + 208 Wy + 208 (1, + 23) Way}, 

Da, = 20 Tg H{— Wi + (2 — 28) Wy + AE ABW + 282 + 244 (T+ 28) Was}. 


(5.8) 


_ Expanding the terms of (5.8), using (5.3) and (5.7), and neglecting terms of 
_ order higher than the first in the e,'s, we find 


Cy = 4{2a, + 4 (2a, + 3.a4) & + (4 + 445 + 1244) (C2 + &s)}; 


. Coe = 4{2a,+ (a, + 4a3 + 1244) (& + 3) + 4 (2a, + 3 ay) ey}, 
Crp = 4 {ay + 2ay + (a + 40g + 1204) & +4 (ay + 3.44) 2 + 

| 
| 


(5.9) 
+ 2(— a — 2a, + 343 + 6a, + as) és}, 
Coy = 4{d, + 2a, + 4 (ag + 3.4) & + (a + 443 + 1204) €y + 
+ 2(— a, — 2a, + 3.3 + 644+ Gs) és}. 
_ From the definition of }, viz 
| b= 215 4{W, + (1, — A = Aa) MW}, 
and (5.7), we find 
Aib = 2{— a+ (— a, + 4a, — 245) e+ (a + 4g — 245) &2 + 
| (3a, + 4a, — 2a, — 245) es}, (5.10) 
Ap bb = 2{— a + (a, + 44 — 2ag) & + (— a + 4a, — 245) & + 
| (3a, + 4a, — 2a3 — 24s) é3}. 
We now evaluate the coefficients in the frequency equation (3.26), viz 
(1 — cy, 0) {4 — (AZ — Ab) Da} AED cop Bk 
= (1 — Ai b x) {[e11 Coe eee LoCo am C1 (Aq a 23) ) b] a — Coo}. 
We note from (5.10) that 
(Ay — ag) b = — 4a, (6, — &). (5.12) 


We may therefore neglect terms involving (Ai — A3)2b2 when expanding (5.11). 

We find 

Web 4c, 5 Cag — C2 n1 — Cn (Ay — Ag) 5}? 8 {Tey 4 Cop — C49 Coy — Coy (At — AG) B]2 + 
+ 222 D coo [11 Co2 — C12 Cos — Cas (At — Aa) b] — 203 b cy3 Coa (Ai — Ag) BF 2+ 
+ {2Cy2[¢11 C22 — S12 C21 — Cay (Aq — Ag) 6] + Ai b co. — 
— 22 b cog [¢r1 + 2(Ai — Aa) b}} o — Coe + Cop Ag bd =0. 


(243) 
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Using (5.9), (5.10) and (5.12), we find the coefficients in (5.13). That of «3 is 
51243 (a + 445) {— a (a + 44) + 

+ (— a2 + 48.43 + 4a, dy — 18a, dg — 48.4, 4g — 244 dg) &y + 

+ (— 3a? + 48.03 + 44, ay — 18a, a, — 48.4, ag — 244, As) €g + (5.14) 

(14 a + 48.45 + 56a, ay — 26.4, dg — 48 ay ag — 


— 104, ds — 24a, ag — 24a» Gs) es}; 
that of «? is 
— 256 {a} (a, + 442) (a + 642) — 
— [4 a, ay (48.05 + 14.4, dg — 29 a, dy — 720, Ay — 2409 a5) + 
+ a? (— at — 2a, ay — 204, a, — 604, a4) ] ey — 
— [4a, ay (48a — 2a, dy — 25 a, a3 — 724, Ay — 244g G3) + (5.15) 
+ 2a}(— 2at — 15 a, a, — 8a, a, — 304, a4) ] ey — 
— [4ay ay (4843+ 664, ay — 39 ay dg — 72.0, 4 — 15 a 4; — 24g Ag — 245 5) + 
+ ai (7 a} + 82a, ag — 28.4, ag — 60a, a, — 84, a5) | 5}; 
that of « is 
64{— 44, a3 (a, + 445) — 
= 24 — 6402 P41 ai as 8a.a, 1 2440, a, — 126, a 
+ 68a, ay a3 +192 a, Ay a4) ey — 
— (— 6443 + 21a} a, + 24a} a,+ 76a, a3 + 32a3 a, + (5.16) 
+ 28 a, ay a3 + 192, ay Ay) 5 
(2a) — 6443 — 8% a, + Sai a, + 244) a, — 80a, a2 + 324.4, 
+ 3203 as + 80a, dy dg + 1924, dy ay + 16.0, ay G5) es}, 


and that of «® is 
— 8 {2a (a, + 4ay) + 
+ (af — 845 + 64, a, + 4a, ag + 124, ay + 36a dg + 96 ay a4) &y + 
+ (56.43 + 10, ay + 12a, dy + 4a, ag + 96a, ay) €y + 
+ (aj — 805 + 2a, ag + 4a, ay + 124, ay + 36.05 ag + 96.5 Ag + 4a ds) eye 


(5.17) 


In solving (5.13) we must bear in mind equation (3.9), viz 
W, + (L, — 42) Wa + Ai As Wy = 0, (5.48) 
which, with the aid of (5.7), yields 
es = — & — {2a9/(a, + 2a) ey. (5.19) 


We eliminate e, from the expressions (5.14), (5.15), (5.16), (5.17). The resulting 
expressions are linear in e, and e,. We denote the coefficient of «3 by@d+4@,@,+ 4,65, 


that of «? by b+, e+ bye, that of « by €+¢,e,+¢,e, and the constant term is 
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denoted by d+d,e,+d,e,. Then equation (5.13) may be written 


(4 + Gy &y + Hy ey) 0. + (b+ dy ey + by ey) 02 + (E +0, e, + Fy ey) w + 


a ee 4 (5.20) 
+4-+ d,¢,+d,¢=0, 


) where, using (5.14), (5.15), (5.16) and (5.17), 


= — 51241 (a + 4a,)?, 
= — 2564} (a, + 4a,) (a, + 645), 
C = — 2564, a, (a, + 4a,), 
d = — 16a,(a,+ 4a,), 
oS 5124} (a, + 4a) (— 12a} — 52a, ay + 8a, ag + 104, as +- 24a, a5), 
by = 256 {4a a, (— 10a} — 52a, ag + 104, a, + 104, a; + 24a, 45) + 
+ at (— 8at — 44a, ay + 8a, a, + 8a, a; + 20a, a;)}, 
C1 = — 644, (19 at + 684, ag — 32a, 4; — 12a, a3 — 164, as), 
dy = — 32 (a, dy — ay 45), (5.21) 
@, = 512a5 (a, + 4a,) (— 3a} — 24a% a, — 18.0? a, — 48.0? a, — 56a, a2 + 
+ 48 a5 a, — 8a, dg ds + 20, ay 5) (a, + 2a), 
by = — 5124, (2a + 2603 a, + 8.03 ay + 3043 ay + 1764, a2 — 9648 a, + 
+ 11607 a} + 38a? a, a,-+144a% a, a, — 8a, a2 a, — 
— 8ai aya; — 604, a3 a5) (a, + 2a,)7, 
Co = — 64(1743 ag + 244? ay + 2484, a3 — 6443 a; + 13449 af + 1249 ay ag + 
+ 1924? a, a, — 72a, a5 a, — 32a, a3 a5) (a, + 2a,)+, 
d= — 8(128a3 + 8a? a, + 12a? a, + 72a, a2 — 64.25 as — Sara; — 
— 4, Ay ag + 90a ay M4) (4, + 24y)™ 


We write 
AH = Wy +  &y + Ho Eg, (5.22) 


where a, %,% are independent of e,, ¢, and are to be determined. Inserting 
(5.22) in (5.20) and expanding, we find, on neglecting terms of order higher 


than the first in ¢, é, that 


aa + ba® +C%+da- 
er 3402 +2bu9+6) +H 054 by 02 + F %q +4} + (5.23) 
OEE one eee ca eer 


Hence, since ¢,, é, are independent, we have 
aos + bot +ea,+d=0; (5.24) 
a4 (3 oR + 2b a» +6) +08 +b, 05 + 21% +4, =0, wes 


iy (3 B02 + 2B oy +2) + Ay 08 + beak + Fp Ot) +d, =0 
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Equation (5.24) is the equation found in classical elasticity theory for the 
frequency of Rayleigh surface waves. This may be seen if we denote 4 (a, +24) 
by ies 2a, by“, where jis jz are the Lamé constants. Then, from (5.21), we have 


a= 1678 (+a) 
b=— i A+i) 34+ 4f), (5.26) 
¢ = 8f(A +m) (A+ 2M), 
d= — (A+ji) (A+ 2p). 


Inserting these values in (5.24) and dividing throughout by (A+j), we find 
16j0° (A +j2) 08 — 877732 + 4 it) 0G + 87 (A + 2jt) % — (A+ 2f)=0. (5.27) 
Taking / %)=%), this equation may be written 


8ae— 8(2 —o)% + 8(1 —o)& — (1—0) =0, (5.28) 


where oa, Poisson’s ratio, is given by 
o =A/2(A+f). 
(5.28) is the equation found in the classical theory [2]. 


6. The Equation for m 


In this section we derive the first approximation to the equation (3.6) for m. 
From (3.6) and (4.18) we obtain 


Ag Com + {(Can + AG b) 0g? — [¢11 Cag + At 4B B® — (¢yp + Ai b)?) p*} m? +- 


: (6.1) 
+ (09? = P?) (CG? = 2, Of?) = O 
Writing 
Pleg=a«, wm/oq?—p?, (6.2) 
(6.1) may be written 
Az b Ce» B4 a {Cos “ile a3 b — [41 Cog + Ai Ay OA (Ges At b)? )?] 0} 62 + (6.3) 


+ (4 — & a) (14 — A? ba) =0. 
We write 
X = Ky + Oy C, + Hy Eo, 
B=Bo+fi%+ Bees, 
Az b Cog = Ry + Rie + Ry eo, 
Cag + dpb = Sot S,e+ Sze, (6.4) 
Ci Cont Ay Ag b® — (on + APD)? = t+ te + Te ee, 
Ay B Cyy = Hy + Hy, + Xp Co, 


2 
C1 + Ap b =@y +O, e+ O25, 


where a%,%,,% are determined by equations (5.24) and (5.25); Bo. Bi, By are 


independent of ¢,,¢, and are to be determined s+ R77S 6 eiGiuele See M1; 


t 
1 


1 R,, So,...,@_ are given by 
i, = — 164,45, 

| So= — 24,+ 845, 

| T = — 324, ag, 

eo = — 104 ae, 
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Mm, = — 24a,+ 8a,, 


4 
8 (a; — 204, a2 + 44, a, + 8a, 45), (6.5) 
8 (aj — 16a, ag + 4a, a3 + 4a, as), 
2(— 6a, + 16a, — 8a, + 245), 
Ry = 32 (— 2a} ay — 3.0} ag — 6.4, a2 + ay dy dg + 207 a,) (a, +- 245) >, 

2 

8 

8 

DP) 


aS 
I 


ay + 3205 + 8a, dy — 2a, dg + 240, Ay — 16g ag + 49 a5) (a, + 245)7, 


3 2 2 2 2 Spa = 
— a, — 64, a, — 4a; ag — 24a a, — 24a, as + 16.42 a;) (a, + 2a,)7, 


as 


— ai + 2at a, — 4a} a, — 1204 a,+ 8a2 a; — 4a, ay a3) (a, + 2a,)>, 


3a + 6a dy + 240, ay + 4 dy ds) (a + 2ay) 


Pn ae 


In deriving (6.5) we have used (5.9), (5.10) and (5.19). Inserting the values (6.4) 
in (6.3) and expanding, with the neglect of terms of order higher than the first 
In ¢,, €, we find 


Ro Bo + (So — To %) Bo + 1 — 9 &y + % HH 4 

E €{(4.Ry Bo + 2 So Bo — 2% % Bo) Br + Ri Bo + (Si — T1 T > %) Bo + 

Hb Hy 1G + 2x o Oty %y — Wy Hy — Wy %q} + (6.6) 
+ €2{(4 Ro Bo + 2 So Bo — 2% % Bo) Bz + Re Bo + (Sz — Te % — Tp %2) Bot 


+ ty 08 + 229 Xp XH — Wy Ay — Wo %} = 0. 


Since ¢,, é, are independent, this equation gives 


Ro Bo (So — 1%) Bo + 1 — Ooo + 40% =O; (6.7) 
(4 Ro BS + 2 So Bo — 2% % Bo) Bx 4 R, Bo + (Sy — TX T) 0%) Bo + (6.8) 
Tt ty 08 + ity Hy Oy — Wy % — WH, %—=0, 
(4 Ry B3 + 2S Bo — 27% Bo) B2 + Re Bo + (Sp — To %q — To Xs) BS + (6.9) 


H Hy 08 + 2x Xp XH — Wy Xe — We Hy =O. 


It is assumed that o,,,%, are known from equations (5.24) and (5.25). 
Having found f, from (6.7) we may then solve (6.8) and (6.9) for B, and fy. 
Equation (6.7), which is the equation found in the classical theory, may 
readily be solved. Inserting the values of Ky; Sp.%--, @) from (6.5) in (6.7) we 
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find 16 ay ay Bt — (— 2a + 8aq + 324, ay &y) By — (6.10) 
— {1+ (2a, — 8a.) % — 16a, a, 06} = 0. 

Denoting the roots by Boz, Box we see that they are given by 
B3=a—F, Bog =H — (A+ 28), (6.14) 


where we have written 44a,+ 84,= 1, —2a,=2. 


7. Results for an Incompressible Body 

In this section the frequency equation and the equation for m, found in 
Sections 5 and 6, for a compressible body on the basis of the second-order theory, 
are used in deducing the corresponding second-order equations for an incom- 
pressible body. The resulting equations are solved in explicit form and ex- 
pressions for the displacement components are then obtained. 

Rivitn [4] has shown that to deduce the second-order theory results for an 
incompressible body from those for a compressible one, the quantity (e,+ é,+ és) 
must tend to zero and ay, a, must tend to infinity in such a way that the product 
ay (e+ +s) and the difference a;— 2a, remain finite. The equation (5.19), viz 


C3 Cle {2a_/(a, + 2a,)} 2, 


yields ¢,+¢,+¢,=0, if we let a, tend to infinity, a, being finite. Thus we see 
that e, and e, may still be regarded as independent. Hence, the frequency 
equation for an incompressible body may be deduced from (5.24) and (5.25) by 
letting a,, a, tend to infinity but keeping a,—2a, finite. In the same way the 
equation for m may be found from (6.7), (6.8) and (6.9). 


Proceeding in this manner, (5.24) and (5.25) become 
128 a3 a8 + 96a} 0G + 164a,%,+1=0; (7.1) 
oy (— 96a} a§ — 48.47 % — 4a,) + 96a? (— 3a, + 245) «8 + 
+ 324a,(— 4a, +345) «8 — (11a, — 845) % =O, 
a (— 96a} ag — 484i % — 4a) + 96a? (— 3a, + 245) a8 + 
+ 32a, (— 5a, + 345) af — (13.4, — 8a;5) % = 0. 


(7.2) 


Equation (7.1) is the frequency equation found in classical elasticity theory 
for an incompressible body. It has a real root given by Ay %y = — 0.54787 ..: 
The remaining two roots are complex and do not give rise to Rayleigh waves [2]. 

The constants a, and a; are related to the constants C, and C, which occur 
in the expression for the strain-energy function for a Mooney material, viz 


W= C(I, — 3) +C(I,— 3), (7.3) 


by 
= — (Ci Gs), a; = — (Ci + 2C,). (7.4) 


Thus the real root a» of (7.1) may be written 


Op 0.54787 (CG at (7.5) 


. 
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Using (7.4) and (7.5), we obtain from (7.2) 


% = — (1.2007C, + 0.1049C,) (Cy + Cy), 


(7.6) 
a = (0.1049 C,=+ 4.2007C,) (C, 42C,)2 
For future reference, we obtain from (7.5) and (7.6), 
Oy he == (2015. Gye 0.1915 C,) (C, = C,) 4, on 


) Op ges (05 015.C, PAO Aas Gaede 
From (6.7), (6.8) and (6.9), we deduce 

2a, Bo — (1+ 44, o%) 63 + o%9 (1+ 24, %&) = 0, (7.8) 
; (2 — 8a, By + 84 %) Bo Py + 4(— a + as) Bo + {403.4 — 245) a + 
+ 4a, 04} BS +4(— 2a, + a5) «2 — 4a, a 90,—%,=0, 
(2 — 8a, Bo + 8a H%) By Bo + 4(— 2a, + a5) Bo + {4 (3. — 245) % + (7.10) 
+ 4a, 9} Bot 4(— a + a5) 09 — 4.4, tg % — 1 = 0. 


§ 


, Equation (7.8) has roots +f 1, +fo2 given by 
Bor = %: Bo2=% +2 a. (7.11) 
| Thus, using (7.4) and (7.5), we obtain 


Boo = 0.2188 (C, + C,)—3 


(7.9) 


(7.12) 


We denote the values of 8,, 8, corresponding to B§,=%, by Bi, Bs respec- 
tively, and those corresponding to Bjy=o%+24;', by Bye, By. respectively. In- 
serting 6? =f%,=a , in (7.9) and (7.10), we find 


2Bo1 Pur — 4 = 9, 2Bo1 B21 — % = 0. (7.13) 
Likewise, inserting 6? =$?,—=a)+a;", in (7.9) and (7.10), we obtain 


2 = 
2.45 Bos Bie — 2.ai X — aya + a,— a, =0, 


: (7.14) 
2.4% Bos Boz + 24; % — 4% + 2a, —a;=0. 
Using (7.1) and (7.2), it follows from (7.14) that 
Bra =3 Bost oe (7.45) 
Boe = % Bor %% - 
By (7.6), (7-7) and (7.12), we obtain from (7.13) 
By, = — (0.81104C, + 0.07086C,) (Cy + Ca)§, fone 
Boy = (0.07086 C, + 0.81104 Cy) (C; + Cy)~? 
and from (7.15), 
Bre pene (0.2397 Cy + 0.0209 C2) (Cy + C2) (7.17) 


Boo = (0.0209 C, + 0.2397C,) (Cy + Ca)? 
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Finally, we obtain expressions for the displacement components. For con- 
venience in the calculations we replace (3.7) by 
u, =tp(m, Ae-™* + m, Be~™”*) expu(p 4+ 42), 
ty = (Cyq + AZ) {0 g? — C41 P? + my Ay b) Ae + (7.18) 
+ (0.9 — %1 P+ my Az b) Be ™™} exp u(p 4 + 97). 
These values of u, and w, satisfy the equations of motion and give rise to the 
same frequency equation as do the values given by (3.7). The limiting value 
of C4 (49+ 420) + is 1, and that of A}b(c,,.+A{b)+ is 0. Hence the limiting form 
of (7.18) is 
u, =tp(m, Ae~™* + m, Be-™*) exp u(p y+ 92), 


(7.19) 
Ue = — 0G? (Oy + % & + Oy y) (A eW™™ + Be-™) exp t(p 4+ 91). 


It is to be understood here that the period 2z/¢ is assumed to be given and that 
it is independent of e, and e,, that #, m, and m, take their limiting values and 
that A and B depend on ¢, and e,. Since the body is incompressible, the density 
o does not change. The dependence of A and B on ¢, and éy is given by (3.16), 


i {(@ 9g? — ¢1 p?) + mi (Ap — Aj) 6 — cye]} A (7 20) 
ex ee {(o g? — ¢11 B®) + m3 [(A§ — Ah) 6 — Ge]} B. 
x it 
SD oA A tial eis Aves |. BEB es Bae ee (7.24) 


where Ay, A,,... B, are independent of e, and e,. Using (7.21), (7.20) may be 
written 


(Ag + Ay & + Ap eo) {(4 — % So — Mo Bor) — &1(H0S1 +05 So + 2% Bor Bir + 
+ 3 Bor) — €2 (%ob2 + &20o + 219 Bor Bas + 12 Bo1)} 
= — (By + By & + Bs es) {(1 — % Fo — Mo Boa) — &1(%oS1 + 0% 9 + 
+ 219 Bo2 Bre + Bos) — €(% Fo + H%2F9 + 2% Boo Boe + 12 fos)}, 


where we have written 
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Cy = Co + the ee, (43 — Aj) 6 — G2=—o— me —N2e,, (7-23) 
Co, 01, +--, being given by 
Ca SOs C1 =4(— a,+ 8a, — 4a,), 
2 = 44, (a, + 4a3 + 1244) (a, + 2a,)+, 
No = 4 (a + 242), ™ = 8 (a + 2a, — as — as), 
No =A(ai + 8ab4 Ody a, + 4a ag + 120, a, — 4a a3 — 4a, 45) (a, + 2a,)7. 


(7.24) 


In obtaining (7.24), we have used (5.9), (5.10) and (5.19). We observe that 
if we let a,, a, tend to infinity while keeping a3 — 2a, finite, then az'f, and 


=i a aA * 
az Np tend to 8, and a3" ¢,, a3 Cy, az, and az"y tend to zero. Hence the limiting 
form of (7.22) is 


(Ap + Ay & + Ap ep) {og + BR, + & (ay + 2Bo1Bi1) + € (% + 2Bo1 Bo)} 


=~ (By + Bye, + Ba) {+ Baa +64 (0s + 2BooPra) + &0(0a+ 2BoaBoa)}. > 
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» Thus, since e, and e, are independent, we have 
(% + Bos) Ao = — (% + Bis) Bo: (7.26) 
(% + B51) Ay ai 2Bo1 B11) Ag = — (% er oa) By, — (0% + 2Bo2Ai2) Bo, 
(%) + Bo:) Ag+ (%+ 2Bo1Be1) Ao = — (% + Bas) By — (%2 + 2Bo2 Bee) Bo- 
Using (7.11), (7.26) may be written 
Ay=—F 0 (% + 62.) By =— (14447! 05") By. (7.28) 
By (7.13) and (7.15), we have 
& + 2Bo1 Bis —20y, 1 + 2Bo2 Bie = a9 (% + Boo), 


(7.27) 


ai 2 (7.29) 
Oe + 2Bo1 Bor = 205, Oe + 2Bo2 Bos = Uy X% (% + Bos) - 
Hence, using (7.28), 
(% + 2Bo1 P11) Ay = 2m, Ap>=— NE (% Pas) By = — (% + 2Bo2Br2) Bo, (7.30) 
(%_ + 2Bo1 Ber) Ay = 24, Ay = — ty 9" (sy i Boo) By =— (%2 + 2Bo2 Boo) By. 
| Thus, from (7.27) and (7.28), we obtain 
Av] — kh B;; A,=—kB,, A,=—k By, (7.31) 
where & is given by 
k=1+ taj! 051 = 0.54369, (7.32) 
by (7.4) and (7.5). 
By (6.2) and (6.4), p2/oeq?=a=a +0%,¢,+%,¢. Hence if we write 
b=bho + bir + Po ee; (7.33) 


we find 
Po=(me9)*, pr=em%(Qg*/%)*, Pe = Fat (9 g*/t%)?. (7-34) 
Similarly, we write 
My = Moy + My1 Oy + Moy Co, My = Moo + Myo &y + Moo Co, (7.35) 
and since m?/0q?=?, we obtain 
my1 = (0 9*)? Bor, M2 = (0 9°)* Bos, M1 = (0 9?)? Bir, (7.36) 
m2 = (0 9°)* Bra, Mz = (0 9*)' Bas, M2 = (0 9?) Boo- 


By (7.44), (7.13), (7.15), (7-33), (7-34), (7-35), and (7.36), we obtain 


m, P = Moy Po + &1 (111 Po + Mor Pi) + &2 (M1 Po + Mo1 Pea) 
= 0 9 Boi 0h + 0 9? & (Bir %§ + $B or 1 Ho *) + 09? €2 (Boi %} + $801 %2 % 4) 
= 0G? (% + Hy & + My lo) = 097%, 

Me P= Moo by + (Me Py + Moz 1) + C2 (M22 bo + Moe bo) 
= 0g Bona + 0 97 % (Br2.%} + $B 2% % #) + 09? eo (Bo2%§ + FBo2 2% *) 
= 0 9? Boo te *(% +% & + Hs €s) = 0G" Box Bor. 


(7-37) 
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Thus, using (7.31) and (7.37), (7.19) may be written 
mh=—togka Ble ee Boo Bork te ™*) exp u(p x, +44), (7.38) 
Uy =o g?ka B(e~™* — kt e—™*) exp (fp % + qt). 
Now, by (7.11), (7-13), (7.34) and (7.36), when % and x, are not large we find 


Bo* B exp [e (by 1 + bo es) X1] EXPL — (M1 1 + Mer €9) X2] 
=1+ é,{B, Bai + tBy1 (0 q’)* x, — Bix (0 q’)? Xa} ag (7.39) 


+ €){By Bo' + ¢Bo1(0 97)? % — Berle gq?) Xo}, 
and 
Bo* B exp [«(p, & + Pe &2) %1] Exp [— (Myo €y + Moo &) X29] 


= 1+ ¢,{B, By + +61 (0 9*)* 4 — Bi2(o q?)* x2} + (7.40) 
+ {Bs Bo" + ¢ Box (0 9°)? x1 — Boo (0 9°)? %9}- 
Hence, from (7.38), using (7.39) and (7.40), we obtain 
tu, E exp[— t(a 0 g?)4 x, — gt] 
= [14 e{B, Bo* + 0% 06° + ¢By1 (0 9°)* 4 — Bir (@ 9°)* Xo} + 
+ €) {By Bo'+%2%0 ++ Bo1(09?)* %1—Bor (0 9?)* x2} ] exp{—(o 9*)} Bor x2} —(7.41) 
— Boz Bork [1+ 4 {B, Bo* + a4 %° + ¢ B11 (0 9")* % — By2(Q 9?)* xo} + 
+ € {By Bo’ +H, % + Bor (0 9?)8 % — Boo (0 9)* x2}] exp{— (0 9?)* Boo x2}, 
and 
uy E exp[— t(% 0 9)? x, — gt] 
=([1+ 4 {B, Bo" + a «9 + ¢By1 (0 92)! % — Bri (0 9?)* x} + 
+ €2{ Bz Bot ty «6 + Be5(0 9°)! 4 —Bo1 (0 9?) ) %}] exp [—(09?)* Boi x2] — (7.42) 
—R*[1 + {By Bo’ + a 46" + B11 (0 9?) )* x4 — Bra (0 9”) ee 
+ €){By Bo’ + % a9°+ tB21(Q 9°)* % —Bo2(o 9°)! %a}] exp [— (0 9?)* Boo Xe], 
where 
E=(p97 x, Bs) (7.43) 


4 % and aap’ are given by (7.7), By: and By, by (7.16), and B;, and By. by 
(7.17). From (7.32) we have 


k+ = 1.8393 (7.44) 
and hence, using (7.12), 


Boo Bor ®t = 0.5438. (7.45) 


In the case when e,=e,=0, (7.41) and (7.42) give, on taking the real part 
of the solutions, 


= |p (e~(e0")* Bor x a 0.5438 ¢~(e2)* Bos m2) Sin (Po mate qt), 
Uy a ka2 (e- (ea) *Bo: a eet 1.8393 ele")? Boe <) Cos (Po xy + q t) P 
as in the classical theory [8]. 


(7.46) 
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8. Love Waves 
A noteworthy feature of Rayleigh surface waves in classical elasticity theory 


is the fact that they do not exhibit any horizontal displacements transverse to 


the direction of propagation. Such horizontal displacements are, however, ob- 


| served on seismographs. To account for these waves which have the properties 


(i) of having transverse horizontal movement and (ii) of being transmitted super- 


ficially, Love [5] proposed to consider the body of the Earth as being covered 
_ by a thin layer of material having mechanical properties different from the material 


underlying it. Using this model, Love showed that such waves can exist and 
found that the essential condition for their existence is that the velocity of 
distortional waves be much less in the layer than in the underlying material. 

Here we shall first examine the possibility of having Love-type waves, pro- 
pagated in a principal direction, in a semi-infinite body which is subjected to 
a pure homogeneous deformation, without assuming the existence of a super- 
ficial layer. It will be seen that this is not possible unless a strong restriction 


is made on the strain-energy function. Then we shall consider the case when 


there is a superficial layer and find a condition essentially equivalent to that 
found by Love for the existence of these waves. 
Let the semi-infinite body be subjected to a pure homogeneous deformation 


_ (2.1) with which the stresses (2.2) are associated. Superimposed on this defor- 
- mation is an infinitesimal deformation (2.3). The equations of motion are given 


by (2.4). The stresses, given by (2.5), may be conveniently written 
3 
Gao =a 2 Cup UB Bp» 
p=1 
Gap = by, (da Mpa t+ Aaa) (B+ y +a), (8.1) 
b, = 2153 (W+ 2, Wr) . 

It may be pointed out here that b,> 0 since Ab, (0), Azo, (0) {4 +B $y #2} 
are each the square of a velocity with which a plane wave may be propagated 
along the x,-direction in an infinite solid of the material when it is subjected 
to the homogeneous deformation (2.1). 

We assume 

Ups te =O, Ug =] (%2) CXpal(p. yg); (8.2) 
where /(x,) is to be determined such that the motion is confined to the neighbor- 
hood of the surface, ¥,=0. The equation of motion not satisfied identically by 
(8.2) is 

22 b, D8 f + (9g — Arb b*) f= 0, (8.3) 
where D=d/dx,. The admissible solution of (8.3) is 
to OTE 
where A is a constant and m? is given by 
22 b, m2 = Az b, p? — 0 @?. (8.4) 


Then 
ug = Ae—™* exp t(p%+ 2). (8.5) 
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In writing (8.5) it is assumed that m?>0, or equivalently by (8.4), 

Li bs f2> OF: (8.6) 
Otherwise A—0, and hence /=0. From the remark made above, we see that 
(8.6) means that the square of the velocity of propagation of the surface wave 
is less than the square of one of the velocities with which a plane wave may be 
propagated along the x, direction in an infinite solid of the same material sub- 
jected to the same deformation. . 

The conditions to be satisfied on the surface of the solid are that it be free 
of traction both in the case of the deformation (2.1) and in the case of the de- 
formation (2.3). Putting /;=(0, —14, 0) in (2.9) and using (8.2) we find /;= (0, 0, 0). 
Hence from (2.7) and (2.10) the conditions to be satisfied on x,=O are 


Oz2=0, O3=0. (8.7) 
The first of these gives 
W, + (ly — 2) e+ BBW, =0. (8.8) 
The second gives 
Azb,m A =0. (8.9) 


If A=O then u4,=u,=u,=0. Also, A,=0 is impossible. Hence we have two 
cases to consider: (i) 6,=0, (ii) m=0. By (8.4) both cases give 


Jib, p= og. (8.10) 
In case (i) m may assume any value. When 0,=0 we have, from (8.1) and (8.7), 
W, = — AW, = AW. (8.11) 

In case (11) 
Us = Aexpt(px%,+q?2). (8.12) 


Equation (8.12) is an inadmissible solution since the motion is not confined to 
the neighborhood of the surface. 

Thus we see that only under the very restrictive conditions (8.11) are waves 
of the required type possible in the material. 

We now consider the case when a layer of material, with physical properties 
different from those of the semi-infinite body, is subjected to a pure homogeneous 
deformation, with principal extension ratios J;,, by forces applied laterally. The 
layer is then set on the semi-infinite solid which has also been subjected to a 
pure homogeneous deformation, with principal extension ratios 4,, by forces 
applied laterally. We assume that the principal directions for the two pure 
homogeneous deformations are the same. We neglect gravitational forces. Let 
the undeformed thickness of the layer be given by H/A,. The interface is given 
by x,.=0 and the free surface by x,=—H. The composite body is now sub- 
jected to an infinitesimal deformation. 

Symbols with primes refer to the layer and symbols without primes to the 
underlying body, and we assume that the displacement components are given 
by (8.2). Then the equations of motion not satisfied identically are 


Job, Df + (0g? Ab. p2)f=0 (x0), 


, t / U 
2” b; D® f + (0' gq? — 24°65 p%) f=0 (—H=S%,S0), oa 
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and the admissible solutions are 


ae (G Cee 08 0) 
f (x2) | (% 2 0) cw 
f(%_) = Acosm'x,+ Bsinm’ x, (-H<x,S0) 
where 
Aj b, m? = 22 b. p? — 0 q?, 
Pil 1%2P eg (8.15) 


A,” b, m’® = — Aj? b, p? + 0! g?, 
A, Bb, C being constants. 


It is essential that A{b,/o> 9q?/p2, in order that the waves do not penetrate 
far into the subjacent material. Without loss of generality we may take m’> 0. 
Note that we have assumed that 4)?3/0’<4q?/p?. Later we shall consider the 
case when A,7b,/0’ > q?/p?. 

At the interface x,=0, the tangential tractions and the displacement com- 
ponents must be continuous. The surface x,=-—H must be free of traction. 
Hence we have the conditions 


- of / 
093 = 3, Uz = Us, One 5103 


(8.16) 


-_/ 


Go, = 0, on My td 


We also have o2,=0 and o,,=0. 
The conditions (8.16) yield 


tan m’ H=— B/A, A2b,mC=—A*b,m'B and A=C. (8.47) 
On eliminating A, B, C, we obtain the equation 
tan m’ H = 72 b, m|(A," b, m’). (8.18) 
From (8.15), we obtain 


Ai by Q A? + 2 Q As? by 19 
= = = 1 yy! 2, 8.4 
tS eae ov ABD, oD) 


m 


In order that m may be real, the right hand side of (8.19) must be positive. 
Assuming this to be the case, we obtain from (8.18), 


13 by, pias ( os p> _ g Agr dy |b (8.20) 


FE le ? > 79 7 A 
Lah eemeaT SDC ALE GORD) oO, 


Then the value of m’ corresponding to any value of # must be such that tanm’H> 0 
A necessary condition that the right hand side of (8.20) be real is 


J? b,/0 > A,? b,/0’. (8.21) 


This condition is analogous to the condition found by Love [5]. Equation (8.20) 
may now be written 


(8.22) 


pie m! HI Q As b2tan?m A, OAs O4 | \" 


AB by (0' AF b.—@ Ay? bs) * (Q! AG ds —@ Ay be 
Thus we see that if (8.21) holds there is always a value of m’ corresponding to 
any given value of f. From (8.22) we see that as m’H increases from 0 to 7/2, 
pH increases from 0 to oo. 
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We see from (8.19) that we have dispersion, as in the classical theory. We 
also have the existence of m nodal planes in the layer if (2n+1)2/2>m'’H>nz. 


So far we have assumed /,7b;/0’<q?/p. If this is not so, then 
f(%_)=Ce"™* (%, 20), 


(8.23) 
}(%_) = A cosh m’ x, + Bsinh m' x, CVE AA 


The boundary conditions (8.16) give 
tanh m’H = B/A, Rb,mC=—A,?b,m'B and A=C. (8.24) 


Hence 
A,” b, m’ tanh m’ H + 23b,m=0. (8.25) 


It is necessary that m> 0, but since m’ tanh m’H is always positive there is no 
relevant solution in this case. 
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Introduction 


The nuclear reactor is based on the fact that by joining two subcritical bodies 
T’ and T” one may obtain a critical or supercritical body T. The neutron 
transport properties of such a 7 as a whole are entirely different from those of 
T’ and IT”, though the properties of 7’ and T” can be considered unchanged. 
The different behavior of such a T is due to the coupling between 7” and T” 
caused by the neutron transfer through the contact surface between 7” and T”’. 
The reverse procedure is possible as well, z.e., every critical or supercritical body 
can be split into subcritical parts. 

It is natural to ask how the properties of a given body T depend on the 
properties of its parts J’ and 7’? This question is the fundamental problem 
of the reactor theory. It is usually answered only indirectly by determining the 
life histories of neutrons in the body T by the neutron transport theory, cf. [J, 2]. 
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As a consequence, such an answer does not give anyteaepers Per of Fi 
expressed explicitly in terms of the analogous properties ol DT sandd: . Such an 
explicit answer, however, would be very desirable from the theoretical point 
of view. This paper presents an attempt to obtain such an explicit connection) 

The neutron transport properties of the body observable from the outside 
are determined by the relation between the distributions of incoming and out- 
going neutron currents. These generalized reflection properties of the body give 
enough information to justify asking how the reflection properties of the body 
T depend on those of its subcritical parts 7’ and 7”. 

The two-group treatment of this problem has been given for the case of 
simple geometries [3], and a similar treatment was recently found useful also 
in the case of electron diffraction, cf. [4]. We shall treat it for arbitrary geometries 
in the frame of the linear transport theory. It is sometimes difficult to under- 
stand the theory if the physical model and the mathematical formalism are not 
separated clearly. Therefore, in Chs. I, I, III, and IV, we shall give the mathe- 
matical description of the physical situation and some of its consequences. Later 
on, in the Appendix, we shall collect a few mathematical derivations. 

In Chapter I the reflection properties of a body JT are described in linear 
approximation by an albedo operator A * and a distribution of independent out- 
going neutron current g. Two cases are treated: (i) the unsteady case in a finite 
time interval (§ Ia), and (ii) the steady case (§ Ib). 

In Chapter II we study the reflection properties A and g of the body T 
composed of two parts, J’ and T’’, which are covered with an absorbing layer 
of negligible thickness. Assuming quite general reflection properties A’, g’ and 
A", gq" of T’ and T” respectively, in § Ila we derive an inhomogeneous equation 
of the form (1—/*A’'A’’)x=y, where x and y are closely related to the out- 
going and incoming neutron currents of T, respectively, and where 4 is the trans- 
parency factor of the boundary layers. This equation is the mathematical 
analogue of the physical situation. By studying its properties we shall learn 
about the reflection properties of the compound body J. The solvability of 
the corresponding homogeneous equation is the condition that T be critical. 
By analysis of the inhomogeneous equation we shall construct addition formulae 
for the reflection properties of 7, T’, and T’’. The results of Chapter II are a 
generalization of those in [3]. 

Studying the reflection properties of T as a function of the transparency 4, 
we shall obtain further information about how the reflection properties of T 
depend on those of 7’ and T’’. In Chapter III, assuming the scattering of 
neutrons to be uniform, we are able to show that T is subcritical, critical or 
supercritical according as A<Ay), A=Ay or A>Ay, where A) is determined ex- 
plicitly by the albedos A’ and A’. The critical condition A=Ay is not directly 
expressed by the solvability of the critical equation. It might have some practical 
use for reactor design in the case of albedos A’ and A” which are determined 
experimentally. For the case AJ), A<d,, we shall derive some simple and 
plain results which are in accordance with the known facts of reactor physics 
regarding the outgoing neutron currents of T. 


* For some applications of the albedo concept see [7], [8]. 
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Nearly all necessary mathematical derivations are collected in the Appendix. 
There are some consequences of the non-negativeness of the neutron current 
distributions which have an important role in these derivations, cf. [2]. There- 
fore, we shall interpret these distributions as elements of a cone K in the Banach 
space L, chosen so as to match as well as possible the logical structure of our 
derivations. Since the inhomogeneous equation is found to be a mathematical 
analogue of the physical problem, the analysis of it is essential for the mathe- 
matical part of this paper. This analysis is carried out in two steps: i) assuming 
the corresponding operator to be additive and to map K into K, and ii) in addition 
assuming it to be some kind of a uniformly positive operator, cf. [2]. Such an 
analysis of the inhomogeneous equation is done here for the first time, although 
the analysis of the corresponding homogeneous equation has been given earlier, 


Eferiic, 0). 
Chapter I. Reflection properties of a single body 


a) Unsteady case in a finite time interval [0, ty] 

For simplicity it is assumed that the body T has a time-independent shape 
and a finite, piecewise smooth, boundary B. This boundary is crossed by neutrons 
of different velocities whose magnitudes and directions we denote by v€ [0, c] 
and by unit vectors s, respectively. At each point r€ B the solid angle Q=4z7 
is divided into two parts, 2. and 2,. The part Q_ is formed by all directions s_ 
of the incoming neutrons, while the 
part 2, is formed by all directions s* 
of the outgoing neutrons (Figure 1). 

Weadopt the following notation: 

i) dB for the neighborhood of 
the point r on the boundary B. 

ii) dQ for the element of solid 
angle around the direction s point- 
ing from Yr. 

iii) dJ for the average number of the neutrons passing through dB during 
the time interval [¢, + dt]. The magnitudes of the velocities of these neutrons 
lie between v and v--dv, and their directions belong to dQ. 

iv) 1=i(v, r, s,t) for the neutron current density, where v¢€[0,c], rcB, 
s€Q, andt€[0,t]. The relation of 1 to dj is given by 


dJ =i(v,r, s,t)dvdB dQ dt. (1.1) 


The distribution z(v, r, s, t) depends on the external conditions and on the 
properties of the body. However, it is always non-negative and such that the 


total number of neutrons 
c 


aff s sa 


to 
fitv,r, s,t)dvdBdQ dt (1.2) 


0 


going through B in the time interval [0, f] is finite. Therefore, the unique 
properties common to all physically possible distributions are assumed to be 
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40 and J<oo. Such distributions form a closed set K in the Banach space 
I=1,(v€[0,c], rEB, s€Q, t€[0,%]) (1.3) 


with the norm ||| =J (V. 4, 2, 3)*. 
We define the incoming neutron current of T by 


(1.4) 


and the corresponding subspace J_ cI by 
Li pera, Bases 
All physically possible distributions 7_ form the set 
Kia ele to 


The set of all non-negative functions belonging to the Banach space I {-) will 
be denoted henceforth by me: 


In an analogous way as for 1_ we define the outgoing neutron current by 


1. eS —— Se 
(1.5) 


= (0), —- Se 


One part of the outgoing neutrons 7, is due to the incoming neutrons 7_, so 
that 7, depends on 7_. In order to have the distribution of the incoming neutrons 
1_ independent of the outgoing neutrons 7,, we suppose that re-entrance of the 
outgoing neutrons 7, is prevented. Hence, the non-convex body T is under- 
stood as submerged in an idealized absolutely absorbing medium. If in the 
body T there are some empty areas which are considered as its outside, they, 
too, are understood as being filled up by an absolutely absorbing medium. 

Let us denote by g€ K, the outgoing neutrons 7, corresponding to the incoming 
current 7_=0. The current q is due to the independent sources inside the body 
T and to the neutrons which came into the body T during the time interval 
t<0. Accordingly, the incoming neutrons 1_ produce the outgoing neutrons 7, — q. 

For a finite time interval [0, fy] any distribution 7_¢€ K_ of incoming neutrons 
is possible. This distribution 7_ is transformed by the body T into a uniquely 
determined distribution 7,€K,. This property of the body T determines, by 
the relations 


Ai=Ai,+Ai_, At,=0 and Ai-=i,—g, (1.6) 


the operator A whose domain is K. The operator A is a generalization of the 
concepts of the albedo and of the transparency, which are known from the 
transport theory, cf. [1]; it will be called the albedo. A and g completely describe 
the reflection properties of the body. 


* In the text we shall denote by (V....) those sections of the Appendix that 
contain the appropriate mathematical tools. 


| 
| 
| 
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As in the linear transport theory we shall neglect the mutual influence of 


_ the neutrons and their influence on the properties of the body. Hence we shall 
| assume: 


i) The distribution A7_ of outgoing neutrons corresponding to any distribution 
7_ is always non-negative. We express this fact by the relation 


AIS (x) 
ii) If 7’, 7’"CK, then 
A(t’ +4") = Ai’ + Ai”, (1.8) 


An operator having the properties (1.7) and (4.8) will be called additive (V. 4). 

In the time interval [0, ¢)], as a consequence of the incoming current 7_, |7_| 
neutrons enter the body, and ||47_| neutrons go out. The ratio ||A7_||/||7_| depends 
on the distribution 7_ but is finite for a real body. Its least upper bound 
is the norm 


|Al = sup [4 ¢ [el = sup [4 +/]- (1.9) 


of the albedo A. Therefore the albedo of a real body is a bounded operator. 
This follows from its additivity alone (V.5). Only in the case when ||A||>1 
can there occur distributions 7_ such as to cause more outgoing than incoming 
neutrons in the time interval [0, 9]. 

Some of the incoming neutrons will proceed straight through the body, while 
the rest will be scattered, making it possible for them to be absorbed or multi- 
pled*. The outgoing neutrons produced directly and indirectly by incoming 
neutrons will be described by additive operators D and S, respectively: 


Ae DEE Ss (1.10) 
In the case of an empty body g=0, S=0 and A=D. 


Let the neutrons 7_ form a narrow beam incident at a very small angle on the 
convex side of the body just after ¢=0. Nearly all of them will have crossed 
straight through the body by the time ¢=¢,, and we have ||Dz_|<||z_|. Hence 
it follows from the meaning of the norm that for the real body 


[ID =1, 
so that from (1.10) we get the estimate 
134]. 
The albedo A and the operators S and D are continuous operators; 7.¢., for 


all 2,,72¢ K we have 
|44— 4%] 5 [4] la — al (1.11) 
and analogously for S and D. We can prove this fact immediately if we extend 
the operator A to the whole space J by the relation A(t, —1,) =At,— At, for 
all 1,,7,€K. 
b) Steady case 
Let the body T have time-independent properties, and let the incoming currents 
i_ be stationary. Some of them produce uniquely determined steady outgoing 


* In this paper it will be convenient to interpret the multiplication of neutrons 
as a scattering event. 
Arch, Rational Mech. Anal., Vol. 8 27 
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currents 7,. Such currents 7 will be called permissible. The situation on the 
boundary B is now described by a distribution 7€ K, 


K={i: 1€L,(v€[0, c], re B, s€Q); +20}, 


as in the unsteady case. 

In order to have i_ independent of 7, we disregard any possible re-entrance 
of outgoing neutrons 7,. In the unsteady case we distinguished only one type 
of real bodies. In the stationary case*, however, there are three types of real 
bodies T. We shall characterize them in accordance with the reactor theory: 


i) The subcritical body T: Each distribution 7_€K of the incoming neutron 

current is permissible. We describe the reflection properties of subcritical bodies 

in the same way as in the non-stationary case. 

However, now the current q is due to the in- 

dependent sources only, and in the case of 

elastic scattering ||Ai_||=||7_| for each 1_€K. 

If the subcritical body T is interpreted 

as compounded of two subcritical parts T”’ 

and 7”, their reflection properties A’, q’, A” 

Fig.2 and q’’ completely determine the reflection 
properties A and g of the T. 

ii) The critical body T: Each critical body can be imagined as compounded 
of two subcritical parts T’ and T’’. For instance, we can obtain subcritical 7” 
and T”’ by slicing up the body T and regarding every second slice as a part 
of T’ (Figure 2). Evidently, if the slices are thin enough, they are subcritical 
by themselves. Consequently, the parts J’ and 7” are also subcritical, since 
neutrons 7’ and 7,’ are not allowed to re-enter T’ and T”’, respectively. 

The properties of the critical body T are not completely determined by the 
properties A’, A’’, q’ and q” of its parts T’ and T”’, because some of the outgoing 
neutrons 7, need not be produced by q’, g’ or i_. The shape of this current, 
but not its intensity, is determined by the properties of the constituent parts. 

ii) The supercritical body T: Some of the distributions 1_€ K are not per- 
missible. Just as in the critical case, we might imagine each supercritical body 
as compounded of two subcritical parts. But the reflection properties of the 
supercritical body are completely determined by the reflection properties of these 
parts. 

Since we interpret the critical and supercritical bodies as compounded of 
two subcritical parts, the properties of such compound bodies grow interesting. 
They will be considered in Chapters II and III. 


Chapter II. Reflection properties of compound bodies 
a) Steady case — General considerations 
We shall consider the body T as compounded of two finite subcritical bodies 
I” and I”. Let their boundaries B’ and B’ be such that their contact surface 
Bo B" can intersect a straight line at most d times. We consider as the boundary 


* It is to be noted that from the physical point of view it would be more adequate 
to consider the steady case as asymptotic to the unsteady case as f, ty +00. 
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surface B of the compound body T the surfaces B’ and B”, excluding the contact 
surface B’n B"”, 


Br Bw Be IB Ouse 


The physical situation on the surfaces B’ and B” is given by the neutron currents 
mek oc!’ and 176k? cl". respectively. In what follows the spaces J’ and I” 
will be interpreted as the subspaces of the Banach space 


I, = {i,} =L,(v€ [0,c], ré B’UB", sEQ), (2.1) 


4.€., 
L— Die 1,=0, rE B” — Bare 


and similarly for I’’. Let P’ and P” be the projectors associated with the sub- 
spaces I’ and I”; 


4! a Pa 4" = (ip and Pp! Pp! = Pp" jal. (2.2) 


We assume that the properties A, and q, of T’ and A, and qg, of T’”’ are known 
and that we do not change them by compounding 7’ and JT” into T. Before 
we put 7’ and T” together we cover them with a layer of negligible thickness *. 
Due to the absorption of this layer the intensity of the neutron current crossing 
the boundary B’ or B”’ is diminished by a factor A€ [0,1]. Therefore the pro- 
perties of TJ’ and T” are changed, and their new properties A’, A”, qg’ and q’’ 
are given by 

MeN Arh RA ALAS ng =a Agr and wag = 1 gah 
4.€., 
t_=q +tA’t_ and t1,=q"+A"%s". (2.3) 
If the boundary layer contains independent neutron sources, additional terms 
are to be added to the q’ and qg”. 

I’ and J’, being interpreted as subspaces of J,, the ranges of A’ and A” 
are contained in K,. It is convenient to extend also their domains K’ and K” 
to the whole set K,. For this reason we multiply the right-hand sides of A’ 


and A” by P’ and P”, respectively. These new albedos will be denoted again 
by A’ and A”. It follows that 


A'=P'A'=A’ je and A”X= PY'A"N=A" P".” (2.4) 


After these changes, which have only formal character, the equations (2.3) 
remain unchanged, and the albedos A’ and A” become additive operators 
having the domain K,. 

On the contact surface B’. B” the outgoing current of one constituent part 
is equal to the incoming current of the other one, due to the transfer of neutrons 
between 7’ and T”; 

peterand ai—4, jor re BlaB”. (2.5) 


As in Chapter I let us assume that on the surface B the outgoing neutrons 7, 
and 7, can not re-enter JT. Consequently, on the boundary B arbitrary distri- 
butions of incoming currents 7” and 7” are possible. They define the incoming 


* To such a layer we attribute absorbing properties only. 
Dap 
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current 12 of I, 4 = i re [Be Be 


i eee (2.6) 


-? 


=n We Dye. 
Similarly the outgoing current 7, of T is defined by 
i,=t,, reB'—B", 
=i’ Pe BB (270 
Ne ASN eos eo 


The set K of all physically permissible distributions 7_ and 7, belongs to the 


subspace 
deo t= 0; re BoB". (2.8) 


Let the projector P correspond to the subspace J. Then the equations (2.6) 
and (2.7) read 
i_= P@_+7) 
and (2.9) 
+, =P, +14,); 
The real body contains no absolutely absorbing material. Hence, for later 
reference, we may assume that for some m 


P(A'A"\"i, £0 or PA(A'A" +O if Pi, +0. (2.10) 


The incoming current 1_ is independent of the outgoing current 7, , due to the 
assumption that the outgoing neutrons 7, are prevented from re-entering T. 
The reflection properties of JT are determined by the dependence of 7, on 7_. 
We shall express them in terms of the properties A’, A’, g’ and q’’ of T’ and 
T”. In order to do this we shall first determine, using (2.3), the dependence 
of 7’, and 2’ on i_. For this purpose we must first calculate 7”. and 7”. 

On B’— B” the current 7”. equals the incoming current 7_; on B’7 B”’ it equals 
the outgoing current 7’. In terms of the projector P’ this gives 


i_=P'i_+#). (2.14) 
Similarly we get 

ho PHe tea Ni (2.12) 
From (2.3), (2.11), (2.12) and (2.4) it follows that 


‘= + Ai +4) 

and ES, aan s: (2.13) 
ee Pla OR ES 

Thence, by use of (2.4), we obtain 


(eA Alp Oe sey (2.14) 
with 
Mts Qu gear) ae ae 
The current 7’, satisfies an analogous equation. It is sufficient to solve one 
of these equations only, since the solution of the other one follows then from (2.13). 
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Equation (2.14) expresses the current 7’, by its distribution P’’7, on B’O B". 
This distribution is determined by the equation 


(1 ss P" A’ A”) Wea Mie PE (2.15) 


which is derived by multiplying (2.14) by P’. Equation (2.15) is of fundamental 
importance for the treatment of the properties of T. Namely, by solving it and 
oe (2.14) and (2.13) we get the outgoing currents 7, and 7’ explicitly. Since 
y€K',, we are justified in interpreting (2.15) as an inhomogeneous equation 
ee with the additive operator P’’ A’ A” P’ which has the domain P” K’, = 

Hook 6 (Ve 42): 
In order to characterize the compound body T, we must answer two questions: 


i) What currents i_ are permissible when the body T is a compound of 7” 
and T”’, z.e., for what currents 7_ does (2.15) have a solution? 


ii) Are the properties of T completely determined by the given properties 
of its parts; 7.e., is the solution of (2.15) unique? 


Regarding these questions we can say: 
1) If (2.15) has a solution at all, the series 


Pl => P(A’ Ay (2.16) 
j=0 


will be a solution also (V. 12). Consequently, only those currents 7_ are permissible 
for which the series (2.16) converges. If the series (2.16) does not converge even 
for 1_=0, then the compound body T is physically impossible. 
ii) The solution of the (2.15) is unique only if the corresponding homogeneous 
equation, 
(RAGA A) Be aaO. eens ee er Py ies, (2.47) 


has no solution (V. 14). If (2.17) has a solution, the body T is critical. Namely, 
due to the assumption (2.10), there can be a part of the outgoing neutrons 17, 
which is not produced by q’, g” or z_. Its distribution is not determined by the 
properties of JT’ and 7” completely, but up to a positive constant factor. There- 
fore (2.17) is the critical equation of the compound body 7. Thus, only in the 
critical case do the properties of the compound body fail to be determined by 
the properties A’, A’’, q’, and q’’ of its parts (see § Ib). 

Let us derive the series (2.16) from the physical picture for the special steady 
case with q’’=0, P’g’=q' andt_=0: The neutrons q’ go through B’  B” into T”’. 
One part of them returns (cf. the footnote on p. 384) through B’ B” into T’, 
and thence one part returns to 7”’ and so on. From the meaning of the albedos 
A’ and A” we infer that those neutrons which originate from qg’ and leave T” 
through B’. B” after being 2m —1 times reflected between TJ’ and T”’ are given 
by P’’(A’A’’\"q’. Consequently, 

PUG PEG GAGA CG SUPr(AY AY) go, 
which is the same as the series (2.16) adopted for the case P’” y=q’. Similarly, 


one can construct the series (2.16) for the general stationary case as for the 
non-stationary case. 
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The series (2.16) gives a clear idea of how the part T" influences the increase 
of the outgoing current P’’7’,. The term P’’(A’ A”)"y of the series (2.16). represents 
the neutrons produced by P’’y (cf. the footnote on p. 384) and entering for the 
(+1) time the part 7’. The norm ||P’’(A’A’)" P’| is the least upper bound 
of the ratio |P”(A’A”’)" y||/|P’’ y]. It has the property that either |P"(A"A")"P [21 
for n=1,2,..., or lim |P’(A’A”)”P’|=0 (V. 6). Of a special interest for us 


will be the latter case, which holds if, and only if, A< A, (V. 6), 
ie dees 
Jy = lim |P"(A, A)" P'[ 4” = lim |P"(A, Ay)" PY **. (2.18) 


The constant A), which is the critical value of transparency, satisfies the in- 
equality 
dy > |P" Ay P' AP’ Ag PY. (2.19) 


Its connection with the critical compound body will be considered in §§ IIb 
and IIIb. 

If the series (2.16) converges for 7 =0, the compound body T is physically 
possible and is subcritical, critical or supercritical. Let us consider these cases. 


b) The subcritical case 


If T is subcritical, (2.15) has a unique solution for each 1_€K, and vice 
versa. Using the series (2.16), we get from (2.15), (2.14), (2.13) and (2.9) the 
connection 

i, =q+Ai_ (2.20) 


between the incoming current 7_ and the outgoing current 7,. This connection 
is expressed in terms of the independent outgoing current of T 


Gc (F pal (ge) (2.21) 
and of the operator 


A=PA'+P(1+A4) (1—A” A)AAL+ AY 


= PA"+ P(A+A”) (4 — AAV Al A), (2.22) 


with 


(= ArAY) 


ime 


(A’ A’)! 
i] 
Operator A in (2.20) is the albedo of the compound body T and is to be inter- 
preted as an additive operator with the domain K. Formulae (2.21) and (2.22) 
are the addition formulae for the reflection properties. They connect the albedo 
and the independent outgoing current of T with the same properties of the 
constituent parts 7’ and 7”. 
It is evident that T is subcritical if the transparency A, determining the 
coupling between 7’ and T”, is small enough. Moreover, we prove that the 
compound body T is surely subcritical if 


Ade (2.23) 
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for in this case (2.15) always has a unique solution (V. 6, 15). Equations (2.19) 
and (2.23) show that the compound body is subcritical if 


A<|P’ A, P"|-*|P" A, P'H. (2.24) 


We shall pause here to give some examples of subcritical compound bodies: 

i) Clearly, the compound body T is subcritical if there is no mutual contact 
between 7’ and T”, 7.¢., if B’~B’ =0. In this case P’ P’ =0, and P=P’+P", 
and equations (2.19), (2.4), (2.22), (4.6), and (2.21) give 


Ay = 00, AeA EA” “and, g=q' +9" 


ii) The neutrons are not multiplied in T’ and T’’, and B’ is concave and 
B" convex at their contact. Therefore, ||P’ A” P’|=1, P’’D’ P’ =0 and|P” S’ P” | 
<1. The compound body T is always subcritical, since the condition (2.24) 
now reads A<||P’’S’ P’’|-#>1 and is automatically satisfied since 2 € [0, 4]. 

ili) Let it be supposed that there is no scattering in T’ and T”, 1.e., A, =D, 
and A,=D,. We shall show that 


Te AD ele = aU eee (2.25) 


where is the smallest integer greater than $(d—1). The relation (2.25) follows 
from the meaning of the operator P’’(D,D,)" P’ which relates (for A=1) the 
neutrons P’2” entering T” through B’O B” with the neutrons P’’(D,D,)" P’ i” 
leaving T’ through B’°B”. Before the neutrons P's” produce the neutrons 
P’'(D,D,)" P’i” they will have passed ” times through 7’ and T’’. Consequently, 
the surface B’~ B” is intersected 2u-+-1 times by their straight trajectories. Since 
it was assumed (p. 386) that this can happen at most d times, we get the relation 
(2.25). From (2.18) and (2.25) it follows that ~Aj=oo. Thus the compound body 
considered is subcritical for all values of 4, and its albedo 


PDE Pie) (eee) UD ce) (2.26) 


is a finite sum of different products of D’ and D”. 
The albedo A of the compound body T can be separated into two terms 
D and S, as has been done in Chapter I, equation (1.10). The neutrons D1_ 
passing straight through the body T are not affected by the scattering properties 
S’ and S” of the constituent parts. Therefore the operator D depends on D’ 
and D” in the same way as the albedo A of the third example. The scattering 
operator S is determined by 
S42), (2227) 


where A and D are given by (2.22) and (2.26), respectively. 

Properties of a subcritical body with an arbitrary number of constituent 
parts can be also determined from the properties of its parts by (2.20), (2.21) 
and (2.22). To do so, we join two arbitrary parts into a new part and proceed 
in this manner until the whole body is built up. 

If the outgoing neutrons 7, do not contribute to the incoming ones 7_, the 
relation between 7_ and 7, is given by (2.20), with A and g determined by (2.21) 
and (2.22). Sometimes we are interested in the relation between 7, and 1_ of a 
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non-convex body TJ, allowing the outgoing neutrons z, to contribute to the 
incoming ones i_. In this case we consider a compound body whose parts are T 
and some suitably chosen empty body. We introduce the empty body to create 
the necessary contribution of the current 7, to the current 7. 
If it is desired, in the first example above, to take into account the transfer 

of neutrons between the separated parts T’ and T” or a possible re-entrance 
of outgoing neutrons, we proceed in a similar way. 


c) The critical case 


If T is critical, (2.17) has one or more linearly independent solutions (P" 1a 
and vice versa. Incoming currents 7_ which satisfy P(A’ A" by za(P'4,); for 
any 7, and some «>0 are not permissible, because the corresponding series 
(2.16) do not converge. If the incoming current 7_ is permissible, from (2.15), 
(2.14), (2.13) and (2.9) we can determine the corresponding outgoing current 1, 
by using the series (2.16). This relation between 7, and z_ can always be written 
in the form of equation (2.20), where A is determined by (2.22) and 


g=(P+A4) (¢ +9") + Daj PUtA)A(P4);. (2.28) 
y] 


The constants «; are positive and characteristic of T but are not determined 
by the properties of JT’ and T’’. The operator A of (2.20) can now be inter- 
preted as the albedo of the critical T; its domain is the set of all those 1¢K 
for which the series (2.16) converges. 


d) The supercritical case 
If the compound body T is supercritical, for some 7_€K (2.15) may have 
a unique solution, and vice versa. If the incoming current 7_ is permissible, the 
outgoing one 2, is given by (2.20), (2.21) and (2.22). The operator A of (2.20) 
can now be interpreted as the albedo of the supercritical T, and its domain is 
those 1_€ K for which the series (2.16) converges. As in the subcritical and the 
critical case this albedo A can be divided into D and S. 


e) Unsteady case in a finite time interval 
The analysis of the unsteady case is analogous to that of the steady one. 
In particular, we get relations which are analogous to those of § IIa. Then it 
follows from the physical picture that the equation analogous to (2.15) must 
always have a unique solution. We shall give a proof of this fact in § IIIc. 


Accordingly, the further analysis of the unsteady case is analogous to that of 
the subcritical case of § IIb. 


Chapter III. The properties of a compound body when the scattering of 
neutrons is uniform in some specified sense 
a) Introduction 
In the foregoing chapters the properties of a compound body have been 
treated under very general assumptions concerning the properties of its parts, 


as given by equations (1.7), (4-8), (1.10) and (2.25). However, these assumptions 
are not sufficient to deduce statements about the real compound body T such as: 
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i) In the unsteady case with a finite time interval each distribution i_¢ K 
is permissible, and the properties of T are completely determined by the prop- 
erties of its parts. 

ii) In the critical case no physically non-trivial distribution of the incoming 
current 7_ is permissible for steady conditions. 

To deduce these and similar statements about a compound body T and to 
get a better understanding of the influence of the transparency 4 on the prop- 
erties of T, some new properties will be attributed to the albedos A’ and A”. 

The scattering of neutrons on atomic nuclei is nearly isotropic, and the 
scattering function of a nucleus is continuous. If the body is large enough, the 
outgoing scattered neutrons have all possible velocities, and their distribution 
is nearly independent of the velocity distribution of the incoming neutrons. Besides, 
the neutrons which arrive together in a body do not leave it at the same time. 
As a consequence, the scattered neutrons are, in general, distributed more uni- 
formly over the boundary surface and over the time interval than the incoming 
neutrons. The supplementary mathematical properties which will be attributed 
to A’ and A” are suggested by just these facts. 


b) The steady case 

Let us assume henceforth that T” is such that the distribution of the scattered 
outgoing neutrons S’i’ =7' (v, r, s) produced by a narrow beam of monoenergetic 
incoming neutrons, 7_=6(v — vp) -0(r —Yy) - 6 (S$ — Sp) is a continuous function of 
v, 1, S, Uy, 7% and S,. Let the same be assumed for 7’. As a consequence: 

i) S’ and S” are integral operators having continuous kernels. 

ii) S’ and S” are compact operators. Because of (1.10), (2.3), and (2.25) 
the operator P’’(A’A’’\™P’ is compact, too. 

ii) The operator S’ satisfies the relation 


Sa AAs) SEK. (3.1) 


The operator S” satisfies an analogous relation. From this and from (1.10), 
(2.3), and (2.25) we get the relation 


PAA y<|P' yA" sy, nem and yeK,, (3.2) 


where s,€ P’’ K’, and depends on the exponent only. 


Since distributions of neutron currents are non-negative, there is always 
a relation opposite to (3.2), namely 


Syke Sa eed uA yy ot Ss Ghee, (3.3) 


where «, is a positive constant chosen separately for each distribution P” y= 0. 
If the compound body T is subcritical, it follows from (2.27), (2.22), (2.26), 
(1.10), (2.25), and (3.1) that its scattering operator S is compact and satisfies 
a relation analogous to (3.1). 
The operator P’’(A’A’’)” P’ contains P’’(A’A"’)"P’ [n/m9] times as a factor. 
Due to (1.10) and (2.25), every factor effects at least one scattering. Consequently, 
the neutrons P’’(A’ A’’)”y have been scattered at least [/m)] times. After each 
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scattering the distribution of the neutrons is smoother and in some sense more 
uniform. Thus we expect that the greater the number [7/9], the more alike will 
be the distributions P’’(A’A”)"y corresponding to different distributions P’’y 
with norms ||P’ y|=1. In what follows we suppose (c/. [2]) that with increasing 7, 
the above distributions P’’(A’A’’)"y become so much alike that it is possible 
to choose an integer m,=m, and a distribution s€ P’ K’, such that the relations 
(3.2) and (3.3) are satisfied for all yEK., if we put s=s,=s, and n=n, (see 
V.18 et seg.). Multiplying the relations (3.2) and (3.3) by P’(A’A”)”, we see 
that analogous relations are valid for all n> m9. 

In the case when 4,<1 the reflection properties of T depend on the trans- 
parency / as follows: 

i) for Ay<ASX1 the compound body is supercritical, 

ii) for A=A, the compound body is critical, 

iii) for OS A<A, the compound body is subcritical. 

i) If A> Ap, only those currents q’, g’’ and 1_ are permissible for which P’” y=0, 
since only if P’”’y=0 does (2.15) have a solution. This solution is P’’7',=0, 
and it is unique (V. 19, 6, 23). With P’” y=0 and P’7,=0 we get from (2.14) 
and (2.13) the relations 
and * ; : (3.4) 

Cae ee ae 

The condition P’’y=0 is satisfied only if the currents qg’, g’’ and z_ do not 
cause any transfer of neutrons through B’ Bb” from T”’ into T”’. As a con- 
sequence, a supercritical real compound body T without empty areas is physi- 
cally possible only if 7. =0, g’=0 and qg’’=0. Indeed, if z_ =£0, since the scattering 
of neutrons is nearly isotropic, some of the incoming neutrons will be scattered 
through Bb’. B” from T’ into T”’. Besides, if a part of such a compound body 
has an independent outgoing neutron current, its distribution is always such 
that there is a current of neutrons through any part of the boundary surface. 
Therefore the domain of the albedo A of such a supercritical T is an empty 
set and g=0 (cf. § IId). Thus, a supercritical compound body with q’ +0, g’’ +0 
or 7_=-0 is not possible unless it contains suitable empty areas. 

ii) If A=A,, only those gq’, qg’’ and i_ are permissible for which P” y=@ 
(V. 19, 6). Further, the critical equation (2.17) has a unique solution (P’’7,), 
(V. 23, 25). With P’’y=0 we get from (2.14) and (2.13) the relations 


t=y+aA' A"(P"G,)y 
= q" +A" +a(P"i1)), 


where « is a positive constant, characteristic of the compound T and not de- 
termined by the properties of 7’ and T”’. 

In the case of a real compound body T without empty areas, the equation 
(3.5) can be simplified by taking into account the physical meaning of the 
condition P’”"y=0. We get for the domain of the albedo A the empty set (cf. 
§ IIc) and for the independent outgoing current g the expression 


q=aP(1+ A’) A"(P'G,)o. (3.6) 


and 


(3.5) 
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in) If A<A,, the compound body T is subcritical, as we have shown in 
§ ILb, where also the addition formulae for the reflection properties were given. 
Nevertheless, the body T under consideration is subcritical if and only if A< A). 
The outgoing current P’’;’ 


, satisfies for all y¢K’, the relations 


Ne—1 


Hp ed mr , 11\9 A ), ana 173% 
Pi > YS P'AA Vy tty ete (P ae 
i—0 ¥, 
and noes A (3.7) 
Prat < 5) P(A Ary! yp GATT Mes (Pi )o, 
j= 140 


where A€ (0, do), and é is the ¢, of the relation (3.3) with y=(P’7,), (V. 24). 

Let us now consider the changes of the reflection properties of the subcritical 
T when 4 tends towards 4). It follows immediately from (3.7) that the outgoing 
current P’’s’, produced by any P’’y-0 becomes arbitrarily great, provided 4 


is close enough to A). Moreover, if AA, its shape P’’7’,/||P’’2’| approaches 
the solution of the critical equation (2.17) (V. 25, 22); 7.e., 

: P''7, ||P’ a] > (P" 7) 

5 [Pred P's, ia 


AA, and P’y=+0. 


The same holds for 7’, and 7”. 

If Ay=1, the first case is impossible. If Ay<1 the second case is impossible 
too, since A€[0,1]. Therefore, the critical condition for the compound body 
considered is 

fp Aeand. A= A: (3.9) 


After the foregoing considerations the operation of a reactor can be interpreted 
as follows: 

Due to spontaneous fission of nuclei the real subcritical body always has 
independent sources q’ and q” of very low intensity. If A is close enough to the 
critical value A,, these sources can still produce an arbitrarily large outgoing 
current 7, even if 7_=0. Its shape 7,/||7,|| can be approximated by the solution 
of the critical equation; 7.e., if ASA», then 


i,[|i| ~ P+ A) AM(P’ 2 )ol|P(A + 4’) AP" oll. (3.10) 


This relation is derived from (2.14), (2.13), (2.9), (3.8), and (2.10), taking account 
of (3.7). 
c) The unsteady case in a fimite time interval 
Let us assume the following properties of the albedos A’ and A”: 


i) The distribution 7, of the outgoing neutrons at a certain moment is not 
affected by the distribution i” of the incoming neutrons at a later time; 7.¢., the 
albedo A’ satisfies the causality condition 

Ali_==0 “forall, 26 [Ost), 
if ” 14 
; i =0 forall #€[0,4]. 5 


The same holds for the albedo A” and for S’, S”, D’ and D”. 
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ii) The scattering operators S’ and S” satisfy relations analogous to (3.1) 


It is evident that any real body satisfies the causality condition (3.41). The 
second property assumed is very general and can be attributed to every rea! 
body without loss of generality. From these two assumptions it follows thai 
the operator P’’A,A,P’ satisfies the causality condition and a relation whick 
is analogous to (3.2)*. As a consequence, Ay=oo (V. 41), and the condition 
(2.23) tells us that, for all values of 4, (i) the properties of T are completely 
determined by the properties of its parts, and (ii) any incoming current 1_€Kis 
permissible. 

Chapter IV. Resumé 


The steady reflection properties of a subcritical body have been described 
by means of a linear albedo operator A and a distribution q of independent 
outgoing neutrons. It has been shown that the boundedness of A follows from 
its linearity and from the non-negativeness of the distribution functions of the 
neutron currents. The same holds for the unsteady case in a finite time interval. 

We have considered some general properties of a compound body whose 
parts T’ and T” have been covered with a negligibly thin layer of absorbing 
material. It has been shown that in the steady case such a compound body is 
completely characterized by the properties A’, A’, g’ and q’’ of its parts, with 
the exception of the critical case. The albedo A of a compound body and its in- 
dependent outgoing neutron current g were expressed by A’, A’”’, g’ and q’”’ in the 
subcritical, critical and supercritical cases. However, in the critical case the 
independent outgoing current of a compound body is not completely determined 
by the properties of J’ and 7’’. The compound body was found to be subcritical 
if the transparency J of the boundary layer of its parts is lower than some 
critical value /,; this value 1s determined by the albedos of the uncovered 
parts T’ and T”. 

Assuming that the scattering of neutrons is uniform in some specified sense, 
we obtained a clearer picture of the relation between the properties of the 
compound body and the transparency / of the boundary layer. It was possible 
to deduce the following statements about the compound body: 


i) In the unsteady case in a finite time interval the compound body is 
completely characterized by the properties of its parts, and any distribution of 
incoming current is permissible. 


ii) In the steady case, there are the necessary and sufficient conditions for 
the compound body to be subcritical, critical or supercritical. They are expressed 
in terms of transparency of the boundary layer and of the albedos of the 
uncovered constituent parts. Thus we obtained the critical condition, which is 
not connected directly with the solvability of a critical equation. 


ii) If, in the steady case, the transparency A of the boundary layer tends 
towards the value A), the compound body is growing critical; i.e., the outgoing 
current becomes infinite with (Aj—A)+, and its distribution approaches the 
critical one. The critical distribution is determined by the critical equation up 
to a positive constant factor. 


x We note that, due to the causality condition, the operator P” A, A, P’ can satisfy 
a relation of the type (3.3) only in the trivial case when SO: 
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Chapter V. Appendix 
Introduction 
The mathematical formalism used in this paper for the description of the 
_ physical phenomena consists of 

i) a set K of all non-negative functions belonging to a space L,, 

ii) linear operators A which map K into K. They will be called additive 
operators. 

The set K will be treated as an abstract closed subset of a Banach space, 
characterized by four properties. If K is a set of all non-negative numbers, 
these properties are: 

i) a non-negative number multiplied by a positive number remains non- 
negative, 

li) and ili) the sum of two non-negative numbers is a non-negative number 

which is never smaller than any of the summands, 

iv) a bounded infinite series of non-negative numbers is convergent. 

The boundedness of an additive operator and some usual properties of its 
powers will be proved first. Then the estimate ||A B]<|A|| ||B|| will be generalized 
to prove that an operator which satisfies the causality condition and the relation 
Ax<\||x||a for some a€K, and for all x€K, is quasi-nilpotent. 

The properties of the inhomogeneous equation (1—AA)x=y, where A>0, 
and x, y€K, are of decisive importance for the physical problem treated. It 
will be shown that: 

i) The solution of the inhomogeneous equation is unique if, and only if, the 
corresponding homogeneous equation has no non-trivial solution. 

u) If there is a solution of the inhomogeneous equation, the Neumann series 
is also a solution. 

iii) The Neumann series converges uniformly if A is smaller than the spectral 
radius A, of A. 

Special attention is given to the spectral properties of a continuous compact 
operator A which, for all x¢€.K and some a€K, satisfies the relation ¢(x)a< 
Ax<||x|| a where e(x)>0 if x=-0. It will be shown that: 

i) The corresponding inhomogeneous equation has a unique solution x for 
all yEK if A<A,. But, there is no solution for any yCK, y+0, if AZ Ap. 

ii) The corresponding homogeneous equation has a non-trivial solution %, 
only for A=A,. This solution is determined up to a positive constant factor. 

iii) The solution x of the inhomogeneous equation, corresponding to any y +0, 
has the property that ||x|>oo and x/||x|—> x9/|xo|| for A>”y, A< Ay. 

1. Definition. Let K be a closed set in a Banach space B. Let the elements 
of this set K have the following properties: 

4. if x€K and if a>0, then ax€k, 

2. if x, yCK, then *+ y€K, 

3. if x, yCK, then ||x+ y|=|x|, 

4. if the subset {x;}9°C K and if lim 

co 


n 


De; 


1=0 


Cc 
<oo, then the series >’ x; con- 
i=0 


verges, 1.¢., >, %,CK. 
i=0 
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2. Lemma. 1. Jf x, —x€K, then x=0. 

2. If x, y, x —yEK, then |x|=Ily]. 

Proof. 1.1f x, —«x€K, it follows from (1.3) that 0=|¥—x|2|x 
x=0. 

2. If x, y,x«—y€K, then x=y+2, where z€K, and we get from (1.3) the 
relation |«||=||y+z|=|ly|- Q-E.D. 

3, Example. Let B be L,(0,1). Then an example of & is the set 


Ky = {r() exp @ p): 1) EL, (0, 1); |e) —«@| SB}, (1, 


where «(¢) and f(¢) are real functions such that 0O<«(#)<2a and OSA (t)Sa/4 
for t€[0, 4]. 

For any value ¢, of the argument ¢, the values of all functions of the set K, 
form the sector | argz—« (¢,)| <A (4) in a complex plane. Therefore, the elements 
of the set K, have the properties (1.1), (1.2) and (1.3). It is easy to prove that 
the condition (1.4) is satisfied too. This can be done by introducing an auxiliary 
coordinate system with the coordinate vectors e,=exp(« (é) — 2/4) and ¢)34 
exp (a (t)-++-2/4) and by making use of the inequality 


(x+y)? SxPty?,  x,yZ0, pei, 


and of the completeness of the space L,. 


> hence 


Analogous examples of K exist in the space /,, too. 


4. Definition. In this paper, the operator A with the domain K is callec 


additive if AKCK (1 
and if A(x+ty)=Ax+Ay, x, yEK. (2) 
5. Lemma. An additive operator A is homogeneous and bounded, i.e., 

AGB =AAS, f= Orn x ore. (1) 
and A(x/||x||)<const., forall x€K. (2 

Thus A has a norm defined by 
[Al =sup | 4 x[/ |x|. (3° 

x€K 


Proof. 1. It follows immediately from (4.2) and (1.1) that the relation (5.4) 
holds for all non-negative rational values A, of A. 


2. There is an infinite sequence 


Man Ay Cis: 


BD |e 3s be 


AGN oe teens (4 


if A is not bounded. Now, considering (4.2), (4.3), (5.1), and (5.4), and using 


the series © 


D %,(WEK, 
n=l 
we have the contradiction 


[4 > x,/n? 
n=1 


Hence any additive operator A is bounded. 


co > 


= lim ||A(x,,/2?)| = lim n = oo. 
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3. Since A is bounded, we have for any A>0 
A(A x) =m {A[(A—A,) x] +4, Ax}=AAx. Q.E.D. 


6. Alternative. If A is an additive operator, the following alternative holds: 
either 


Ath efor all Pn Ao, (1) 
or 
lim ||A”| =o. (2) 
Furthermore, the limit rig 
2 1 
Jy =lim |A*|-* 2 [AJ > 0 6) 


exists, and the operator 1A (A=0) satisfies (6.1) tf AZApo, or (6.2) if A<dy. We 
shall call the number i, the spectral radius of the operator A. 
Let A and B be additive operators. Then 


lim |(4 B)"|-* = lim |(B.A)"|-*. (4) 


Proof. 1. Let be ||A%|=qg<1. Then we choose an integer N such that 
q\||A —1]. Thence we get for the power A” (n=n,(N-+m) +7, 
w—0), 1, 2,....) the-estimate 


|A"| SA" |e" [Alf < 9” 


having the limit (6.2) as a consequence. 


2. The operator JA satisfies the relation (6.1) if, and only if, A= oy al m 
It is known [6] that 


sup [A"|-* lim |A"|-* 
3. The relation (6.4) follows from the estimate 
1 1 EN 
[4 By = (ANB Bay Oe 
and from the one we get by interchanging A and B. 


7. Definitions. The direct sum of N Banach spaces E;=({x;}, 1=1, 2,..., N, 
is denoted by 


By = {42 K=(%y, Xo, -.+, Hy) fs (1) 
The norm of x is given by 
1p 
Ikl=(SbvP) 224 Q) 
The projectors P, are defined by 
P, x = (0j1 %, O52 %2) +++» O;n %y)- (3) 
We have N 
ae (4) 
i=1 


A mapping of BY into a space 


1=(zr" ow 


eae x! = (4X1, XB, ++, Xy)s x;€ (— 00, 0) 
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is defined such that to each *€ Be corresponds an x CN having the components 
x; = |B, «|. (6) 


This mapping is isometric; 7.¢., 
lel = (7) 


Let KY be a closed set in By which has the properties (1.1), (1.2), (4.3), (1.4) and 
POE AES ae ee Mee (8) 


Further, let A be an additive operator with the domain K}¥. Then an associated 
matrix A’ is defined by 


Aj,,=|RAP|, 154,75N. (9) 


A’ is to be interpreted as a linear operator with the domain J. 
8. Theorem. Let A, (k=1, 2, ..., m) be additive operators with the domain KY, 
and let Aj, be the associated operators with the domain l¥. Then, for any integer m, 


we have the estimate 
™ 


IT Aj}. 


k=1 


™m 


[] A, 


k=1 


Proof. From (7.6), (7.4) and (7.9) we get the estimate 
N 


(1) 


= 


N 
| <Y Al xt =(A'x');. 
jaa 


(A x),= 


RAPP x) 
j=l 
By its use we get from (7.7) and (7.5) the estimate 
m m—1 ' 
IT A, < |4s( IT 4. x) 
k=1 k=1 


m—2 


sas 41-( Fas) 


k= 


Nas 


= 


laa 


™m 
TT Ane’ <| 
k=1 


ih 
This proves (8.1). - 


9. Lemma. Let T be the linear matrix operator with the domain LX and having 
the elements 


AeA OSes 
=0, j>t, 18%, jSN. . 
Then its n™ power, T”, is given by 
- Ca = il Ra eeie 
(Tj =| ers ) f= 4; 
i= 4 
=0 74 2) 
and satisfies the relation 
|Z" |S Nmax (T),;< (w+). (3) 


Proof. 1. When =2, the relation (9.2) can be checked directly. For an 
eaten! integer ” we prove (9.2) by complete induction. It follows from (9.1) 
that 

N i 
("),7=(T I"); == LL yea 


=1 k=1 
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From this expression and from (9.2) with the exponent —1, the relation (9.2) 
with the exponent » follows immediately, provided 7>7. For 7<7, the same 


follows from 
Ca) ce oem em 


k=j ge 2 n—1 


2. The estimate (9.3) follows from (9.2) and from Hoelder’s inequality. Q.E.D. 
10. Theorem. Let A be an additive operator with the domain KX, whose associated 
operator A’ has the properties 
and 


(1) 


max Aj;=q<1. 
v7 
Then 
lim | A" =0. 


Proof. From (8.1), (10.4) and (9.3) follows 
lim |A”|| < lim |(A’)*| Slim gq” |T”| Slim g” (n + N)* =0. 


11, Example. The Banach space L,(0, 1)={x(¢)} can be interpreted as a 
space apy which is a direct sum of an arbitrary number N of the spaces 


E,=1,(53 <)> petuT a aENG (1) 


In this case, the set K, given by (3.1) is an example of the set kK}. As an example 
for Theorem 10, we consider the additive operator S having the domain Kk, and 
satisfying: 

i) the causality condition; 7.e., if 


w= 0 ator all” elicit ClO, Tike Vek, (2) 
then 
Sa—0 forall ve lO. and 
il) the relation 
|els —S*eK, (3) 
for all x € K, and some Sek, 
Namely, if N is large enough, the associated operator S’ has the properties (10.1). 


Proof. 1. It follows from (7.9), (7.3), (14.1) and (11.2) that 
Si=|RSBI=0,  4<7. 
2. From (7.9), (5.3), (7.3), (7-2), (2.2), and (11.3) we get the estimate 
Si;= sup |F SF yl xl Slhs 


, 


and with |s|€L,(0, 1) it follows from (10.1) and (11.1) the relation 


N>o 


iN Up 
lim qs lim max ||P s|| = lim max ( h |s(@)|? i) = ANG OF Dens) 
ay 
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Let D be an additive and nilpotent operator, 
D* =05. tte, (5) 


which satisfies the causality condition (11.2) and has the domain K,. Then 
the operator A=S-+D is quasi-nilpotent; 7.é., 

lim ||4” A” || =0 
for each A€ [0, 00). 

Proof. Since S and D satisfy the causality condition, the operator (4A)” 
satisfies this condition, too. From (11.5) and (11.3) we infer that the operator 
(A.A)" satisfies a relation which is analogous to (11.3). As a consequence there 
exists, for each 2€[0, c), an integer N such that ((4.A)”)’ has the properties 
(10.1). Hence |(AA)"”||<1 if the integer m is great enough. Therefore, in 
accordance with Alternative 6, lim \(A.A)"|=0 for any value of A. Q.E.D. 


12. Lemma. Let 
(1—A)*=y, yeK, (4) 


be an inhomogeneous equation associated with an additive operator A having the 
domain K. Then the Neumann series 


Xap = af i (2) 


as a solution of (12.1) af, and only if, x,,€K, cf. [6]. 

Proof. 1. It is evident that this condition is necessary, since K is the domain 
of A. 

2. ane, 


LAty— AD Aty—y| 5 (1+[Al) lim] 3 4%y|=0. QED. 


13. Theorem. Let x,¢K be a solution of (12.1). Then x,,€K and is, there- 
fore, a solution of (12.1). The difference 


4% — Xqp = lim A" x%,€K 
n 
is a solution of the homogeneous equation 


(1—A)x=0, (1) 
corresponding to (12.1). 


Proof. 1. Writing the solution x, in the form 


HE VE Aly a ARIS AR yeh tS AL ean Bey (2) 
and taking account of (2.2), we get 
lim 2 4éy| < |x, |< 0. 


Thence and from (1.4) it follows that x, »€4, and as a consequence of Lemma 12, © 
the series x,, is a solution of (12.1). | 
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2. From (13.2) we get the relation 
m+n—-1 


2 4's 


which shows, due to the convergence of the series x, p» that the sequence {A” *,}9° 
is a Cauchy sequence. Since K is a closed set, it follows that lim Aine ds, 
n 


Using (13.2), we prove that 


WAae a, = A” x,|| —_ | 


*1 — %a» — lim A’ x, | <lim | A” x, —lim A‘ x, | + lim 


DA‘y| Sale 


It follows from x,—x,,€K and from 


y = (1 — A) x = (1 = A) Xap + (1 — A) (44 — Map) = ¥ + (1 — A) (4 — Xap) 
that the difference x,—x,, is a solution of (13.1). Q.E.D. 


14. Theorem. If the inhomogeneous equation (12.1) has a solution, then this 
solution 1s unique if, and only if, the corresponding homogeneous equation (13.1) 
has no non-trivial solution. 

Proof. 1. It is evident that the stated condition is necessary. 


2. Let the homogeneous equation (13.1) have no non-trivial solution, and 
let x, be a solution of (12.1). Then, as a consequence of Theorem 13, the series 
%g» 18 a solution of (12.1), and the difference x,—x,, is a solution of (13.1). 
Therefore %, — x,,=0, which proves that the stated condition is sufficient. Q.E.D. 


15. Theorem. The series x, 1s uniformly convergent if it satisfies the Cauchy 
condition of convergence uniformly for all elements of the set {y: yCK; |y| =1}. 
The condition lim ||A"||=0 is necessary and sufficient for the uniform convergence 


of the series x,,. If this condition is satisfied, the series x,, 1s the unique solution 
OindeA)a (cf. (6 |). 

Proof. 1. It follows from Alternative 6 and from the definition of the uniform 
convergence that the stated condition is necessary. 

omit lim |A”||=0, we have 


n+m d 
[Sa4iy 


|SbIZI41< bl Poni sees 
Pais 7=(|n/n,| 1= 


sb 5 Sla-par=blo(3 


1 
a 
4=0 


Af) /(4 —q), m>0, n>0, 


where 7, is chosen such that 
q=|A"| <1. 


This estimate tells us that the series x,, satisfies uniformly the Cauchy test of 
convergence. According to Lemma 12, the series x,, is a solution of (12.1). 


3. Let x be a solution of (13.1). Then 
4—A%\x=0, =0,1,2,3,+5- 
If lim || A”|| =0, it follows that ||x||<||x||lim | A" =0, and according to Theorem 14 


the series x,, is the unique solution of 2-4)2 


28* 
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16. Theorem. Let A be an additive operator with the domain K, and let (1—AA), 


A=0, be an operator with the domain eee fel AA) x € K}. Then when A< Ag 
there is an inverse additive operator (1—AA)+ with the domain K, which 1s given 
by the uniformly convergent serves 


(1-24) => aay (1) 


The same holds if (12.1) has a unique solution for any yCK; however, the serves 
(16.1) need not be uniformly convergent. 

The proof follows immediately from Definition 4, from Theorem 15, and from 
Alternative 6. 

17. Example. As an example for Theorems 13 and 15, let us consider the 
additive operator A given by 


A x(t) = x(@2), O22 gS 41. 


Let its domain be the set K, with «(¢)=0 and B(¢)=0 (viz, 3.1); 7.e., the set | 


of all non-negative functions belonging to L, (0, 1). 


t 


| 


Let x(t)€K, be the positive non-decreasing function determined by the | 


equation x(t) — x(at)= y(t). Since y(¢) CK,, it follows from Theorem 13 that 
x (t) = %q, (2) + lim x(a" t). 
From Theorem 15 and (6.3) we infer that the equation 
x(t) —Ax(at)=y(t), Aza, 


for all yCK,, has a unique solution belonging to K,. This solution is given by 
the uniformly convergent series %,,. 


18. Definition. An additive operator A having the domain K will be called | 


uniform (cf. [2, 5]) if there is an integer 7, and an element s€K such that 
i) for all xCK the relation 
|x||s — A*®xEK (1) 
is satisfied; and 
ii) for each xCK, x=£0, there is a positive constant ¢ such that 
A™x%—es€K. (2) 


19. Theorem. Equation (12.1), which is associated with a uniform operator A, 


has a solution for a y,CK, y,=-0, if, and only if, lim||A"|=0. This solution is 


n 


unique. 
Proof. 1. It follows from Theorem 15 that the stated condition is sufficient. 


2. In order to prove that the stated condition is necessary, let us assume 
that (12.1) has a solution for some y,€ K, y,=-0. Then it follows from Theorem 13 
and (12.2) that lim |A"y,]=0. From (18.1), (18.2) and (2.2) we get, for any 


y CK, the estimate 
A” y| = 4%" A" |< | yA" 5] S |] 4” yall/er- (1) 
Hence we have 
lim || A” Ss lim ||_A” Vy |e =—=(() OED, (2) 
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20. Lemma. Let the spectral radius A), of a uniform operator A be finite, and 
let yEK, y+0. Then A" y+0 for n=1, 2,.... 


Proof. Let A” y,=0 for n=n’' and y,€K, y,--0. Then it follows from (19.1) 
that A”—0 for all n=7’, i.e., 44 =00, which proves the lemma. 


21. Theorem. Let the spectral radius hy of a uniform operator A be finite, 
and let yEK, y+-0. Then 


lim  ||A"(4—4A)+y|=0, tgs eae 
A<hy, A> Ay 


and the approach to the limit is monotone. 
n+m 


2 (Ay A)*y 


for any yCK, y=+0. For, in the opposite case it would follow from (1.4), from 
Lemma 12, and from Theorem 19, that lim |(4)A)”|| =0, which would contradict 


the definition of A). Thus we prove the theorem asserted, using Theorem 16 and 
the estimate 


ts 


Proof. The norm 


and the integer m tend to infinity together 


[oe) 


24) y] = 
Namely, here the right side is as great as desired, provided m is great enough 
and A is close enough to A), and, due to (2.2), the left side grows monotonely 
with A. 

22. Lemma. Let the spectral radius i, of a uniform operator A be finite. Then 
the elements 


{x(n, A): x(n, A) =|A*(1—1A)+y|4A4*(1—14)4 yCK, 


n+m 


D (AlAo)' (A 


1=Nn 


n+m 


x (Ao A)'y 


=n 


A)iy|| = alae 


fateh ator 


1 
VEOe = 0, Ae 2 ony Ae oe (1) 

satisfy the relations 
jim (1 — 9 A) x(n, A) =0 and |x(n, A) =1. (2) 


Proof. The existence of elements x(n, A) follows from Theorem 16 and from 
Lemma 20. We prove the limit (22.2), using Theorem 21, as follows: 


dim |(1 — Ap A) #(m, A)|| SA” [|] Jim |" 1 — 2.4) y | + [Al im (2, — 4) = 0. 


23. Theorem. Let the spectral vai Ay of a uniform Reni A be finite. 
Then the homogeneous equation 


((— AA) x=0, KEK Six S15 AEH &, 00), (1) 


may have a solution only for A=Ap (cf. [5]). If there is a solution xq, it is unique 
and satisfies the relations 
Apes — %y EK (2) 
and 
—MregseK, e&>0. (3) 


Proof. 1. It follows from (23.1) and (2.1) that the eigenvalue A must be 
positive. If 4, is an eigenvalue with an eigenvector %,, then it follows from 
(23.1) that lim |(A,.A)”%||=1, and as a consequence of Alternative 6 we get 
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4,24). But if A, were greater than A), then from the limit 


Lim fay A)" x] = lim (Ay/4)" [2 4)" aa = 0, 


and from (19.2) we should get lim (A) A)”|| =0, which contradicts the definition 


of Ay. Consequently, the real eigenvalue can only be 4). Relations (23.2) and ; 
(23.3) follow immediately from (18.1), (18.2) and the fact that 2)> 0. 


2. It remains to prove that %,=%, if x, and %, are the eigenvectors of the 
eigenvalue J). To this purpose we need the supremum of all real constants a 
which satisfy the condition x,—ax,€K. All such constants satisfy the relation 
a<e51, which follows from (23.2), (23.3) and (2.1). Since a=O satisfies the 
stated condition, the supremum exists, and, K being closed, 

%, — (supa) x,E€K. 
Applying the relation (18.2) to the difference x,— (supa) x,, we get 
x, — (sup a) %,— esEK. 
From this relation and from (23.2) it follows that 
%— (supa+t ely”) x.EK. 
Comparing this relation with the definition of sup a, we infer that «=0, and as 
a consequence of Definition 18 we have x, — (supa) x,.=0; 1e., x,=%,. Q.E.D. 


24. Theorem. Let the spectral radius i, of a uniform operator A be finite, 
and let x» be tts eigenvector. Then, for all yEK, OSA</,, and m=O, 


[e.e) 


i 7) Xo Mam 
2 Aly s Wn EK (1) 
and 
AJ Ag) *™ |y||/e ~ ; 
(A/ a 0 X> > (A A) yeK, (2) 
w i=n,+m 


where &, equals the ¢ occurring in (18.2), with x=y. 


Proof. Taking account of (23.2) and (23.3), we derive from the relations 
(18.1) and (18.2) the analogous ones 


|| Ao” eo" %y — A* xEK (3) 
and 
A™% — elo x EKR. (4) 


Using them, we can immediately prove the theorem stated. 


25. Theorem. Let the spectral radius Ay of a uniform operator A be finite, 
let A be continuous (i.e., for all x, yEK we have |Ax—Ay|—0 if |«— y||—0), 
and let the n," power of A be a compact operator. Then dy is an eigenvalue (cf. [5]), 
and 


%> = jim x(n, EK, Ws OF42 TS (1) 


1s the corresponding eigenvector, the norm of which equals 1. 
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Proof. First we prove the assertion for N=MN3, considering the infinite 
sequences 


X; = {% (mg, A;,): R=0,1,2,...5 lim 2;,= Age 


Since the operator A is continuous and K is closed, it follows from (22.2) and from 
Theorem 23 that the limit x;,= lim %(n3, 4;,) of each convergent sequence X; 


is the eigenvector x) belonging to the eigenvalue A. The existence of one con- 
vergent sequence X;, 1.¢., the existence of x), follows from (22.1) and from the 
compactness of the operator A”. 


To complete the first part of the proof, we have to show that any sequence 
X; converges. If a sequence X; does not converge to xy, it contains a subsequence 
X; whose elements satisfy the relation ||x(x3, A;,) — x9|>6>0. But, due to the 
compactness of A”, the sequence X; has to contain a subsequence converging 
to %), which contradicts the definition of the sequence X;. Hence any sequence 
X, converges to %), proving the theorem for n=7,. 


Deriving the relation 


jim le(7, A) — £0, A) == 0, 5 st 0,152,000, 


from Theorems 16 and 21, and from (22.1) and (1.3), we prove the limit (25.1) 
for all integer x. 
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Lineare Operatoren mit positiver Inversen 


JOHANN SCHRODER 


Vorgelegt von L. COLLATZ 


Bei vielen Abschatzungen méchte man von Ordnungsbeziehungen zwischen 
bekannten ,,Ausgabewerten‘‘ auf entsprechende Beziehungen zwischen den zuge- 
horigen ,,Eingabewerten‘‘ schlieBen. In welchen Fallen ist dies méglich ? Abstrakt 
formuliert: Es sei M ein Operator mit halbgeordnetem Definitionsbereich und 
halbgeordnetem Wertebereich; wann folgt wv aus Mu<xMv? Einen Operator | 
M mit dieser Eigenschaft nennt L. Cottatz [4] ,,von monotoner Art“. 


Bei Gleichungen der Form 
Mu=r (0.1) 


lassen sich z.B. oft Schranken fiir eine Lésung w* angeben, falls diese Aufgabe 
(0.1) von monotoner Art, d.h. falls der Operator M von monotoner Art ist. 
(Siehe L. Cotratz [4] und [6], S. 42f.) Insbesondere erhalt man dann Fehler- 
abschatzungen fiir eine Naherungslosung u). Dabei mu8 vorausgesetzt werden, 
daB w* existiert. Mit Hilfe von Fixpunktsatzen kann man ferner in vielen Fallen 
die Existenz einer Lésung u* von (0.1) beweisen und Schranken dafiir herleiten, 
falls nur ein M benachbarter Operator M, von monotoner Art ist [/6]. 

Wir betrachten hier im wesentlichen nur lineare Operatoren M. Der obigen 
Frage entspricht dann: Wann folgt wo aus Mu=Zo? oder: Wann besitzt M@ 
eine ,,positive’’ Inverse? Die monotone Art nichtlinearer Operatoren laBt sich 
vielfach aus der monotonen Art linearer Operatoren folgern (s. Abschnitt 6.1). 
In einer anschlieBenden Arbeit sollen jedoch nichtlineare Operatoren direkt 
behandelt werden. 

Probleme einer monotonen Art, welche sich nicht genau in der genannten 
Weise abstrakt beschreiben la8t, treten bei parabolischen Differentialgleichungen 
ante ($.7B.112), 118,15), [13]). Auch bei Matrizen arbeitet man zum Teil 
mit anderen Relationen (z.B. [7]). Wir verallgemeinern daher die Fragestellung. 

In einem archimedischen halbgeordneten linearen Raum fiihrt die gegebene 
Ordnungsrelation uo zu weiteren Beziehungen u>o0, w>o und uo. Fiir 
endlichdimensionale Vektoren wu = (w') z.B. bedeutet bei der iiblichen =-Definition 


u>o: alle u' sind =0, aber nicht alle sind =0, 
u>o: alle uv’ sind >0. 


Mit diesen verschiedenen Ordnungsbeziehungen definieren wir acht Klassen der 
monotonen Art linearer Operatoren. Wir zeigen, daB nur zwei davon ,,wesentlich 
verschieden™ sind. Eine dieser zwei Klassen enthalt die oben genannte monotone 
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Art, die andere wird durch die starkere Forderung ,,4>o folgt aus Mu>o‘ 
gekennzeichnet. Die letzte Eigenschaft findet man z.B. bei Randwertproblemen 
(Abschnitt 4). Die bekannte monotone Art bei parabolischen Differentialglei- 
chungen dagegen gehért zur erstgenannten Klasse. 


In Abschnitt 1 werden ferner hinreichende Kriterien fiir die monotone Art 
von Operatoren M bewiesen. Sie sind von der Form: M ist (in bestimmtem 
Sinne) von monotoner Art, wenn ein ,,gréBerer’’ Operator A, d.h. ein Operator 
A=M-+B mit positivem B, (in bestimmtem, eventuell anderem Sinne) von 
monotoner Art ist und ein Element z mit z>o und MzZo (bzw. Mz>o oder 
Mz > 0) existiert. Die Beweise sind elementar und sehr einfach. 


Praktisch hat man also eine monotone Art fiir einen gréBeren und einfacheren 
Operator A direkt nachzuweisen und ein geeignetes z zu konstruieren. Mit dieser 
Methode erhalten wir in den Abschnitten 2 bis 5 und 6.3 Bedingungen fiir die 
monotone Art bei verschiedenen Problemtypen. Bereits bekannte Resultate, zu 
deren Beweis man bisher zum Teil sehr verschiedenartige Mittel verwendete, 
werden so auf einheitliche und einfache Weise hergeleitet. Zum Beispiel kann 
man elliptische und parabolische Differentialgleichungen weitgehend gemeinsam 
behandeln. Dariiber hinaus erhalten wir eine Reihe neuer Ergebnisse. Zum 
Beispiel werden teilweise allgemeinere Problemtypen behandelt, Voraussetzungen 
abgeschwacht, starkere Behauptungen bewiesen. Bei Aufgaben mit parabolischen 
Differentialgleichungen wird gezeigt (wie bereits erwahnt), daB sie auch im Sinne 
der urspriinglichen Definition von L. CoLLatz [4] von monotoner Art sind. 

Bei der Fehlerabschatzung fiir eine Aufgabe (0.1) nach der Methode von 
L. CoLtratz [4] mu8 man ein Element z der gewiinschten Art ohnehin konstruieren. 
Praktisch bedeutet dies: Man braucht nur die monotone Art eines gréBeren 
Operators A, nicht die von M, zu beweisen; ist dann der Formalismus der Ab- 
schatzung durchfiihrbar, sind die erhaltenen Schranken auch giiltig (Satz 1.2). 
Hat man umgekehrt ein geeignetes z, 1aBt sich auch der Formalismus durchfihren. 


1. Ergebnisse ftir abstrakte Operatoren 
1.1. R= {u, v, ...} sei ein reeller halbgeordneter linearer Raum. Dieser Raum 
sei ferner als ,,archimedisch“ vorausgesetzt, d.h.: 
Aus nu Sv fiir alle n=1,2,... folge wo, (1.4) 
wenn u, v feste Elemente €R bedeuten und o das Nullelement ist. 


Wir schreiben: 
u>o, wenn uo und u=0, 
u>o, wenn zu jedem vER eine natiirliche Zahl n mit nu>v existiert, 
uo, wenn u>o oder u=o. 
Offenbar folgt w>o0 aus u>o. 
Es sei z€R ein Element >o und v irgendein Element €R, welches nicht 
<o ist. Dann gibt es eine Zahl « mit 


VEO 


und fiir 
%& =inf{a: —o<a< 0, v<az} 
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gilt Oy = 0. (1.2) 
AuBerdem hat man 


vS (+7), d.h. n(v—%2)S2 (n= 452, 3y)e 


Aus der letzten Ungleichung folgt mit (1.1) v—a%2zS0, also 
% = min{a: —o<a< co, vS a3}. (1.3) 


Es seien nun R und S zwei solche archimedischen reellen halbgeordneten 
linearen Raume. Wir verwenden die gleichen Symbole =, >, >, = fiir die 
Ordnungsrelationen in R und S. Selbst wenn R und S die gleichen Elemente 
haben, kénnen diese Symbole jedoch in R und S etwas Verschiedenes bedeuten. 
Mit A, B, M bezeichnen wir im folgenden lineare Operatoren, die R in S abbilden. 

Wir interessieren uns fiir Operatoren M, welche eine der folgenden Eigen- 
schaften (a), ..., (d), (a’), ..., (d’) besitzen: 


(a) aus Mu>o folgt u>o, (c) aus Mu>o folgt u>o, (1.4) | 
(b) aus Mu>o folgt u>o, (d) aus Mu>o folgt u>o, ; 
(a’) aus Mu=o folgt u=o, (c’) aus Mu=o folgt uZo, 


(b’) aus Mu=o folgt w=o, (d’) aus Mu>o folgt u=o. (1.5) 

L. CotLatz [4] nannte M einen ,,Operator monotoner Art“, falls M die Eigen- 

schaft (c’) besitzt. Wir betrachten hier also nun verschiedene Sorten von ,,Opera- 

toren monotoner Art“. Eine etwa (b) entsprechende Eigenschaft wurde bereits 

bei parabolischen Differentialgleichungen benutzt ([72], [18], [5], [13]). Die» 

folgenden zwei Hilfssatze erlautern die Beziehungen zwischen den verschiedenen 
Eigenschaften in (1.4) und (1.5). 


Hilfssatz 1.1. 1) Jede der Eigenschaften in (1.4) folgt aus der entsprechenden 
Exgenschaft in (1.5). Wenn M™ existiert, sind die Eigenschaften in (1.4) und (4.5) 
paarweise dquivalent. 


2) Wenn M eine der Eigenschaften (1.5) besitzt, existiert M+. 
3) M habe eine der Eigenschaften (1.4), und es gebe ein Element zC R, derart dap 


Mz>o im Falle (a) oder (c), 6 
Mz> 0 im Falle (b) oder (d). (1.00 
Dann existiert M. 

Beweis. 1) ist klar. Zum Beweis von 2) und 3) sei p€R mit My=o. Es 
ist jeweils p=o zu zeigen. 


2) Es gilt auch M(—q)=o. Mit jeder der Eigenschaften (1.5) folgt daher 
sowohl y=o als auch = 20d. Naa: 


3) Die Elemente z+n@p und z—ng (n=1, 2,...) geniigen derselben Unglei- | 
chung wie z in (1.6). In jedem der Falle folgt daraus 


z+ngo>o und z—ngp>o (wad, 2a 
a. hi 


n(—qg)<z und np<z (tA) ee Te 
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Da R archimedisch ist, erhalt man hiermit wiederum —gSXo und ySo, d.h. 
y=o. 

In praktischer Hinsicht ist es keine wesentliche Einschrankung, zu fordern, 
daB ein Element z mit (1.6) existiert. Dies bedeutet nur, da8 es tiberhaupt einen 
Sinn hat, mit der entsprechenden Eigenschaft in (1.4) zu arbeiten. Wir werden 
im folgenden bei der theoretischen Untersuchung hauptsdchlich die Falle (1.4) 
betrachten. 

Mit den verschiedenen Ordnungsdefinitionen kann man auch verschiedene 
Arten der Positivitat von linearen Operatoren C definieren, welche eine Teil- 
menge Dec S in R abbilden: 


(a) aus v>0 (vEDe 
(b’’) aus v>o (vEDe 
(c’’) aus v>0 (vEDe 
(d’’) aus v>o (vEDe 


folgt Cv>o, 
folgt Cu>o, 
folgt Cu>o, 
folet Cu>o0 


und entsprechende Eigenschaften mit dem Gleichheitszeichen. 


Falls M*™ existiert, hat M genau dann eine der Eigenschaften (1.4), wenn 
C=M~ die entsprechende Eigenschaft in (1.7) besitzt. 

Der folgende Hilfssatz gibt Beziehungen zwischen den verschiedenen Eigen- 
schaften in (1.4) an. 


Hilfssatz 1.2. 1) Aus (a) folgen (b), (c) und (d). Sowohl (b) als auch (c) haben 
(d) zur Folge. 


2) Wenn es ein 2,€R mit 
Mz>o (1.8) 
gibt, folgt (c) aus (d). 
3) Wenn es ein 2,€ R mit 
Zz, > 0 (1.9) 
gibt, folgt (b) aus (d). 
4) Wenn es sowohl ein z,€R mit (1.8) als auch ein 2,€R mit (1.9) gibt, sind 
die Eigenschaften (b), (c) und (d) adquivalent. 
Beweis. 1) ist klar. 
2) M habe die Eigenschaft (d), und es sei Muw>o. Dann gilt: 
M(nu+%)=nMu+Mz,>0 (4 =A, 2,830 3) 
und daher wegen (d) 
nutzgz>0, dh. na(—u<xy (AG D3 cisie ie 
Daraus folgt (—u)<o, d.h. u>o (Mu>o schlieBt u=o aus). 
3) M habe wiederum die Eigenschaft (d), es sei aber nun Mw>o. Dann ist 
M (nu — %) =nMu—M2g> 0 
fiir ein gentigend groBes m und daher wegen (d) 
nu—z,>0, adh. Wu > — a> 0. 


4) folgt unmittelbar aus 2) und 3). 
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Die Bedingung, daB Elemente z, und z, mit (4.8) bzw. (1.9) existieren, ist 
wiederum keine wesentliche Einschrankung in praktischer Sicht. Der folgende 
Satz gibt eine Antwort auf die Frage: Unter welchen Bedingungen besitzt 
M=A— B eine der Eigenschaften in (1.4) oder (1.5), wenn A eine dieser Eigen- 
schaften hat ? 


Satz 1.1. A und B seien lineare Operatoren, die R in S abbilden, und folgende ' 


Voraussetzungen seren erfullt: 


i) aus u=o folgt Buzo (u€R), (1.10) 
ii) es gibt en zER mt 

z>o und (A—B)zZo. (1.14) 
Dann gilt: 


1) Wenn A die Eigenschaft (a) besitzt, hat A—B die Eigenschaften (a), (b), 
(c) und (d). 
2) Wenn A die Eigenschaft (b) besitzt, hat A — B die Eigenschajten (b) und (qd). 


3) Ist 
(A — B)z>0 (1.12) 


t 


und besitzt A eine der Eigenschaften (a) oder (a’), so hat A —B alle Eigenschaften 


in (1.4) und (1.5) und (A — B)* existiert. 
4) Ist 
(A — B)z>o (1.13) 
und besitzt A eine der Eigenschajften (b), (c), (d), (b’), (c’) oder (d’), so hat A—B 
alle diese Eigenschaften und (A — B)* existiert. 


Beweis. A habe die Eigenschaft (a) oder die Eigenschaft (b) (Fall (a) bzw. | 


(b)), und es sei w ein Element €R mit 


= Bie >o0 im Falle (a), (1.14) 
> o im Falle (b). 
Wir wollen in beiden Fallen 
u>o (1.15) 
beweisen und nehmen an, diese Ungleichung sei nicht richtig. Dann ist v=—w _ 


nicht <o, und man hat (1.2) fiir die durch (1.3) gegebene Konstante %). Aus 
(1.14), und (1.11) folgt 
im Falle (a), 


Aa} =| Besa Beng ds ; 
~< im Falle (b), 


also 
>o im Falle (a), 


A (% 2 — v) ( 
> o im Falle (b). 


In beiden Fallen ergibt sich daraus 


CHS = Oe 
Dies hat zur Folge, daB 


W(%oz2—v)>z, dh. v<(a— <2 
n 
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fiir ein geniigend groBes n, im Widerspruch zur Definition von a». Deshalb 
gilt (4.15). 

. Dies beweist, daB A—B die Eigenschaft (a) bzw. (b) hat, wenn A diese 
Eigenschaft besitzt. Die iibrigen Behauptungen folgen mit den Hilfssatzen 1.1 
und 1.2 (auf A — B bzw. A angewendet). 


1.2. Es sei nun eine Gleichung 
(A— B)u=s (1.16) 


gegeben. A und B bilden R in S ab, s sei ein festes Element €S. Fiir diese 
Gleichung sei eine Naherungslésung “, bekannt; 


d[u])=—(A—B)u +s (Clea) 


bezeichnet ihren Defekt. Mit Satz 1 erhalt man fiir dieses Problem eine Fehler- 
abschatzung : 


Satz 1.2. i) Der Operator A sei von der Art (a) oder (b) in (1.4) (Fall (a) 
bzw. (b)). 

ul) Fir den Operator B gelte (1.10). 

ill) Es sec sC(A—B)R, ah. die Gleichung (1.16) besitze eine Lésung u* ER. 

iv) Es gebe Elemente v, wCR, derart daB 


v~<w (1.18) 
und 
(A — B)vsxd[u,|S(A—B)w im Falle (a), 


(A — B)vXd[u]X(4—B)w im Falle (b). (1.19) 


Unter diesen Voraussetzungen gilt: Entweder sind beide Funktionen ugtv und 
Uy+w Losungen von (1.16), oder (1.16) besitzt nur die Léisung u*. Im letzten 
Faille ist 

vs u*¥—uXv, (1.20) 


wobei links bzw. rechts das Gleichheitszeichen genau dann richtig ist, wenn im (1.19) 
das Gleichheitszeichen auf der entsprechenden Seite gilt. 


Beweis. Es sind genau dann beide Funktionen u)+-v und u)-+w Lésungen 
von (1.16), wenn in (1.19) auf beiden Seiten das Gleichheitszeichen gilt. Es 
werde nun angenommen, daB dieser Fall nicht vorliegt. Dann hat man 

(A — B)z>o0 im Falle (a), 

(A — B)z>o im Falle (b) 
mit 

Z=W—VUP oO. 
Nach Satz 1.1 existiert daher (A — B)+, und die Lésung u* ist eindeutig. 
Der Defekt kann in der Form 

dug] = (A — B) (w* — 4) (1.21) 
geschrieben werden. Wir betrachten jeweils die rechte Ungleichung in (1.19). 
Wenn dort das Gleichheitszeichen richtig ist, hat man -— (A — B) (uw+w)+s=o, 
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und die eindeutige Lésung ist daher u*=u,+w, so daB in (1.20) w*—um=w 
gilt. Im anderen Falle (< bzw. <) ergibt Satz 1.1 die Aussage u* —u ~<w, 


wenn man (1.21) beachtet. 
Die linke Ungleichung in (1.19) kann in derselben Weise behandelt werden. 


Bei Anwendung des Satzes 2 konstruiert man also mit v und w auch ein 
Element z=w—v, wie es in Satz 1.1 bendtigt wird. Kennt man andererseits 


ein z€R mit z>o und (A—B)z>o, 
so geniigen die Elemente v=az, w=fz 


der Voraussetzung (1.19), wenn man die Konstanten « und # geeignet wahlt. 


Anwendungen 


Die abstrakten Ergebnisse lassen sich in mannigfacher Weise auf konkrete | 


Probleme anwenden. Die folgenden Abschnitte geben einige Beispiele. In diesen 
Beispielen sind die Elemente der Raume R und S Vektoren mit skalaren 
Komponenten, Funktionen, Vektoren mit Funktionen als Komponenten oder 
Kompositionen solcher GréBen. Die Definition wo ist dabei — falls nichts 
anderes bemerkt wird — immer punkt- und koordinatenweise zu verstehen. 


Dies bezeichnen wir als die iibliche Ordnungsdefinition. Es habe w€R (oder S) - 


z.B. die Form u=<@(t),@> mit einem Vektor p(t) = (p* () pat, a0 und einem 
Vektor w= (w');-1,2,..., m- Die p* (t) seien auf einer Menge G erklarte Funktionen, 
die w’ skalar. Dann bedeutet wo im Sinne der iiblichen Ordnungsdefinition!?: 


o' () 20 Mh S42 (6G), 20, Ga 


Alle vorkommenden Gr6éBen werden als reell vorausgesetzt. 


2. Matrizen 


Die Elemente von R und S seien u-dimensionale Vektoren «= (u') mit (reellen) 
Komponenten u’ (¢=1, 2,...,”). Die linearen Operatoren von R in S und 
ebenso die von S in R kénnen dann mit den x xn-Matrizen A = (a;;), B=(8;;), ... 
identifiziert werden. 


Bei der tiblichen Ordnungsdefinition in R und S bedeutet: 
u>o fir wCR oderS: u'>o (sack Re ee id (2.1) 
Eine » xn-Matrix C hat eine der Eigenschaften in (1.7) bei Do=S, wenn 
C720 (t,4 = 12 eet) 
ist und auBerdem bezichungsweise gilt: 
(a’’) alle c,; sind >0, 
’’) es gibt in jeder Zeile mindestens ein c;;>0, 


( 
(c’’) es gibt in jeder Spalte mindestens ein c;;>0, 
(d’’) es gibt mindestens ein C;;>0. 


1 Sind in einer Ungleichung die Argumente einer Funktion hingeschrieben, so ist 
diese Ungleichung punktweise zu verstehen. Sind die Argumente weggelassen, ist 
die Funktion als Element des betrachteten halbgeordneten Raumes anzusehen ‘und 
die Ungleichung als Ordnungsbeziehung in diesem Raum aufzufassen. 
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Hier sieht man daher leicht unmittelbar, daB die Eigenschaften (b’’), (c’’) und 
(q”) aquivalent sind, falls C1=M existiert. In diesem Falle sind jedoch auch 
die entsprechenden Eigenschaften fiir M in (1.4) und fiir C=M7 in (1.7) aqui- 
valent. Wegen dieser Beziehungen ergibt Satz 1.1: 
Satz 2.1. A und B seien zwet nxn-Matrizen mit folgenden Eigenschaften: 
i) A= (a,;) existiert, 
li) es ist entweder a;,;>0 (i 2,.+.,”) (Fall (a)), 
oder o,;2=0 (t,7=4,2,...,”) (Fall (b)), 
iii) b,,20 (,7=1,2,...,n), (2.2) 
iv) es gibt einen Vektor z mit 


Z>0 
yaa (A —B)z>0 im Faille (a), 
Dap 

(A — B)z> 0 tm Falle (b) 23 

(tm Sinne von (2.1)). 

Dann existiert (A —B)*=(u,;), und es gilt 

4; >0 (t,7=1,2,...,n) 1m Falle (a), (2.4) 
i720 (t,7=1,2,...,) tm Falle (b). (2.5) 


Im Spezialfall einer Diagonalmatrix A und einer nicht zerlegbaren Matrix B 
1aBt sich mehr aussagen, als in Satz 2.1, Fall (b) enthalten ist. 
Satz 2.2. A und B seien nxn-Matrizen mit folgenden Eigenschajften: 
1) A 1st eine Diagonalmatrix mit Diagonalelementen a; ;> 0. 
ii) B ist nicht zerlegbar und geniigt (2.2), 
ill) es gibt einen Vektor z>0 mit 
(A — B)z>o0. (2.6) 
Dann existiert (A — B)1=(;;) und es gilt (2.4). 


Beweis. A und B geniigen den Voraussetzungen i), ) und ii) des Satzes 2.1, 
Fall (b). Mit z beginnend kann man einen Vektor z’>o konstruieren, derart 
daB (A — B) z'>o gilt. (Dies Verfahren beschreibt HOUSEHOLDER in [9], S.319.) 
Satz 2.1 sagt daher aus, daB (A — B)* existiert und (2.5) gilt. (2.5) bedeutet: aus 


(A —B)u>o (2.7) 


folgt u>o. Wir zeigen nun, daB (2.7) sogar w>o zur Folge hat, welche Eigen- 


schaft (2.4) 4quivalent ist. 
Es werde angenommen, daB w>o, aber nicht w>o gilt. Dann existiert eine 


Nummer / mit 1<l<n 
und =o. fire Buia, ha 
Tee COP MTU Gd = A ap ey Pe 


Es gibt mindestens ein Element 0;,,>0 mit 1<pSl, 14+1<5qSn, da B nicht 
zerlegbar ist. Daraus folgt, daB die 1,-te Komponente von (A — B) u im Wider- 
spruch zu (2.7) <0 ist. 
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Beispiel. Es sei M=(m;,,) eine n <m-Matrix mit 

| Gs TU yeaa, 

m,;=\—4 {Hi ees UN tafe, 
| 0 sonst, 


wobei die c; Konstanten bedeuten, derart daB 


7 <6. tur alles ; 
POS ee De ae A ae ‘ sd) DE Aa 8) 
m+1|<c, fiir mindestens ein 7 


M kann in der Form 


geschrieben werden. A und B geniigen den Voraussetzungen 1) und u) des 
Satzes 2.2, und die Voraussetzung iii) ist mit 


i je ISOEE yas 
= = st AB Oae 
Z=— (3), 18 in (7 ) 


ebenfalls erfiillt wegen 


(A— B)z=(2I—B)z+(A 21) 2=2(1 cos Ta) e+ ((6 2) z') 


Deshalb existiert M+ und hat nur Elemente > 0. 


G. Scuurz [17] zeigte mit geometrischen Uberlegungen, daB (2.5) gilt, wenn 4 
eine Diagonalmatrix mit Diagonalelementen > 0 ist, 0;;>0 (1, 4=1,2,...,”) statt 
(2.2) gilt, Voraussetzung iv), Fall (b), von Satz 2.1 erfiillt ist und (A —B) + existiert. 
In [4] leitete L. Cortatz das Ergebnis des Satzes 2.1, Fall (b), her mit der Ein- 
schrankung, da8 A diagonal und (2.3) mit z*=1 (i=1, 2,..., 2) erfiillt sein soll. 
Ferner bewies CoLLatz in [4] den Satz 2.2 mit der Einschrankung, daB (2.6) 
wieder mit z'=1 (1=1, 2,..., ) gefordert und (2.5) statt (2.4) behauptet wird. 
Das Ergebnis des Satzes 2.1, Fall (b), fiir 4 =J findet man bei Ky Fan [7]. Die 
allgemeine Behauptung des Falles (b) lat sich leicht aus diesem Spezialfall her- | 
leiten. Beziehungen zu diesem Satz, Fall (b), haben auch determinantentheoreti- 
sche Arbeiten von A. OstRowskI [/4] und H. SCHNEIDER [15]. L. CoLLatz zeigte 
in [6] (S. 44f.), daB Satz 2.1, Fall (b), fiir den Fall einer Diagonalmatrix A richtig — 
ist, und er bewies ferner Satz 2.2 mit der Einschrankung, daB® (2.5) statt (2.4) © 
behauptet wird. 


Man kann auch allgemeinere Ordnungsdefinitionen benutzen. Es sei z.B.— 
H=(h;,) eine nicht singulare Matrix, und 


n 
“= 0 bedeute >i h,,;w=0 (eek Oss ncee 
=I 


Mit dieser Ordnungsdefinition in R und einer entsprechenden Definition u=o 
K 


in S lassen sich entsprechende Resultate herleiten. Man hat z.B. die Matrix 
KA H™ an Stelie von A zu betrachten. Es sei — ohne Beweis — nur ein einfaches 
Beispiel genannt: 
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Beispiel. Die Matrix F = (f;;) mit 


ey larie = %, 
fig=\4 fiir ¢=7-+1 und +=7 —1, 
QO sonst 


besitzt eine Inverse F4=(g;;), welche die Eigenschaft (a”’) beziiglich der Ord- 
nungsdefinition uzo in R und S mit H=((—1)'*76,,) hat. Das bedeutet 


(— 1)'*7 9, >0 (aA; De cot) 
3. Anfangswertaufgaben bei gewéhnlichen Differentialgleichungen 
erster Ordnung 
Wir betrachten (formale) Differentialoperatoren der Form 


dui ie : 
ee Te 
fiir » Funktionen wu’ (t) im Intervall [0, T]. Diese Funktionen werden zu einem 
Vektor u(t) = (w’ (t)) zusammengefaBt. Die a; ; (¢) seien stetig in [0, T], und es gelte 


(3-1) 


Im folgenden wird u(¢) immer als stetig differenzierbar vorausgesetzt. Der ab- 
strakte Satz 1.4 laBt sich in verschiedener Weise anwenden. Wir beweisen als 
Beispiel: 


G20 fit te) (,7=1,27..5%;, OSTST)- 


Satz 3.1. Unter der Voraussetzung (3.1) gilt beziehungsweise: Aus 


= 0) =0 
L,[w] (t) = (OStST) und u'(0) “in (G24 oF! 2) 
=) 0 
= 0 =0 
folgt . (3.2) 
u' (t) >0 (0O<t<T) und uw(0o)>0 
hes Ors) ete 1) (eed) ob, WM 
w(i)>o0  (StST) 
1 (b) 20W* BOS EST) 
Beweis. Mit C,,[0, 7] bezeichnen wir die Menge der Vektoren p= (¢'(é)) 


mit m auf [0, T] stetigen Komponenten g(t). 
4) Zum Beweis der Behauptungen in den ersten beiden Zeilen von (3.2) 

setzen wir 

R (we, (0; Pie a (0) tstetig Pow (0) 05 5G SA ,12).5: 


und definieren A und B auf R durch? 


, n)}, S= C,, [0, Dy 


n \ 
Die a; ; ud) * 


Au=(#—a,;u), Bu=( 


2 Der Strich an der Summe bedeutet: 
Arch. Rational Mech. Anal., Vol. 8 


j=1 


es wird nur tiber 7 ==7 summiert. 
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Bei der tiblichen Ordnungsdefinition bedeutet 


u>o fir u€R: w(t)>0 (O<tST) und u' (0) > 0 (4, 2, ee 


) 
o>o fir veS: #WM>0° OSf=7)- 
Der Operator A hat die Eigenschaft (b), wie man durch explizite Losung der | 
Gleichung Au=s erkennt. Offenbar geniigt ferner B der Ungleichung (1.10). ' 
Fiir das Element 
z=()ER mit #() =te™' 
und 


1a;;,(t)<N (Ost 7 gs et) (3.3) 


is 


7 
gilt schlieBlich z> 0 und (1.13). Daher besitzt A —B nach Satz 1.1 die Eigen- 
schaften (b) und (c’). Das ist die Behauptung der ersten beiden Zeilen in (3.2). 


2) Zum Beweis der restlichen Behauptungen definieren wir: 
R= {ueC, [0, T]: uw’ (é) stetig (¢ =1,2,...,)}, | 


S ist die Menge der Kompositionen v= <p, w> mit »€C,,[0, T] und einem 
reellen m-dimensionalen Vektor w= (w’), 


Au=<(w —a,,;u'), (u'(ty))>, 
B = Seer ines 
w= {( 2 aH!) 0) 
Bei der tiblichen Ordnungsdefinition in R und S bedeutet 
u>o fir uE€R: uw(t)>0 OSts ft: t= 28 ieee 
iO fr V=<MOores: v'(t) >0 (OStST) und w'>0 (4.25% t. cee 


Wieder erkennt man durch explizite Lésung der Aufgabe Au=s, daB A die 
Eigenschaft (b) hat. Fiir B gilt (1.10), und 


z= (2) mit 2() =e"'—und (3.3) 


erfiillt (1.11) und (1.13). Die Eigenschaften (b) und (c’) von A—B sind in > 
diesem Falle mit den letzten beiden Behauptungen in (3.2) identisch. | 

Fiir den Fall konstanter a;; (i, 7=1, 2,..., 2) folgen die in der zweiten und 
vierten Zeile von (3.2) enthaltenen Ergebnisse aus einem Satz bei R. BELLMAN [3] 
(S.172) tiber die Positivitat der Matrix e4*. 

Mit Satz 3.1 lassen sich auch leicht Aussagen fiir Differentialgleichungen — 
hoherer Ordnung herleiten, indem man diese in der iiblichen Weise auf Diffe-— 
rentialgleichungen erster Ordnung zuriickfiihrt. Zum Beispiel kann man so ein 


Ergebnis von L. CoLratz [4] (S. 372) unter schwacheren Voraussetzungen be- 
weilsen. 


4. Randwertaufgaben bei gewéhnlichen Differentialgleichungen 
zweiter Ordnung 
4.1. Es sei 


L[u] = L[u] &) = — po(t) i — py (t) + po(t) u 
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ein im Intervall 0<¢<1 gegebener Differentialausdruck fiir eine Funktion u(t) 
mit stetigen #;(¢) ({=0, 1 
Pz (t) > 0 (O74); 
Daneben betrachten wir die Randausdriicke? 
Ulu] =au(0) —B%(0), Viu] =pu(t) + x(t). 
a, B, y und 6 seien gegebene Konstanten mit 
SOS "020: 
Alle vorkommenden Funktionen u(¢) werden in diesem Abschnitt 4 als zweimal 


stetig differenzierbar auf [0, 1] vorausgesetzt. 
4.2. Es werden eine stetige Funktion P,(t) und Konstanten a,c gewahlt, 
derart daB 
RWZh(), BO>0 (StS), 
Cc 


a= max (a, 0), =max(y,0), a+tB>o0, c+é>0. 


Wir definieren dann 
Lalu]l=—pou—pyut+ yu, 
Us[u] =au(0)—Bu(o), Vylu]=cu(1) —du(1). 
Fiir den Augenblick sei nun u (¢) eine fest gewahlte (zweimal stetig differenzier- 


bare) Funktion, derart daB 
I4f¥JQ@20 (S#S1) (4.1) 


gilt. In einem Punkte ¢,€ [0,1] nehme x(t) ein relatives Minimum an. 
Ist 0<¢)<1, so kann nicht 
u (ty) <0 (4.2) 


sein, da diese Ungleichung zusammen mit #(f)=0 und u(t,)=0 der Voraus- 
setzung (4.1) widerspricht. 
Gilt neben (4.1) 
U,[u] 20, (4.3) 


so kann (4.2) auch im Falle ¢;=0 nicht richtig sein. (4.2) schlieBt den Fall B=0 
aus. Im Falle B>0, a>0 ergeben (4.2) und (4.3) #(0)<0, so daB dann kleinere 
Werte als (0) in jeder (rechtsseitigen) Umgebung von ¢)=0 existieren. Fiir 
B>0, a=0 kann man nur auf “(0)=0 schlieBen, jedoch ergibt (4.1) dann 
#%(0)<0, und «(¢) kann daher in ¢=0 ebenfalls kein relatives Minimum besitzen. 

Ebenso zeigt man, daB (4.1) und V, [wv] 0 den Fall ¢) =1, w(1) <0 ausschlieBen. 


4.8. Es sei nun zunachst 


TZ ONE Oto pS Oy Upeeoe 0 (4.4) 


und 
Gn Cy 


8 Sind die p;(t) periodisch (mit der Periode 1), lassen sich in ahnlicher Weise 
auch Ergebnisse fiir die Randausdriicke U[u]=u(0)—w(1), V[w]=u(0)—w(1) her- 
leiten. 

29* 


420 JOHANN SCHRODER: 


vorausgesetzt. Wir definieren dann 
R={uEC[0, 1]: w(é) stetig, Ulu] = Vw] =o}, 
S= Clon 
und verwenden in R und S die iibliche Ordnungsdefinition. Dann bedeutet 
se (=F 1) 


) 
u>o fir wER: ,u(0)>0 fir f=0, u(0) > oO” fir p=0; . (4.59 
eres fiir, 6=50;0—e4(1)>0 fur 0=05 
o> o. fir, wes: v(t) >0 (OS?s1). 

A und B seien auf R durch 


Au=Ly[u], Bu=(P(t) —po(t)) u 
erklart, so daB 
(A — B)u=L[u]. : 
Wie die Uberlegungen unter 4.2 zeigen, folgt aus 
Au=o (fiir wER), 


daB wo gilt. Das heiBt der Operator A hat die Eigenschaft (c’). 


4.4. Zu zeigen, daB A auch die Eigenschaft (a) besitzt, erfordert zusatzliche 
Uberlegungen. Es seien g(t) und p(t) zweimal stetig differenzierbare Funktionen, 
derart daB 


Lalg]=0, Uslp]=0, (4.6) 
o(0)<0 fir B>O, (0)<0 fir a>O0, (4.7) | 
Iylyl=0, Valy]=0, 
yi)<0 fir d>0, —y(1)<0 fir y>O, 


Wie wir unter 4.2 bewiesen haben, kann (¢) wegen (4.6) in einem Punkte fy 
mit 0S4,<1 kein relatives Minimum annehmen, fiir das g(t))<0 ist. Mit (4.7) 
folgt daraus aber 


PaO falisape aut fir Ost =A etalls ! 
(<0 BO vice’ HAM ORs t Sal vier Oa 
fiir0< f= 1 wfalls BO, | fir 0O<#S1, falls «=O. 

Ebenso sieht man ein: 
fir Os ts 4, /falls 060, fii = eS 
vd<0] (i) OUNT TOE oie See aise 
fir 0<t< 1, falls 6=0, fir 0S 0-24, falls.) 0: 
Daraus ergibt sich zunachst 
WHY=OvO—-ePOPOH>0 (0StS1). (4.10) 


Die Funktion 
‘)y(s) fir OStSs<1 
y(s)p() fir OSsSt<1 
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ist dann die zum Operator A gehérige Greensche Funktion. A~ existiert also und 
hat die Form 


EERO DS Giees (4.11) 


: Aus (4.8) bis (4.11) folgt, daB jede Funktion w= A~v die Eigenschaften (4.5) 
hat, wenn v(t) 20, =0 (OSS) ist. Das bedeutet aber: A+ besitzt die Eigen- 
schaft (a’’) und A daher die Eigenschaft (a). 

4.5. Da B offensichtlich die Forderung (1.10) erfiillt, kann man mit Satz 1.1 
Aussagen tiber A —B machen, wenn ein z€ R existiert, fiir das (1.41) und eventuell 
(1.12) gilt. Im Falle #)(¢)>0 (0<#<1) kann man P—f, setzen und hat 
A — B=A, erhilt also nichts Neues. Es interessiert jetzt, ob A—B auch von 
monotoner Art ist, wenn #,(¢) auch den Wert 0 oder negative Werte annimmt. 

Unter der Voraussetzung (4.4) und 


ata 0 (4.12) 
erhalt man neue Aussagen z.B. mit Hilfe der Funktion 
a(t) =+ {8 501) +ay[O(1) PW — OW P(t) +ByLO() — OW]+ 
+a 6[Q(1) P(t) — P(t) OJ}, 


wobel ? ; 
Pil) Ea pis) as: Q(t) Sa ARN (4.13) 
_ f pal) 
p2(s) t 
Dt)'= exe , l= ols\ass 
A= By +ay P(t) +a52(1) 
ist und @(¢) eine stetige Funktion mit 
O20, 220. (Os r= 1) (4.14) 


bedeutet. 
Man berechnet 


(A—B)4()=—p.4) PHo0O+P a), Ul4l=O und Vi[z]=0. (4.45) 
Sei nun zunachst o(t)=1 (OS¢S1) und 
po) =0 (0StS1). (4.16) 


Dann ergibt eine elementare Betrachtung, da8B z,>o0 ist im Sinne von (4.5), 
und es gilt ferner [(A —B)z]()>0 (OStS1). Mit Satz1.1 folgt daher, daB 
im Falle (4.4), (4.12), (4.16) 4 —B die Eigenschaft (a) hat. Bedeutet also @ (:) 
nun wieder irgendeine stetige Funktion mit (4.14), so ist die Lésung der Aufgabe 


—pri—p%=h.()P He, Ulul=0, V[uJ=o0 


+o im Sinne von (4.5). Diese Lésung ist jedoch gerade z,. Im Hinblick auf 
(4.15) erfiillt also z=z, aus (4.13), (4.14) die Bedingungen (1.11) und (1.12), 
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wenn 


p()PHoeW+pWaQ20, #0 (0S¢S1) (4.17) 


gilt. 
Im Falle 
Pil) =0 (Oe 724) es hee ae 
geht (4.17) z.B. iiber in 
b(t) +P t1—NZ0, #O (OS#S1), falls p()=1 (0StS1), 
(P2(0) + holt) —-)sinat=0, +0 (0<t<1), falls 9()=sinat. 


Insgesamt liefert Satz 1.1 das folgende Ergebnis. 


Satz 4.1a. Es gelte (4.4) und es gebe eine (zweimal stetig differenzerbare) 
Funktion z(t), welche >o0 ist im Sinne von (4.5) und die Bedingungen 


L[z\@ 20,. $0 (083121); Ujzj = Viz|=0 (4.18) 
er fiillt. 
Dann folgt aus 
Llu@20 ((Stan5 Clas Vigi=p; 


daB entweder u=0 (OStS1) ist oder u die Eigenschaften (4.5) besvtzt. 

Im Falle (4.12) erfiillt z=z, aus (4.13), (4.14) die genannten V oraussetzungen, 
falls (4.17) galt. 

4.6. Man kann auch Aussagen tiber Funktionen « machen, welche die Rand- 
bedingungen U[u]=V[u]=0 nicht erfiillen. Dabei braucht dann auch (4.4) 
nicht vorausgesetzt zu werden. Wir definieren jetzt: 

TG Olt, 


S ist die Menge der Kompositionen v= ¢@ (¢), €, 7> mit »€C[0, 1] und reellen | 


Zahlen S. UB 
Au=<L,[u), Us [u],V4[u]>, 


Bu = (E(t) — polt)) u, (a — a) w(0), (c — y) w(1)>. 
Es bedeutet nun 
u>-o fir wER: u(é)>0 Of 1)3 
v>o fir v=<p(t),émeS: p)>0 (St#S1), E>0, n>O0. 


Aus 4.2 folgt, daB A die Eigenschaft (c’) besitzt. Eigenschaft (a) 1a8t sich ! 


dagegen nicht allgemein nachweisen. B erfiillt auch hier (1.10). 
Satz 1.2 liefert daher jetzt 


Satz 4.1b. Es gebe eine (zweimal stetig differenzierbare) Funktion 2(t), derart 
dap 


z0@>0  (oS% 


I 


1)3 (4.19) 
1), (4.20) 
Ulz]>0, V[z]>0 (4.21) 


Li) 0. oes 


I 
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gut. Dann folgt aus 


L[uJQ)20 OS#S1), UlujJzo0, ViujJ=o, 
dap 
u()=0 (o<#S<1) 
ast. 


Man kann hier z.B. mit der folgenden Funktion z=z, arbeiten: 
ee “t fy(P(1) —P@) + 6 P(t)} + “2 {a PW) +} 
mit den Bezeichnungen (4.13), (4.14) und Konstanten ¢,, ¢.>0. Man hat 


L{z»] t) = pot) P@ et) + hot) 22 (2) (4.22) 
Ulz.]=e, V[2q] = é. 


Im Falle (4.4) erfiillt z=z, die Bedingungen (4.19) und (4.21), und iiber (4.20) 
gibt (4.22) Auskunft. Aber auch wenn (4.4) nicht gilt, kann z, die gewiinschten 
Eigenschaften haben. Schwierigkeiten macht hier vor allem die Forderung (4.19). 
Bei e,= €,=¢ und o(t) = 1 sind z.B. folgende Bedingungen hinreichend fiir (4.19) : 


A()/=0 (0S¢S1), 
(yee Ay 
CS =I SSE pes? =O. 


L. CoOLLATZ zeigt in [4] und [6] im wesentlichen, daB L[w](¢)=0, U[u]=0, 
V[u|]=0 die Aussage u(t)=0 gestattet, falls )(t) 20, «20, y=O, aber nicht 
bo (t) =9, «=0 und y=0 ist. Die Beweismethode entspricht der des Abschnittes 
4.2. Fiir selbstadjungierte Aufgaben beweisen einerseits ARONZAJN und SMITH [7] 
und andererseits R. BELLMAN [2] solche Ergebnisse, wobei f)(¢) auch negative 
Werte annehmen kann. ARONSZAJN und SMITH arbeiten mit der Methode der 
reproduzierenden Kerne, BELLMAN verwendet Mittel der Variationsrechnung. — 
Auch H. GLorHsTEHN kann in seiner Arbeit [8] negative f(#) zulassen. In seinen 
Ergebnissen kommt eine z(¢) entsprechende Funktion vor. Die Resultate in [8] 
werden mit einer anderen, aber ebenfalls elementaren Methode bewiesen. 


5. Elliptische und parabolische Differentialgleichungen 

5.1. Es sei G ein offenes, beschranktes (nicht notwendig einfach zusammen- 
hangendes) Gebiet des -dimensionalen reellen Euklidischen Raumes E”. Die 
Elemente von E” werden mit x=(x}, x*,..., x”), y,... bezeichnet. J” bedeute 
den Rand von G und G die abgeschlossene Mannigfaltigkeit G=GUTI. I sei die 
Vereinigung zweier punktfremder Bestandteile J; und J;. J, kann leer sein. 
Andernfalls sei J; ein offenes Gebiet einer (nm —1)-dimensionalen Hyperebene 
x" =const. Fiir jedes x)»€ J} sollen dann auBerdem die Punkte 


x= 4%, —t(0,0,0,...,0,1) mit 0<tS e(%) (5.1) 


bei gentigend kleinem ¢(%) >0 in G liegen. 
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Ableitungen einer auf G erklarten Funktion u(x) nach den Koordinaten x’ 
werden durch ein Komma und nachfolgende Indizes gekennzeichnet, z.B. 
Pulox'dd=u ,,. Fir x Cy ist u,,(%) als einseitige Ableitung aufzufassen. 

Wir betrachten den Differentialausdruck 


vp (x) n 
L{u]=— > 4,;% 4, + D0, (x)u;+e(x)u fir xeGul, 
El 


i,j=1 


mit auf GUJj definierten gleichmaBig beschrankten Koeffizientenfunktionen. 
v(x) bedeute eine (eventuell von x abhangige) nattirliche Zahl, derart daB 


<n fiir «EG, 
<n—1 fir «Ek. 


y (x) 


Fiir jedes feste x€GUI} sei die v(x) xv(x)-Matrix U(x) = (a;;(*)) symmetrisch 
und positiv semidefinit. Ist 3 nicht leer, so gelte 


b,()>0 (xe h). (5.2) 
Ferner benutzen wir Randausdriicke der Form 


Ulu] (x) = a(x) w(x) B(x) 24 (2) fir vel] (5.3) 


mit auf J; definierten gleichmaBig beschrankten Funktionen « (x), B(x), derart daB 
B(x)20 (x€f). 


Oul0o (x) fiir x),Cl] bedeute die einseitige Ableitung in Richtung eines Halb- 
strahles 
L=% +04, (0Sa0< &) (5.4) 

mit einem Vektor yo(%9) der Lange 1. Fiir 0So<e(x,) bei geniigend kleinem 
é(%) >0 sollen die Punkte dieses Strahles ganz in G=GUT liegen. 

Wir definieren R, als die Menge der auf G erklarten Funktionen u(x) mit 
folgenden Eigenschaften: 

a) u(x) ist stetig auf G. 

b) Fiir x CGT) existieren die Ableitungen 


u(x) (7 =1,2,...,m) und 4 4.(x) (4,7 =4, 2, ..., »(2)) (5.5) 
und sind stetig. L[u](x) ist gleichmaBig beschriinkt auf GuoT}. 

c) In jedem Punkte x € I, mit B (x)=+ 0 existiert die in (5.3) vorkommende 
Ableitung 0ujda(x). Ulu] (x) ast gleichmapig beschrankt auf I). 

Jede der in diesem Abschnitt 5 vorkommenden Funktionen (x) soll € Ry sein. 

Wir wahlen eine auf GUJ} definierte gleichmaBig beschrankte Funktion C (x) 
und eine auf J; definierte gleichmaBig beschrankte Funktion a(x) mit 


x Cees) (x€ GUD) 


)20 = 
a(x)20, a(x)Za(x), a(x) +B(x)>0 (x€Q), (5.6) 


Lineare Operatoren mit positiver Inversen 425 


und wir definieren dann 
v (x) n 


L,{w] Sa 2 Mat scth ae se Qe 


Us [u] (x) = a(x) u — B(x) w. 


Die hier in 5.1 genannten Voraussetzungen seien im ganzen Abschnitt 5 
erfiillt. 


5.2. Wir beweisen zuniachst: 


Satz 5.la. Es sei u(x) CR, eine Funktion, derart da 
Lalu|(x)>0 (*EGT), (5-7) 


Us{uJ(x)>0 (xe). ss fs8) 
Dann gilt: 
u(x) >0 (x€G). 
Beweis. Die stetige Funktion «(x) nehme in einem Punkte x)€G ihr Mini- 
mum an. Wir setzen 


u(%) SO (5.9) 
voraus und leiten daraus einen Widerspruch her. 
Zunachst sei %€ GUI; angenommen. Dann ist 


H (8) =04 1 fg ty oy aoe falls x€G, (5.10) 
; Mir p41, 2,2... ier falls” 9,61) 
Ferner ist die Matrix 
U(%) = (4 2; (%0))25=1, 9, we, ¥ (Xo) 

positiv semidefinit. Dies hat 

¥ (%o) 

ii %) u ;;(%) 20 (5.44) 
pe 


zur Folge, denn diese Summe ist die Spur der ebenfalls positiv semidefiniten 
Matrix %(%) W(x»). 

Im Falle x%)€G widersprechen (5.9), (5.10) und (5.11) unmittelbar der Voraus- 
setzung (5.7). Im Falle x,€J} ergeben (5.9), (5.10), (5.11), (5.7) und (5.2) die 


Aussage [du/o(— x")] (%) <0, 


im Widerspruch dazu, daB w(x) das Minimum von u(x) sein soll. 

Es ist aber auch x,€ J; ausgeschlossen. Fiir B(%))=0 ist dies wegen (5.8) 
unmittelbar klar. Fiir B(x) >0 sagt (5.8) aus, daB du/dc(x%))<0 ist, so daB 
u(x») nicht das Minimum sein kann. 

5.8. Um von Satz 5.1a ausgehend mit Satz 1.1 weitere Ergebnisse zu erhalten, 
definieren wir: 

Itty; 

S ist die Menge der Kompositionen v=<qy,@> mit auf GuUJ, definierter 
gleichmaBig beschrankter Funktion g(x) und auf J; definierter gleichmaBig be- 
schrankter Funktion  (*). 
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Bei der iiblichen Ordnungsdefinition bedeutet 


u>o fir uCR: u(x) >0 (xEG), (5.12) 
v>o fir v=(p,0)€S: inf{p(x): £€Gul}}>0, inf {a (x): x€IT}>0. 


Ferner seien A und B auf R definiert durch 
Au=<L,[u],U,[u]>, Bu=<(C(x) —e(x)) u, (a(*) —a(%)) 4). 


Beide Operatoren bilden in S ab. 
Satz 5.1a impliziert offenbar, daB der Operator A die Eigenschaft (b) hat. 
Satz 1.1 liefert daher: 


Satz 5.1b. Es gebe eine Funktion z(x)€Rpo (s. (5.5)), derart dap 
z(x)>0 (x€G), inf{L[z](*): «€GUG}>0, inf{U[z] (x): xe T}>0. (5.13) 
Dann folgt aus 
L[u](x)20 (xeGuh), Ulu|()20 (xEh) (5.14) 
fiir u(x)ERy, dap 
u(x) 20 (x€G) (5.15) 
gilt. 
Beispiel mit n=2: 


Gee a) == FS GS Sy ae Ee ee he Meet p 


ww, fir A=1, 0< <1, 
Uluj= (=u, fir 47 == 1, Opopean, 
uU fin 07 Se 


a(x, x*) = 10 — O(a)? + 2(x*)? erfillt (5.43). 

5.4. In den Abschnitten 5.4 und 5.5 geben wir fiir allgemeinere Falle Funk- 
tionen z an, welche (5.13) geniigen. Dieser Abschnitt 5.4 bezieht sich im wesent- 
lichen auf elliptische Differentialgleichungen. Wir fordern 


ieee (5.16) 
und 
inf {a (x): xEI}>0 (5.17) 
und setzen 
v(x)=n (x€G). (5.18) 


A,(x) bezeichne den kleinsten Eigenwert der »xn-Matrix (a; ; (x). Mit (5.17) 
beschranken wir uns auf Randoperatoren erster und dritter Art. 

Es gibt eine abgeschlossene n-dimensionale Hyperkugel K, derart daB GCK. 
vy sel ihr Radius und 7 der Abstand von ihrem Mittelpunkt *=(x'). Man hat 
ee ‘typ mit 


Deiat. (5.19) 
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Wir definieren 


(t) dt 


(4) ==25 (Fe De fs Pee tds (5.20) 


mit einer Konstanten 6 und einer auf 0<7<7, stetigen Funktion o(r), und 
zwar wahlen wir diese GréBen derart, daB 


o()=0 (0Sr<n) (5.21) 
und 


inf {a (x) 6 — B(x): xEL}>0 (5.22) 
Das ist wegen (5.17) immer méglich. 


Es gilt dann z,>0 im Sinne von (5.12). AuBerdem ergibt sich in x€ Jj (mit 

dem zugehérigen 7) 
Ula] (x) = a(x) (x) + B(x) rexp|—f'0(s) ds} c0s (r(x), 0 (2) 
= a(x) by) — B(x) 7%, 
woraus mit (5.22) 
inf {U[z] (x): x€T}>0 

folgt. 

SchlieBlich berechnet man 


La] ]=exp|— Sol at|[ S as(2) + +7(o@) 2 Y aie wl +> byi)| toa. 


a 


Hierin kann man mit (5.19) abschatzen: 


n 


> ai;( a(nyyZia(y, Daays Sa (x) (x€G), 


i j=l 
erhalt also 
Llz](*%)20(x) — (x€G) 
mit 
n f a 
(2) exp|— fo ad] { Sait) +7(00) Al) — |/ DH) realm. 6.2) 
= = 
Satz 5.1b laBt sich daher erganzen durch 
Satz 5.1c. Es gelte (5.16) und (5.17). Ferner gebe es eine Konstante b und 
eine stetige Funktion o(r) mit den Eigenschaften (5.21), (5.22), derart da fiir dre 
durch (5.23), (5.20) definierte Funktion 3 (x) 


inf{#(x): «EG\>0 (5.24) 
gilt. 
Dann erfiillt z=2z, aus (5.20) die Voraussetzungen des Satzes 5.1b. 
Unter den weiter einschrankenden Voraussetzungen 


e()=0  (xEG) (5.25) 
und 
inti (a) ZG) 0 (5.26) 
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kann man bei gegebener Konstanten b mit (5.22) offenbar immer eine Funktion 
o(r) so wahlen, daB (5.21) und (5.24) gelten. Damit folgt 


Satz 5.1d. Es gelte (5.16), (5.17), (5.25) und (5.26), und es bedeute we Ry 
eine Funktion mit der Eigenschaft 


L[w](x)20 (x€G), Ulw](x)20 (x€t). 
Dann gilt - 
u(x) 20 (xEG). 
Die Voraussetzung (5.26) bedeutet, daB L[w] auf G gleichmaBig elliptisch ist. 
Diese Voraussetzung wird gar nicht einmal immer bendtigt, z.B. nicht im Falle 


int { Y.a,,(x): x€0} >0, b(x)=0 (EG, 7=1,2,...,0). 
t=1 


Satz 5.1d entspricht etwa einem Ergebnis, das L. CoLLatz in [4] und [6] 
mit Hilfe des Randmaximumsatzes beweist. (Siehe auch A. G. MEYER [11].) — 
Fiir selbstadjungierte Differentialgleichungen hat man auSerdem wieder Satze 
von ARONZAJN und SmitTH [7] sowie R. BELLMAN [2], die von ahnlicher Art sind 
und mit entsprechenden Methoden bewiesen wurden, wie die in Abschnitt 4 
erwahnten Ergebnisse bei gewoéhnlichen Differentialgleichungen. — H. GLoI- 
STEHN [8] erhalt ebenfalls Ergebnisse fiir elliptische Differentialgleichungen, 
welche den in Abschnitt 4 erwahnten entsprechen. 

5.5. Der folgende Satz bezieht sich auf parabolische Differentialgleichungen. 
Wir schreiben ¢= x”. 


Satz 5.le. Es sez 
y(x)Sn—1 (xEGUK), (5.27) 


und es gebe eine Konstante N>0O mit der Eigenschaft 
inf {b, (x) N + c(x): xeEGuly}>0, (5.28) 
inf {a (x) — B(x) cos (¢, o(x)) N: «€I}> 0. (5.29) 
Dann folgt aus 
ve Llw](2)20 (xeGR), Uluj(jy=Zo0 (xeEN), 
a 


u(x) 20 (x€G) 
gilt. 


Die Richtigkeit der Behauptung ergibt sich daraus, daB z—e™' die Voraus- 
setzungen des Satzes 5.1 e erfiillt. 


Die Forderung (5.29) enthalt die notwendige Bedingung 
inf {a (x): +E R\>0. (5.30) 
Ist diese Bedingung erfiillt und gilt auBerdem 
inf {6,, (x) : x€Gul}>0, (5.34) 


so gibt es in den folgenden zwei Fallen immer eine Konstante N , fiir welche 
(5.28) und (5.29) gelten. 
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1) Im Falle 
c(x) 20 (x€ Gu) 


cos(t,o(%)) <0 falls B(x)+0 (xEQ), (5.32) 


so ist (5.29) fiir alle N>0 erfiillt und (5.28) fiir alle geniigend groBen N>0O. 

Die Bedingung (5.32) sagt aus, daB der Randoperator U keine Ableitung 
langs eines Halbstrahles enthalt, welcher eine positive Komponente in ¢-Richtung 
besitzt. Da der Differentialoperator L[w] nur die erste Ableitung u,, nicht aber 
die zweite Ableitung u,, enthalt, erscheint diese Forderung natiirlich. 

Fiir parabolische Differentialgleichungen bewiesen K.NaAGumo[12] und 
H. WESTPHAL [18] einen Satz, der sich auf den Fall »—1=1 und die erste 
Randwertaufgabe bezieht. L.CoLratz [5] leitete ein entsprechendes Ergebnis 
fiir beliebiges m und allgemeinere Randbedingungen auf der ,,Mantelflache“‘ eines 
zylindrischen Gebietes her. (Siehe auch K. Nicxer [/3] und L. Corzarz [6], 
S. 329f.; in [6] findet man weitere Literaturangaben.) Die Beweismethode dieser 
Satze ist elementar, benutzt jedoch spezielle Eigenschaften parabolischer Diffe- 
rentialgleichungen. Es wird dabei eine monotone Art bewiesen, welche der 
Eigenschaft (b) entspricht. In Satz 5.1e haben wir die bequemere Eigenschaft (c’) 
verwendet, jedoch sagen unsere Ergebnisse aus, daB (c’) und (b) aquivalent sind. 

5.6. Wir skizzieren erganzend noch einige weitere Uberlegungen, ohne auf 
Einzelheiten einzugehen. 

5.61. Wir denken uns Satz 5.1b auf L, und U, an Stelle von L und U ange- 
wendet. Hat man dann ein z mit den gewiinschten Eigenschaften gefunden, so 
weiB man, da der Operator A die Eigenschaft (c’) besitzt. Da jedoch der Fall 


U[ul(x)=0 = (xEN) (5.33) 
eingeschlossen ist, sagt Satz 5.1b dann auch aus: 
Wenn R ={wE Ry: Ulu] (x) =0 (xEh)}, 
S die Menge der auf GUI} gleichmabig beschrankten Funktionen, 
Au=L,{u] und Bu= (C—c)u 
ist, hat A die Eigenschaft (c’). 

Man kann dann Satz 1.1 auf diese Operatoren A und B anwenden, um weitere 
Ergebnisse zu bekommen. Dabei kann man mit anderen Funktionen z arbeiten, 
solchen nimlich, welche die Randbedingungen (5.33) erfiillen, also z.B. Eigen- 
funktionen und dgl. 

Natiirlich kann man die Randbedingungen auch etwa nur auf einem Teil 
von J; fordern. 


5.62. Die Ergebnisse der vorangehenden Nummern sind z.B. von Nutzen bei 
Abschatzungen fiir Randwertaufgaben der Form 


L{u](x)=r(x) (x€Gulj), Ulu} (x) =y(%) (xEf). 
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Bei manchen Randwertaufgaben wird nicht in allen Punkten «¢€J; eine Rand- 
bedingung gegeben sein, z. B. nicht in gewissen Unstetigkeitspunkten (etwa Ecken) 
des Randes. Jedoch werden die Randbedingungen in der Umgebung einer solchen 
Unstetigkeitsstelle x) oft eine Randbedingung der gewiinschten Form in Xp» selbst 
implizieren. Besitzt u(x) in % z.B. eine Tangentialebene und sind alle vor- 
kommenden Funktionen in einer Umgebung von x, geniigend glatt, so werden | 
die Randbedingungen in der Umgebung von %, oft eine lineare Beziehung zwischen 
u(x) und der Stellung der Tangentialebene festlegen, und damit auch eine lineare 
Beziehung zwischen u(x ) und einer Ableitung @u/éo(x9) in einer geeigneten 
Richtung. Entsprechende Differenzierbarkeitseigenschaften sind dann zusatzlich 
von den Funktionen u€ Ry zu fordern. 


5.63. Wir nehmen 
C(x) = clas (*6Gub)se ae) =ae(een) 


an. Dann sagt Satz 5.1b aus, daB A die Eigenschaft (c’) hat, wenn man ein 
geeignetes z finden kann. Um ein solches z zu konstruieren, muBten wir vielfach 
Randoperatoren zweiter Art ausschlieBen. Wir zeigen nun, wie man die Eigen- 
schaft (c’) auch ohne die Voraussetzung a(x) >0 unter einigen zusatzlichen An- 
nahmen direkt nachweisen kann. Diese Uberlegungen haben mehr theoretisches 
Interesse, da man bei der Fehlerabschatzung gemaB Satz 1.2 praktisch ohnehin | 
eine Funktion z ermitteln muB. 
Es sei jetzt 
inf{C(x): xEG}>0. (5.34) 


Dann 1a8t sich sehr einfach beweisen, daB w(x) in GUJ} kein negatives Minimum 
haben kann, wenn 


L4\a|\(e) 220 (x€ GI) 


ist, und daB auch in einem Punkte x)€J{] kein negatives Minimum vorliegen 
kann, falls 
Us [uw] (%) 20 und «(x)>0 


gilt. Die Beweismethode ist ganz ahnlich wie in Abschnitt 5.2. 
Wir zeigen nun unter geeigneten Voraussetzungen: Im Falle 


Liu] (x) 20 (xE€GUT), (5.35) 
U[u] (x) 20 fiirein x€, a(x) =0, (5.36) 
U(X) <0 (5.37) 


kann u(%) nicht das Minimum von u(x) sein. Zum Beweis nehmen wir an, 
u(x) sei das Minimum. Dann folgt mit (5.36) sofort 


U(X) = 0. (5.38) 
Um einen Widerspruch herzuleiten, bendtigen wir neben (5.34) weitere ein- 
schrankende Voraussetzungen. &1, £2)... , —"~1/g mégen die Koordinaten eines 


rechtwinkligen Koordinatensystems mit Ursprung in %) bedeuten. Es gebe Kon- 
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stanten e, d6>0, derart daB die durch 
e&oa—d(6)?2=0, O<oSd mit (€)2?= (G)24 (E)2 +--+ (12 


beschriebene Punktmenge P ganz in G liegt. Wenn J" und w(x) (in einer Um- 
gebung von %») geniigend glatt sind, kann man folgende Schliisse ziehen. 


Wegen (5.34), (5.35) und (5.37) gilt 


v (x) n 
2 Hi (%) 45 — 215 j(%) u,SC(x)uS—3N (x€G) (5.39) 
i,j= j= 


mit einer Konstanten N>0. 


Wenn u(x) als Funktion auf J] in x, differenzierbar ist, verschwinden in x, 
alle Ableitungen von uw in Richtung von J{. AuBerdem gilt aber (5.38), wobei o 
ins Innere von G gezahlt wird. Verhalten sich die Ableitungen w ,(x) (j=1, 2, ..., 7) 
bei Annaherung «> x, stetig, so gilt also bei gentigend klein gewahltem 6 


Daj (x)u | SN («€P) (5.40) 
p= 
und damit 
v(x) 
), 4; (%) 44; S—2N<0 (xEP). (5.41) 


1,j=1 
Wir definieren nun 
w(x) = u(x) + u(x) — u(x), 
wobei 
v(x) =x ((é)?6 — oe?) 
ist mit einer Konstanten x>0. 
6 sei so klein gewahlt, daB (5.40) gilt. Bei hinreichend kleinem x ist 


v (x) 
|= N (xEP), 


also wegen (5.41) 
Dd 4%; 4;,5—-N<0 (x€P). (5.42) 


Ce Peete} 
i,j=1 
Ferner gilt u(x) >u(x9) fiir x€G, da w(x) fir x€G kein negatives Minimum 
annehmen kann. Bei geniigend kleinem x hat man daher w(x) 20 auf dem ganzen 
Rande von P. AuBerdem berechnete man mit (5.38) 


0 dv 
oo (4) = a (Xoo 0. 


w nimmt im Innern von P also negative Werte an und mu daher dort ein 
Minimum besitzen. Dies widerspricht jedoch (5.42). 


5.64. Zum SchluB des Abschnittes 5 sei noch eine gewisse Erganzung des 
Randmaximumsatzes (bzw. Randminimumsatzes) bewiesen. Sie bezieht sich auf 
Differentialgleichungen der allgemeinen Form in Abschnitt 5.1, also sowohl auf 
elliptische als auch auf parabolische Differentialgleichungen. 
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Satz 5.2. Es existiere eine Funktion z(x)€Ro, welche den Bezehungen 
inf {z (x): x€T}>0, inf {L[z] (x): x€GIT}>0 
geniigt. wu bedeute eine Funktion CRo, derart dap 
L[u](x)20 (xEGuI;) 


und 
pe C(x) SO (x€GuI)) 
mut 
p= int {u(y): ye Ti}. 
Dann ist 


u(x)=m (x€G). (5.43) 
Beweis. Fiir die Funktion v(x) =wu(«) —mw hat man 


L{v] (x) = L[u] (x) — cu = Liu] (x) 20 (x€GuI}) 
und 
v(x)20 (Et). 


Satz 5.1b ergibt daher die Aussage v(x) 20 (x€G), welche (5.43) aquivalent ist. 


6. Ergénzungen und Verallgemeinerungen 

6.1. Mit der Verallgemeinerung einer Methode, welche L. CoLLatz in [4] bei 
Differentialgleichungen anwendete, lassen sich aus den Ergebnissen fiir lineare 
Probleme Resultate fiir nichtlineare Aufgaben herleiten. 

Es sei T ein nichtlinearer Operator, der aus einem linearen Raum FR in einen 
linearen Raum S abbildet. Fiir bestimmte w, v lasse sich die Differenz Tu— Tv 
in der Form 

Tu—Tv=M(u, v) (u—v) (6.1) 


darstellen. Dabei bedeute M(u,v) einen (von uw und v abhangigen) linearen 


Operator, welcher eine lineare Teilmenge ROR in S abbildet. R und S mégen 
die in Abschnitt 1 geforderten Eigenschaften besitzen. Ist M(wu, v) fiir bestimmte 
u,v von monotoner Art im Sinne einer der Eigenschaften in (1.4) und (4.5), 
so folgt aus (6.1), daB T im Bereich dieser Elemente in entsprechendem Sinne 
von monotoner Art ist. 

Im einfachsten Falle ist M von w und v unabhangig. Wir nennen T dann 
affin. Dieser Fall liegt z.B. bei ,,linearen‘‘ inhomogenen Randwertaufgaben vor. 
Fiir einen affinen Operator T lat sich mit Satz 1.1 leicht ein Satz 1.2 ent- 
sprechendes Ergebnis herleiten, das die Fehlerabschatzung fiir eine Gleichung 
Lu=r direkt mit dem Beweis der monotonen Art von T verkniipft. 

6.2. In Abschnitt 1 haben wir mit den Ordnungsbeziehungen =, >, > und 
= gearbeitet. Sind die Elemente v€ S (oder auch die w€R) z.B. Kompositionen 
v=(v},...,v°) mit Elementen v‘ aus archimedischen halbgeordneten linearen 
Raumen S;, so kann man Ordnungsrelationen betrachten, welche fiir die ein- 
zelnen Elemente v' verschiedene Beziehungen benutzen. Zum Beispiel wird im 
Falle =:2 durch 

Oe, Ue 
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eine bestimmte Art von ,,Positivitat’‘ des Elementes v=(v!, v?) beschrieben. Es 
lassen sich dann Satz 1.1 und Satz 1.2 entsprechende Ergebnisse gewinnen. So 
kann man vorgehen, um einige von H. GLOISTEHN [8] fiir Differentialgleichungen 
bewiesene Resultate auf abstraktem Wege herzuleiten. Darauf sei hier jedoch 
nicht naher eingegangen. 


6.3. Fiir hyperbolische Differentialgleichungen geben wir hier nur ein spezielles 
Beispiel. Und zwar betrachten wir den Differentialausdruck 


Lha| = Uyy — a(x, y) u, — b(x, y) Uy — (X,Y) u 
im Gebiet 


G={(x,y): OS *S1, OS ySi4}. 
a, b und ¢ seien auf G stetige Funktionen, derart daB 
G(s WV) On D(x. y) 2.0. 0(%, V0 Tur ; (ay) eG. (6.2) 


Alle im folgenden vorkommenden Funktionen w(x, y) seien mitsamt ihren 
Ableitungen u,,, x, und uw, als stetig auf G vorausgesetzt. 


Satz 6.1. Aus 
Ltwl(x,y)20 (0Sx, yS1) 


USO) 00 NO Sry Se 4) 5 nd. UO, )e=0. Oey S&S 1) (6.3) 
folgt 
u(x,y)20, 4u,(%,y)20 und u,(x,y) 20 (Oras eye) i (G4) 


Beweis. Es sei R die Menge der auf G definierten Funktionen u(x, y), welche 
die erwahnten Differenzierbarkeitseigenschaften besitzen und den (charakteristi- 
schen) Anfangsbedingungen (6.3) gentigen. S bedeute den Raum der auf G 
stetigen Funktionen. In S verwenden wir die iibliche Ordnungsdefinition, in R 
dagegen erklaren wir wo durch (6.4). 

Der durch Au=u,, auf R definierte Operator hat dann die Eigenschaft (c’), 
und fiir den durch Bu=a(x, y) U,+ d(x, V)u,+-c(x, y)w gegebenen Operator B 
gilt (4.10) wegen (6.2). Die Funktion z(x, y) =x ye" is > 0 als Element ER 
und erfiillt (1.13) fiir eine geniigend groBe Konstante N. Daher folgt Satz 6.1 
mit Satz 1.1. 
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A Mixed Boundary-Value Problem 
for Hyperbolic First-Order Systems with Derivatives 
in the Boundary Conditions 


VIDAR THOMEE 


Commumicated by W. MAGNUS 


This paper* is concerned with a modification of the customary boundary- 
value problem of mixed type for first-order hyperbolic systems of partial differ- 
ential equations in two independent variables. The modification consists in that 
first-order derivatives, and not only the unknown functions themselves, appear 
in the boundary conditions. Boundary-value problems of this type occur in gas 
dynamics, as was pointed out to the author by GLENN LEwIs. 

This paper is to be considered as an immediate continuation of the previous 
paper [/]. The notation is not re-explained here, and arguments coinciding with 
arguments in [/] are not repeated in detail. 


1. The boundary-value problem 
We consider the linear hyperbolic operator 


ou ou 
(1.1) Lu= aE Q 
in the domain W = {(x, t)| x20, 120, x+#t/AS1} in Euclidean 2-space. Here 
QZ is a real diagonal matrix CG ,.»,(V) of the form 
Ge 
‘ Q= 
3 2) 


where D*C Gor yine(V) and DCECp-yn-(W) are positive and negative in ¥, 
respectively, and YC@p,.,,(V). We assume 


+ Gu, LE&(E(Y), BV) 


0<AS|D|FZ. 
Corresponding to the partitioning (1.2) of Z, we can write (1.1) as 
ou* Gut 
yp + + 
AGO) a Q ae 1+Gty, 
a BROS Be yee = 
EO i) = Pare G pa +Gu 


(w= (ut, u>), Wi COLMA) ne CoP CC ae SF ); G CCnxm-(V)). 


* This paper was written while the author was a Temporary Member at the 
Institute of Mathematical Sciences, New York University. 
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With 0° and 9 being the restriction operators to Pay att=0) and 
LS il anvil we have the initial-value operator 


Su=o'u, Ferx (Cr(V), fod | 


We partition u~ into (wr s u’), uU- oe oe ),2=1, 2. On SF we have the boundary 
operator f-u= (fu, Fu) Eg (Cn €p,(S-)),. where 


fu = 07 aw +Biou, $€E€ 


ot 
fu=Bo-u, £ 
Here 
BEE lS), B= (GB TCO lS), 


where J* is the identity matrix of type mx me. 
The boundary-value problem we consider is the following: Given F, £°, # , 
find w such that 


Sa 
(1.3) fu = f°, 
US. 


2. The discrete boundary-value problem 


In the same way as in [1], we approximate # by means of the difference 
operators 


LPs A,ut ce Dea ln oe Gtu, LAE 2 (€,(V), Ent (Y)) , 
L-u=A,u-—D-Agu+Gu, L-€&(G,(V), G-(Y)) 
(for notation, cf. [7)). 
The discrete initial-value operator is 


Pu=ou, PER (G,(V), G,(S%). 
Corresponding to Y~ we introduce the notation 


S* = {(0, th)| 7 =0,..., 2}, 
S? = {(0, th)|t=1,..., 2}, 
So Vn ore teed, 2 
and the discrete boundary operators 
Fuso Aw+Bou, FE&(G(V),6(S)), 
Pu= Bou, FER (G,(V), @ s(S>)), 
where Bi=r Z’, 7 and ot being the restriction operators from Y~ and V to Sh. 
In addition to the boundary-value problem (1.3), we shall consider the following 
discrete boundary-value problem: Given F=(¥*, F-)C@(V), PCC AF», 
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1 2 


£=(6, £)C@(S)), find we, (V) such that 


ie 
ame ae 

(2.1) Piaf, 
rust, 

Pu f, 


In the same way as in [/] we have the following theorem: 


Where Fe=k, 7", 


where 
where 
where 


where 


F-2 RF 
0— 0 70, 
1 dees: 
pay 
Perf. 
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Theorem 2.1. There is one and only one solution u€&,,(V) of the discrete 


boundary-value problem (2.1). 


The solution of (2.1) can be determined recursively from the following explicit 
scheme, which is obtained directly from the definitions 


==) 0 
Ugg = Pug, 


C= Oy Shows oy TE 


Me 41 = ADE, Mia t (LA De) Ue eo 
= hGe Ue Lae, 


ine) 
ho 
Soe! 


—RGz Ue +RL u, 7, 


ik Bf ah if 
Uo, 1+1 = M0,r —R Bim, +kl 4; 


2 


&€=071,...,%—t—1; 
Up cpr = (E+ AD57) Wg a BD ge gee 
E=1,;....n—t—1;3 


= 2 2+ + gi 1 
U0, 741 = ) as Baa Uo, 741 — Ber U0, r+1- 


3. Estimates 
Theorem 3.1. There is a constant C, independent of n and of ucG,,(V), such 


that 


July SCf|Lt aly, + [lye [Oe 


Proof. Put 
a? = max (|u*|sr, 
ae = | u-|s¢, 


ie — MAG) 


K=|Lt uly, +|L-u|ly-+ [Pulse t [Pulse 


u"|se), 


50+ | 2] 51+ |i? | 52}. 


Then, because of the boundedness of the coefficients in the scheme (2.2) and the 


condition 0<AS< |Z|7, we obtain 


Gn az +n1C, (a,+ Ky) 


G54, max{az+n1C,(a,+ K); Cs(az.1+ K)}, 


The theorem therefore follows from Lemma 4.2 in [/]. 


In this way we can also estimate difference quotients and obtain the follow- 


ing theorem. 
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Theorem 3.2. Assume that DG, GY, B', BEG". Then there 1s a constant Ge 
independent of n and of uC&,,(V), such that 


2 — 
ia C{|L* alr, Ping oe uly; gt |Pulse,g+ IF u| 5 at |Fulst.gt|% leek: 


Just as in [7] we can obtain the following theorem in the limit. 


Theorem 3.3. Assume that G,G, B, B*CE!; then there is a constant C, 
independent of ucGi**(V), such that 


| Hl rg S CLL uly gt [Ol gt [Eula 


4. Compatibility conditions 
We say that FCG" (V), PEGS), pees ae, £ECMS ) satisfy 
compatibility conditions of order g at the origin (for the boundary-value problem 
(1.3)) if there is a wC@/ (VY) such that 


6* (Lu) (0,0) = 8*F(0,0), |a|<Sq—1, 


| d° (ton) (0) =a’ £90), Bq (a = a) ; 

(4.4 il 1 
#0) =f), Bsq-1 (d=), 
ea) =#f0), psq (d=). 


It is of course a necessary condition for the existence of a solution in @ (¥) of the 


boundary-value problem (1.3) that F, £°, £-, & satisfy compatibility conditions 
of order g at the origin. We have the following theorem: : 


Theorem 4.1. There are ei sgais Mi ze and Ni,, which depend only on GZ, GY, 


yur 
B, and B*, such that if F, £°, Jas and ie satisfy compatibility conditions of order | 
q at the origin, then for ySq—1 and yXq im the cases 1=1 and 2, respectively, 


a’ £(0 = 2 M},d°£°(0) + Ni, &F(0, 0), 
(4.2) “ 
2 (0) ~~ Y M2, d®£9(0) + »' N2,eF(0, 0), 
B=0 ja| Sy—1 


and such that if on the other hand (4.2) as satisfied for F, £°, o. and = then we 
can find a polynomial u=1 (x, t) of degree q such that the relations (4.1) are satisfied. 


In particular, for g=1 the matrices M and N are given in the following 
theorem: 


Theorem 4.2. Assume that B*C@1. Necessary conditions that the boundary- 
value problem (1.3) have a solution in €.(V) are the relations 


(43)  £°(0) =B*(0, 0) d#o*(0) — G (0, 0) £9(0) + F*(0, 0) + B(0) £°(0), 
(4.4)  £°(0) = 4(0) £9(0) 


) 
(4.5) d£° (0) = F(0) [ZO, apo F(0, 0) £°(0) + F(0, 0)] + dB? (0) £°(0). 
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In order to see that conditions (4.3), (4.4), and (4.5) are also essentially 
sufficient compatibility conditions for the boundary-value problem (1.3) we prove 
the following lemma: 


Lemma 4.1. Assume that G,G,B, F, ECO: BP’, CC, and that the 
conditions (4.3), (4.4), and (4.5) ave satisfied. Then there is a constant C such 


that if a is the solution of the discrete boundary-value problem (2.1), then, for 
SC; 


Proof. We first remark that A? TRANS k=O, 1,2, are bounded since 
°C¢2, We have - . 4 
4 Ante o— Di Ache. ese Goofot Kivo, 
1 1 
A.M 0=fo — Bo tos 


and thus A, “9 9 and A, Us, 0 are bounded. We have 


2 2 Leet 
ee eee 2+ + Sy — 
U,1 =f; — B; U0,, — B? Uo,» t=0,1, 2. 


In case t=0, this follows from (4.4); in case t=1, 2, from (2.2). Therefore 
A, o,0 =A, (fz — BY* ,, — BE up, -)y 
is bounded because of the regularity assumptions and the boundedness of 4, 19,9 
and Avie 6: already proved. We have 
Az A, Uo,0 = As (Dio Az Ui,0 — Gé,0 Ue,0+ Feo)o, 
A? u6,0 = AMPS Agu, — Go,2%,2+Forzo: 
A? 5,0 =A, (f;— Bi uy Jo, 
Atiso —=Ai(fi— Bi us,,— BY ,.)o, 
and we may conclude succesively that A; 4,u6.0, 47 “6,0, Ly and A? U6, 
remain bounded. For A; 4, 6,0 and A; A, TN after simple computations, we 
obtain 1 1 z 1 
A, A, 5,0= 1{F*(0,0) df" (0) — $*(0,0) £4(0) + F*(0, 0) + 
+ B (0) £°(0) ose (0)} + O41) when noo, 
A; A, 46,9 = "{B* (0) [B(0, 0) 44°(0) — (0, 0) £°(0) + F(0, 0)] + 
+ 4B*(0) £9(0) — af*(0)+ 
+ B* (0) [f° (0) eee 0) df= (0) + F (0, 0) £°(0) — 
— F*(0, 0) — B(0) £9(0) J} + O(1) when n>oo. 


According to (4.3) and (4.5) this implies the boundedness of | wzlo.2: 1—Ae 2, 
and the proof of Lemma 4.1 is completed. 


Lemma 4.2. I} F, ase 2, £ CC, and if 
F(0,0)=0, |e|S¢—14, 

a (0) =0, BSq, 
70) te eh = eal 

a e(0) =0, pq, 


(4.6) 
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then there is a constant C such that if u, is the solution of the discrete boundary- 
value problem (2.1), then for all n 


(4.7) | “nlo,gt1 = C- 


Proof. According to (2.2) and (4.6) we have the following expression, in , 
which the coefficients F, ,, Q,,R:, and T, remain bounded when noo: | 


ar| |ot| 
sont SRR DOA te SRA DEA 


€4+7S|a|—1 


= O(n'|-2-?) = O11) whenn—+oo for |a| Sq+1; 
thus (4.7) holds. 


5, Existence and uniqueness theorems 


The following uniqueness and existence theorems are proved in the same | 
way as are the corresponding theorems in [/]. | 


Theorem 5.1. If F, #°, £- €@®, then there 1s at most one solution ucEs(V) 
of the boundary-value problem (1.3). 


Theorem 5.2. Assume that DZ, G, B, F, £°, £- €C* and that the conditions (4.3), 
(4.4), and (4.5) are satisfied. Then there is a solution uC G,(V) of the boundary- 
value problem (1.3), and if u,, 1s the solution of the corresponding discrete boundary- 
value problem (2.1), then 

EA*u, > &%u, Ja) <1 


uniformly in ¥ when n—>oo. (E ts the linear interpolation operator.) 


Theorem 5.3. Assume that g=1, that D,G, B, F, £°, Pi Ai Bee and that F, 
f°, f satisfy compatibility conditions of order gq at the origin. Then there is a 
solution uC Ci (V) of the boundary-value problem (1.3). 


We also have the following theorem concerning the rate of convergence: 
Theorem 5.4. Assume that D, G, B, F, £°,£ €C and that the compatibility 


relations of order 2 are satisfied at the origin. If u and u, are the solutions of the 
problems (1.3) and (2.1), then | 


(5.1) |E u, —Uly= O(n) when n>. 


Proof. Using Theorem 3.1, we get (5.1) in the same way as in the proof 
of Theorem 6.5 in [7], if we only add the following estimate (OS0S1): 


ies FRu,| = | (ck) -- FRu,| 


a i 
= |S (0, th) + B(x h) w(0, cA) — A, (0, t&) — BU(ck) w(0, A) 


1 
Cu 


k 
oF (0, («4 8) R)| <0 [ya = O(n) when nu>oo. 


“B 
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6. The quasi-linear case 
In this section we shall give local existence and uniqueness theorems for the 
quasi-linear boundary-value problem 


49" _F in ¥, G=Dutuw, FF (2,10, 
w=~0 on SF, f= $%(x), 
aE: OD P= Gi. 
Wee 7 OURS a; Yeo ur, wu). 
Vv, SF, and S-, with 4 to be specified shortly, have the same meaning as before. 
DDOUL OD, F, 7°, oA Vi we make the following assumptions: 
1) £°CGz(S). There are closed sets 
Q= {(x,t,uw)|OS*S1; 0OSt56,; |w—£2(x)| Se}, 
Y= {(t,w)|OStS 4; |u—£9(0)| Se}, 
O = {(t, ut, w)|0 StS 4; |u* —£*(0)| Se: 


— £%*(0)| Se}, 
such that DE C2..m(Q), FEC (Q), £CCA(Q), £°CC%x(Q"). D can be written 


Cpe ie 
“(o 2} 
where D*>0 and D-<0 in Q. 
2) The following compatibility conditions are satisfied at the origin: 


a) £70, (0) = B*(0, 0, £9(0)) *. £0"(0) + F*(0, 0, £9(0)), 
B) A OAH Aad) at Ak 
Y) oF" (0,£9*(0), £*(0)) + 
tact , £°*(0), £°(0)) |F* (0, 0, £°(0)) -F?* (0) + F* (0, 0, £9(0)) | + 
+ on Z°(0, f9*(0). £°*(0)) [B*(0, 0. £°(0)) F-#"* (0) + F*(0, 0, £°(0))| 


d 2 2 
=F (0, 0, £°(0)) $ £%* (0) + F7(0, 0, £°(0)). 
We assume that 4 is a positive number <|Q|5', and we consider the discrete 
boundary-value problem 
A,ui,—Di,4iu,=Fe, on XV, 
A, ue ,-— Dp, Aiue.=—Fe, on Vi, 


(6.2) Wee ae ici) fon 5°. 
7A Ug = f on S, 
2 2 
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Theorem 6.1. There is one and only one solution u€G,,(V) of the discrete 
boundary-value problem (6.2). 

The solution can be defined by induction over t by means of the following 
scheme: 


pg =I, Siz Oetes 5 
Ui 41 = ADE Wey et I —ADE) Meet Rly, §=0,..."-T—1, 
Hig a= (THAD) Vag — Devine ee ee ee 


ak 1 1 
UO, 1+1 = Uo,2+ k Tes 


2 2 
Uo, 7+1 — (ae 


Theorem 6.2. Assume that D, F, £°, ta Be satisfy 1) and 2). Then there 
exists a 6>0 such that the boundary-value problem (6.1) has a solution in ¥°. 
The proof will not be carried through. It is analogous to the corresponding 
one in [/] and consists in proving that there is a 6 > 0 such that, independently 


of n, |u|yi.<C. The modification needed is the same as that in Section 3; 
1.¢., to obtain the estimates, u is split into (w*, uw) and w and not into u* and u, 


as was the case in [J]. 


We have the following uniqueness theorem. The corresponding theorem in 
[1] was formulated a little less strongly. Therefore we give the proof of the 
theorem here and observe that as a special case follows a theorem stronger 
than that in [7]. 


Theorem 6.3. Assume that uC@! is a solution in ¥° of the boundary-value 
problem (6.1) and that 9, FE@(Q), £EO(Q), HCO(Q) where 


Q = {(x,t,v)|0S*S1, OStS6; |v —u(z,d| Se}, 


PF ={(t,rv)|0<t<4, |v —u(0,8)| Se}, 


2 


LO Var v)|O<t< 6, |v? — u*(0,4)| Se; vu u- (0, t)| Se}. 


Assume also that D is a diagonal matrix of the form 


+ 
en DEO 
On Z- 
where D*>0 and D-<0 in Q. Assume that v€@ is also a solution of (6.1) in V°. 
Then u=v in V°. 


Proof. We have u(x, 0) = v(x, 0)=£°(x). We put 


(6.3) = sup 7 
w=v in yl 
nso 
and claim that 6;=6. We shall prove that the assumption 6,< 6 is contradictory. 
Assume that 6,<6, and let 6, be a number with 6,<06,<6 such that |v—u| Se 
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in VY, In ¥®, w=u—v satisfies 


ow 5 ow 
a (x, t) — D(x, t, u) om (x, t) 
= F (x,t, u) — F(x, t,v) + [D(x, t,u) — D(x, t, v)] = (4st), 
(6.4) 3 i 


Fe 0.)=£ tu) — £60), 


w (0, 2) =a (t, u*, uw) — Up (é,0*, v) 


The right-hand sides in (6.4) can be written in the forms Yw, Zw, Bt wt + Bw, 


respectively, with GY, #1, Z?*, BCC since ¢.g.-F (x, t, u) —F (x, t, v) can be 
written 
1 


F (x, t, u) —F(x,t,0) = | F (x,t,0+4 Ow) dd-w. 
0 


Therefore w is a solution of the linear boundary-value problem 


(Gu) = G(x, t, u(x,t), B= (— B*, — B*, T°) 


ow ow fee . 5p 
a Gu) ae Gy (0), My 
ow =0, 
B e 
o 5 -—Few=0, on (F)%, 


Bro w=0, on (F-)*. 
Because of Theorem 5.1, w=0 in Y®, which contradicts the assumption (6.3). 
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Bemerkung 


zur Arbeit CLAUS MULLER & RotF LEIs, 
Zur Theorie der Approximationen vermittels stetiger Operatoren, 
Vol. 7, p. 160—167, 1961 


In unserer kiirzlich erschienenen Arbeit ,,Zur Theorie der Approximationen 
vermittels stetiger Operatoren“ [1] wurde die Differenzierbarkeit von Operatoren 
aus Annahmen iiber die Differenzenquotienten geschlossen. Die dort gegebenen 
allgemeinen Formulierungen miissen noch verscharft werden, wenn die ange- 
kiindigten Anwendungen [/, 2] auf die Potentialtheorie gefolgert werden sollen. 

Der Grundgedanke der benutzten SchluBweise besteht darin, aus der Stetig- 
keit und Kompaktheit der Differenzenquotienten 


R(x,h) g == {T(x+h)—Tihe 
die Existenz der Ableitung 
tire ' (Th) — T(0)}g 


zu schlieBen. Dieses Prinzip bendtigt die folgenden Voraussetzungen: 
1. T(x) g ist stetig fir OS xa (lim |Z («) g — T(x») g|| =0). 


2. Es gibt eine Folge 4, mit h,—0, so daB in 0S 4,5 %5%,5a R(x, h,)g 
gleichmaBig beziiglich x konvergiert. 
Daraus folgt [7] mit @ (x) = lim R (x, h,)g fir x, xS x, daB T(x.) g— T(x.) g 


= Sols) dx ist. 


Fir die Ziele der Approximationstheorie ist es entscheidend, aus der Be- 
schranktheit des Differenzenquotienten auf die Existenz einer Nullfolge im Sinne 
der Forderung 2 schlieBen zu kénnen. Dies ist beispielsweise der Fall, wenn die 
zugrunde gelegten Raume reflexiv sind und die Konvergenz im schwachen Sinne 
verstanden wird. 

Ist andererseits 7’(x) vollstetig fiir x >0 und besitzt T(x) die Halbgruppen- 
eigenschaft, so folgt auch fiir nichtreflexive Raume aus der Beschranktheit von 


1 
A(h) g= -{T(h) — Te 
die Existenz einer Nullfolge h, derart, daB der Differenzenquotient 
K(x, hy) § = T(x) A(h,) g 


im Intervall 0< %,<*< x, den Voraussetzungen geniigt und T(x) somit fiir alle 
X= x, differenzierbar ist. 
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Diese SchluBweise laBt sich jedoch nicht auf die Differenzierbarkeit im Null- 
punkt tibertragen, wenn 7(0)=TJ ist, so daB die Voraussetzung 2 nur unter 
weiteren Annahmen gefolgert werden kann. Dies ist beispielsweise méglich, wenn 
die Menge A(h)g nicht nur beschrankt, sondern auch kompakt ist. Legen wir 
die Raume L? (1< p< cc) zugrunde, so kann die Kompaktheit aus der Beschrankt- 
heit gefolgert werden. Im Falle nichtreflexiver Raume ist die Beschranktheit 
des Differenzenquotienten jedoch nicht ausreichend fiir die Existenz der Ab- 
leitung, wie das folgende Beispiel fiir den Raum der stetigen Funktionen auf 
dem Rande des Einheitskreises zeigt (C-Norm): 


Die Funktion /(z)=arc sin (i se) bildet das Innere des Einheitskreises so 


auf den Streifen — a <Rew< a Im w >0 ab, daB der Punkt z=1 in oo iiber- 


geht. e/® ist daher fiir alle z aus dem Einheitskreis gleichmaBig beschrankt. 
Bilden wir 

g (2) = (1— 2) (1+ 24) 
so folgt 


g'(2) =F {1422-32 +i)2 14 2} 


g(e’”) ist fiir alle » stetig, insbesondere ist g(1)=0. Wahlen wir nun U(r, g)= 
Re g(z), so ist U(1, g) stetig in g, und es gilt U(1, 0)=0. Fiir festes g erhalten 
wir mit r<14<1 
Uleo)-U@) 2.0U: 
Veal 


2 = =nVU folgt daraus 


U(r, y) —U (1, 9) 
y¥—1 


Wegen (V U)?=|¢’ 


gleichmabig fiir alle gy. Fiir y=0 besitzt jedoch der Differenzenquotient 


U(r, 0) —U(1, 0) 
y¥—1 


V2(1+7?) 1+?) +1-+7 
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= — (1+77?) cos log 


fiir y->1—0 keinen Grenzwert. 

Die in der voraufgegangenen Arbeit benutzten Methoden fiihren daher dann 
zu den gefolgerten Anwendungen auf das Randverhalten der Potentialfunktionen 
bei Benutzung der Raume C, L! und L®, wenn auBer der Beschranktheit der 
Differenzenquotienten auch deren Kompaktheit verlangt wird. 
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